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Abstract 

We review the general features of nonsingular universes (i.e. those that go from 
an era of accelerated collapse to an expanding era without displaying a singularity) as 
well as cyclic universes. We discuss the mechanisms behind the bounce, and analyze 
examples of solutions that implement these mechanisms. Observational consequences 
of such regular cosmologies are also considered, with emphasis in the behavior of the 
perturbations. 
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The world, an entity out of everything, 
was created by neither gods nor men, 
but was, is and will be eternally living fire, 
regularly becoming ignited 
and regularly becoming extinguished. 

Heraclitu^. 



The proof of this assertion - which is still missing - was left by Heraclitus to future generations. 
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1 Introduction 



The standard cosmological model (SCM) (see for instance [T] for an updated review) furnishes 
an accurate description of the evolution of the universe, which spans approximately 13.7 
billion years. The main hypothesis on which the model is based are the following: 

1. Gravity is described by General Relativity. 

2. The universe obeys the Cosmological Principle |281j . As a consequence, all the relevant 
quantities depend only on global Gaussian time. 

3. Above a certain scale, the matter content of the model is described by a continuous 
distribution of matter/energy, which is described by a perfect fluid. 

In spite of its success, the SCM suffers from a series of problems such as the initial singularity, 
the cosmological horizon^he flatness problem, the baryon asymmetry, and the nature of dark 
matter and dark energy □, Although inflation (which for many is currently a part of the 
SCM) partially or totally answers some of these, it does not solve the crucial problem of the 
initial sing ularity [29] B The existence of an initial singularity is disturbing: a singularity can 
be naturally considered as a source of lawlessness |282j . because the spacetime description 
breaks down "there", and physical laws presuppose spacetime. Regardless of the fact that 
several scenarios have been developed to deal with the singularity issue, the breakdown of 
physical laws continues to be a conundrum after almost a hundred years of the discovery 
of the FLRW solution (which inevitably displays a past singularity, or in the words of 
Friedmann |256] , a beginning of the world) . 

In this review, we shall concentrate precisely on the issue of the initial singularity B We 
will see that non-singular universes have been recurrently present in the scientific literature. 
In spite of the fact that the idea of a cosmological bounce is rather old, the first explicit 
solutions for a bouncing geometry were obtained by Novello and Salim [177j . and Melnikov 
and Orlov [H] in the late 70's. It is legitimate to ask why these solutions did not attract 
the attention of the community then. In the beginning of the 80 's, it was clear that the 
SCM was in crisis (due to the problems mentioned above, to which we may add the creation 
of topological defects, and the lack of a process capable of producing the initial spectrum 

^There are even claims that standard cosmology does not predict the value of the present CMBR tem- 
perature |265j . 

•^Some "open questions" may be added to this list, such as why the Weyl tensor is null, and what the 
future of the universe is. 

^Inflation also presents some problems of its own, such as the identification of the inflaton with a definite 
field of some high-energy theory, the functional form of the potential V in terms of the inflaton |284j , and 
the need of particular initial conditions [285j . See also [286j 

^This acronym refers to the authors that presented for the first time the solution of EE that describes 
a universe with zero pressure (Friedmann |256p and nonzero pressure (Lemaitre |257) ). and to those who 
studied its general mathematical properties and took it to its current form (Robertson [258] and Walker 
P59j ). For historical details, see [MS] 

^We shall not analyze the existence of future singularities, such as the so-called sudden future singularities 
[290] or the "Big Rip" [TSg]. 
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of perturbations, necessary for structure formation). On the other hand, at that time the 
singularity theorems were taken as the last word about the existence of a singularity in 
"reasonable" cosmological models. The appearance of the inflationary theory gave an answer 
to some of the issues in a relatively economical way, and opened the door for an explanation 
of the origin of the spectrum of primordial fluctuations. Faced to these developments, and 
taking into account the status of the singularity theorems at that time, the issue of the 
initial singularity was not pressing anymore, and was temporally abandoned in the hope 
that quantum gravity would properly address it. At the end of the 90 's, the discovery of 
the acceleration of the universe brought back to the front the idea that p + 3p could be 
negative, which is precisely one of the conditions needed for a cosmological bounce in GR, 
and contributed to the revival of nonsingular universes. Bouncing models even made it to the 
headlines in the late 90 's and early XXI century, since some models in principle embedded 
in string theory seemed to suggest that a bouncing geometry could also take care of the 
problems solved by inflation. 

Perhaps the main motivation for nonsingular universes is the avoidance of lawlessness, as 
mentioned above 0. Also, since we do not know how to handle infinite quantities, we would 
like to have at our disposal solutions that do not entail divergencies. As be seen in this 
review, this can be achieved at a classical level, and also by quantum modifications. On a 
historical vein, this situation calls for a parallel with the status of the classical theory of the 
electron by the end of the 19th century. The divergence of the field on the world line of the 
electron led to a deep analysis of Maxwell's theory, including the acceptance of a cooperative 
influence of retarded and advanced fields [227] [^, and the related causality issues. However, 
this divergence is milder than that of some solutions of General Relativity, since it can be 
removed by the interaction of the electron with the environment. Clearly, this is not an 
option when the singularity is that of a cosmological model. 

Another motivation for the elimination of the initial singularity is related to the Cauchy 
problem. In the SCM, the structure of spacetime has a natural foliation (if no closed timelike 
curves are present), from which a global Gaussian coordinate system can be constructed, with 
Qoo = 1, Qoi = 0, in such a way that 

ds'^ = dt^ — Qij dx^dx-' . 

The existence of a global coordinate system allows a rigorous setting for the Cauchy prob- 
lem of initial data. However, it is the gravitational field that diverges on a given spatial 
hypersurface t = const, (denoted by S) at the singularity in the SCM. Hence, the Cauchy 
problem cannot be well formulated on such a surface: we cannot pose on E the initial values 
for the field to evolve. 

^It is worth noting that Einstein was weU aware of the problem of singularities in GR [4 , and he made 
several attempts to regularize some solutions of his theory, such as the so-called Einstein- Rosen bridge, in the 
early 30s. Indeed, he wrote "The theory (GR) is based on a separation of the concepts of the gravitational 
field and matter. While this may be a valid approximation for weak fields, it may presumably be quite 
inadequate for very high densities of matter. One may not therefore assume the validity of the equations for 
very high densities and it is just possible that in a unified theory there would be no such singularity" [5^. 

*In fact, it can be said that the problem of the singularity of the classical theory of the electron was 
transcended, if not resolved, by the quantization of the EM field. 
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There are more arguments that suggest that the singularity should be absent in an 
appropriate cosmological model. According to [TT], the second law of thermodynamics is to 
be supplemented with a limit on the entropy of a system of largest linear dimension R and 
proper energy E, given by 

S 2nR 
— < . 

E - he 

Currently this bound is known to be satisfied in several physical systems |250j . It was 
shown in [10] that the bound is violated as the putative singularity is approached in the 
radiation-dominated FLRW model (taking as R the particle horizon size). The restriction 
to FLRW models was lifted in [250] . where it was shown, independently of the spacetime 
model, and under the assumptions that (1) causality and the strong energy condition (SEC, 
see Appendix) hold, (2) for a given energy density, the matter entropy is always bounded 
from above by the radiation entropy, that the existence of a singularity is inconsistent with 
the entropy bound: a violation occurs at time scales of the order of Planck's time [^. 

From the point of view of quantum mechanics, we could ask if it is possible to repeat 
in gravitation what was done to eliminate the singularity in the classical theory of the 
electron. Namely, can the initial singularity be smoothed via quantum theory of gravity? 
The absence of the initial singularity in a quantum setting is to be expected on qualitative 
grounds. There exists only one quantity with dimensions of length that can be constructed 
from Newton's constant G, the velocity of light c, and Planck's constant h (namely Planck's 
length ipi = ^/Gh/cP). This quantity would play in quantum gravity a role analogous to 
that of the energy of the ground state of the hydrogen atom (which is the only quantity with 
dimensions of energy that can be built with fundamental constants) [8|. As in the hydrogen 
atom, ipi would imply some kind of discreteness, and a spectrum bounded from below, hence 
avoiding the singularity 0. Also, since it is generally assumed that ipi sets the scale for the 
quantum gravity effects, geometries in which curvature can become larger than ipf or can 
vary very rapidly on this scale would be highly improbable. 

Yet another argument that suggests that quantum effects may tame a singularity is 
given by the Rayleigh- Jeans spectrum. According to classical physics, the spectral energy 
distribution of radiation in thermal equilibrium diverges like at high frequencies, but when 
quantum corrections are taken into account, this classical singularity is regularized and the 
Planck distribution applies [9]. We may expect that QG effects would regularize the initial 
singularity. 

As a consequence of all these arguments indicating that the initial singularity may be ab- 
sent in realistic descriptions of the universe, many cosmological solutions displaying a bounce 
were examined in the last decades. In fact, the pattern in scientific cosmologies somehow 
parallels that of the cosmogonic myths in diverse civilizations, which can be classified in 
two broad classes. In one of them, the universe emerges in a single instant of creation (as 
in the Jewish-Christian and the Brazilian Carajas cosmogonies |307j ). In the second class, 

^For an updated discussion of the several types of entropy bounds in the literature, see |268j . 
""^"This expectation has received support from the proof that the spectrum of the volume operator in LQG 
is discrete, see for instance |405j . 
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the universe is eternal, consisting of an infinite series of cycles (as in the cosmogonies of the 
Babylonians and Egyptians) |3U8] . 

We have seen that there are reasons to assume that the initial singularity is not a feature 
of our universe. Quite naturally, the idea of a non-singular universe has been extended to 
encompass cyclic cosmologies, which display phases of expansion and contraction. The first 
scientific account of cyclic universes is in the papers of Friedmann [278] , Einstein [7] , Tolman 
[H], and Lemaitre [3] and his Phoenix universe, all published in the 1930's. A long path has 
been trodden since those days up to recent realizations of these ideas (as for instance |91j . 
see Sect J10.2T4l) . We shall see in Ch{TO]that some cyclic models could potentially solve the 
problems of the standard cosmological model, with the interesting addition that they do not 
need to address the issue of the initial conditions. 

Another motivation to consider bouncing universes comes from the recognition that a 
phase of accelerated contraction can solve some of the problems of the SCM in a manner 
similar to infiation. Let us take for instance the fiatness problem (see also Sect JTOj) . Present 
observations imply that the spatial curvature term, if not negligible, is at least non-dominant 
wrt the curvature term: 

r2 = 1^1 < 1 

but during a phase of standard, decelerated expansion, r grows with time. Indeed, if a ~ t^, 
then r ~ t^~^ . So we need an impressive fine-tuning at, say, the GUT scale, to get the 
observed value of r 0. This problem can be solved by introducing an early phase during 
which the value of r, initially of order 1, decreases so much in time that its subsequent growth 
during FLRW evolution keeps it still below 1 today. This can be achieved by [H] power- 
law infiation (a ~ t^, (3 > 1), pole infiation (a ~ (— ^)^, /? < 0,t — 0_), and accelerated 
contraction (0 < /3 < 1, t — 0_) |397] . Thus, an era of accelerated contraction may solve the 
fiatness problem (and the other kinematical issues of the SCM |9T]). This property helps 
in the construction of a scenario for the creation of the initial spectrum of cosmological 
perturbations in non-singular models (see SectJTTl). 

The main goal of this review is to present some of the many non-singular solutions avail- 
able in the literature, exhibit the mechanism by which they avoid the singularity, and discuss 
what observational consequences follow from these solutions and may be taken (hopefully) 
as an unmistakable evidence of a bounce. We shall not pretend to produce an exhaustive 
list, but we intend to include at least an explicit form for the time evolution of a represen- 
tative member of each type of solution 0. The models examined here will be restricted to 
those close or identical to the FLRW geometry Ef]. Although theories other than GR will 
be examined, we shall not consider multidimensional theories (exception made for models 
derived from string theory, see Sect l3.3.2l) or theories with torsion. 

"'^"'^But notice that the flatness problem may actually not be a problem at all if gravity is not described by 
GR, see SectlO 

^^The issue of singularities in cosmology has been previously dealt with in 

^•^Notice however the solutions given in [20]. These are non-singular but do not display the symmetries of 
the observed universe, although they are very useful as checks of general theorems. 
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We shall start in Sect ll.ll by stating a working definition of nonsingular universe, and 
giving a brief account of the criteria that can be used to determine whether a certain model 
is singular or not. It will suffice for our purposes in this review to define a singularity as the 
region where a physical property of the matter source or the curvature "blows up" ^/T]. In 
fact, since we shall be dealing almost exclusively with geometries of the Friedmann ty pe, the 
singularity is always associated to the divergence of some functional of the curvature o 

Let us remark at this point that there are at least two different types of nonsingular 
universes: (a) bouncing universes (in which the scale factor attains a minimum), and (b) 
"eternal universes" , which are past infinity and ever expanding, and exist forever. Class (a) 
includes cyclic universes. The focus of this review are those models in class (a), although we 
shall review a few examples of models in class (b) in SectJH 

1.0.1 Notation, conventions, etc 

Throughout this report, the Einstein's equations (EE) are given by 

1 

2' 

where A is the cosmological constant, and k = SttG/c"^, which we shall set equal to 1, unless 
stated otherwise, while the metric of the FLRW model is 



ds^ = dr - a\t) 



1 — er^ 



where e = —1, 0, +1. The 3-dimensional surface of homogeneity t =constant is orthogonal 
to a fundamental class of observers endowed with a four- velocity vector field = 5^^. In the 
case of a perfect fluid with energy density p and pressure p, EE take the form 

P + 3(p + p)- = 0, (2) 
a 

a 1 , , A 

- = --(p + 3p) + -, (3) 
a 3 

in which A is the cosmological constant, and the dot denotes the derivative w.r.t. cosmolog- 
ical time. These equations admit a first integral given by the so-called Friedmann equation: 

a\ e 



oP=[-] +- (4) 
3 \a J 

The energy-momentum tensor of a theory specified by Lagrangian C is given by 

_ 2 5{^gC) 



where g = det{g^ 



-"^"^But notice that not all types of singularities have large curvature, and diverging curvature is not the 
basic mechanism behind singularity theorems. If we consider the problem of singularities in a broad sense, 
we seem to be "treating a symptom rather than the cause" when addressing exclusively unbounded curvature 
[301]. 
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1.1 Singularities, bounces, and energy conditions 

The issue of the initial singularity of the FRLW solution was debated for a long time, 
since it was not clear if this singular state was an inherent trace of the universe or just 
a consequence of the high degree of symmetry of the model. This question was discussed 
firstly in an analytical manner by Lifshitz and collaborators in [28], where geometries that 
are solutions of EE with a maximum number of allowed functions were analyzed. The results 
wrongly suggested that the singularity was not unavoidable, but a consequence of the special 
symmetries of the FLRW solution 0. 

From a completely different point of view. Hawking, Penrose, Geroch and others devel- 
oped theorems that give global conditions under which time and null geodesies cannot be 
extended beyond a certain (singular) point [282j . The goal in this case was not about proving 
the existence of a region of spacetime in which some functional of the metric is divergent. 
Instead, the issue of the singularity was considered from a wider perspective, characterizing 
a spacetime as a whole, by way of its global properties, such as the abrupt termination of 
some geodesies in the manifold. Let us present a typical example of these theorems [2]: 

Theorem: The following requirements on space-time 9Jl are mutually inconsistent: 

1. There exists a compact spacelike hypersurface (without boundary) i^; 

2. The divergence 9 of the unit normals to is positive at every point of 

3. R^vV^v'^ < for every non-spacelike vector 

4. IDT is geodesically complete in past timelike directions. 

Notice that the link of this theorem with physics comes through condition (3) via EE, yielding 
a statement about the energy-momentum tensor: 

- I > 0, (6) 

called the strong energy condition (SEC), see the Appendix. Notice also that, although not 
explicitly mentioned, this theorem assume the absence of closed timelike curves [282] . With 
hindsight 0, it can be said that the strength of these theorems is the generality of their 
assumptions (at the time they were conceived), while their weakness is that they give little 
information about how the approach to the singularity is described in terms of the dynamics 
of the theory or about the nature of the singularity. In any case, if we assume that the 
universe is nonsingular, a positive attitude regarding the singularity theorems is to consider 
that they show the limits of applicability of "reasonable" hypothesis (such as GR or the 
energy conditions, see the Appendix) [3Ulj. 

^^For a reappraisal of the work in [28], see for instance |283| and references therein. 

^^From a mathematical point of view, a negative energy could also allow for a bounce. We will not examine 
this possibility in the present paper. 
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A local definition of a bounce can also be given, in the GR framework, in terms of the 
so-called Tolman wormhole fTZl IT3] (see below). Both in this case and in that of the above 
mentioned theorems, the non-singular behavior in GR is only possible when the SEC is 
violated. The assumption of such a condition seemed reasonable in the early seventies, but 
several situations have been examined in the literature that may be relevant in some epoch 
of the evolution of the universe, for which SEC is not fulfilled, such as curvature-coupled 
scalar fields and cosmological inflation [121 1237] . 

Next we shall examine in some detail how the singularity can be avoided. In the following, 
we shall use a simple form of the singularity theorems 0. Let us first introduce some 
definitions (following [H]). The covariant derivative of the 4- velocity of the fluid that 
generates the geometry can be decomposed as follows 

= ^0 h^^ + a^,y + + v^Vy, (7) 

where 6 = v^^^ is the expansion, h^^ = g^^ — v^Vu, the trace- free symmetric shear tensor is 
denoted by a^^, and u^i, is the vorticity tensor (see Eqns. (13631) and (I364p ). Defining 5* byEfl 

5-3' 

the Raychaudhuri equation [15], which follows from Eqn.(I7I) can be written as0 

sf + 2{a' - u') - v^^ = -l{p + Sp) + A, (9) 
where = v^v^-^i, = i)^ is the acceleration. 

Theorem[16j: In a universe where 3p > is valid, A < 0, and = cj'' = at all times, 
at any instant when if = |^ > 0, there must have been a time to < ^IH such that S" — as 
t ^ to. A space-time singularity occurs at t = to? in such a way that p and the temperature 
T diverge. 

Several remarks are in order. First, EE were used to obtain Eqn.(l9]). Hence, the conse- 
quences of the theorem are only valid in the realm of GR. Second, the singularity implied 
in the theorem is universal: any past-directed causal curve ends at it with a finite proper 
length, in line with a coherent definition of a cosmological singularity 0. Third, since there 
is no restriction on the symmetries of the geometry, 6 is in principle a function of all the 
coordinates, so that the theorem applies not only to Friedmann-Lemaitre-Robertson- Walker 
(FLRW) models, but also to most of the spatially homogeneous, and to some inhomogeneous 
models (see examples in [20]). Fourth, as we mentioned before, the condition p + 3p > 0, or 
more generally, SEC, is violated even at the classical level, for instance by the massive scalar 



^^This will suffice for our goals, more refined formulations can be found in [20] . 
^^S corresponds to the scale factor a in the case of the FLRW universe. 
^^This equation was independently obtained by A. Komar [274] . 
^"See \2Q, and [25^ for a classification of singularities. 
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field, and also at the quantum level (as in the Casimir effect o). So it would be desirable to 
have singularity theorems founded on more general energy conditions, but this goal has not 
been achieved yet (see PU]). 

Notice that in the general case, acceleration and/or rotation could in principle avoid the 
singularity [20], but high pressure cannot prevent the initial singularity in the FLRW model. 
Rather, it accelerates the collapse. This can be seen as follows. The conservation equations 
T^^;; = give 

^''P,/. + (P + p¥ = 0, 

Since p^i = in the FLRW, there is no acceleration. Furthermore, the pressure contributes 
to the the active gravitational mass p + 3p. Finally, not even a large and positive A can 
prevent the singularity in the context of the theorem [16| . 

As mentioned before, a bounce can also be defined locally. The minimal conditions from 
a local point of view for a bounce to happen in the case of a FLRW universe were analyzed in 
p!2] , where a Tolman wormhole was defined as a universe that undergoes a collapse, attains 
a minimum radius, and subsequently expands. Adopting in what follows the metric Eqn.([l]), 
to have a bounce it is necessary that = 0, and dh > 0. For this to be a true minimum of 
the scale factor (conventionally located at t = 0) there must exists a time t such that a > 
for all t E (— t, 0) U (0, i). From EE in the FLRW universe we get 

/a e 
\ a 

From these, the combinations relevant for the energy conditions (see Sect J1.3p are: 

d^ In a e 

2 



a 

p + 3p = —6 —. 

a 



From these conditions and d;, = 0, and db > it follows that [12] 

3 bounce and e = — 1 ^ NEC violated, E§ 
3 bounce and (e = 0; d^ > 0) ^ NEC violated. 



^-•^In fact, it has been shown in [27] that the Casimir effect associated to a massive scalar field coupled to 
the Ricci scalar in a closed universe can lead to a bounce. 
^^For the energy conditions, see the Appendix. 
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3 bounce and (e = 1; > ^) ^ NEC violated, 
The definition of p and p and d > imply tliat: 

p + p< 2 ^, 
e 

p-p > 2 — , 
a"' 

p - 3j9 < 0. 

It follows that 

3 bounce and e 7^ 1 ^ NEC violated, 
3 bounce ^ SEC violated. 
The case that minimizes the violations of the energy conditions can be stated as 

3 bounce and (e = +1; d^ < a'^) NEC, WEC, DEC satisfied; SEC violated. 

This result may be expected since the curvature term with e = +1 acts like a negative energy 
density in Friedmann's equation. Notice that in this analysis, only Einstein's equations and 
the point-wise energy conditions were used, without assuming any particular equation of 
state. In a certain sense, this is the inverse of the theorem stated earlier, which assumed the 
validity of the SEC B 

The restriction to a FLRW model was lifted in a subsequent paper [13], and the analysis 
in a general case was done following standard techniques taken from the ordinary wormhole 
case [18] . It was found that even in the case of a geometry with no particular symmetries, 
the SEC must be violated if there is to be a bounce in GR. Consequently, one can conclude 
that the singularity theorems that assume that SEC is valid cannot be improved. A highlight 
in these analysis is that only the local geometrical structure of the bounce was needed; no 
assumptions about asymptotic or topology were required, in contrast with the Hawking- 
Penrose singularity theorems [23]. Equally important is the fact that, as mentioned above, 
SEC may not be such a fundamental physical restriction. 

To summarize what was discussed up to now, we can say that there is a "window of 
opportunity" to avoid the initial singularity in FLRW models at a classical level by one or a 
combination of the following assumptions cj: 

1. Violating SEC in the realm of GrH 

2. Working with a new gravitational theory, as for instance those that add scalar degrees 
of freedom to gravity (Brans-Dicke theory being the paradigmatic example of this type, 
see Sect|3]), or by adopting an action built with higher-order invariants (see Sectjl]). 

^^An analysis along the same lines but with a more general parametrization for the scale factor was carried 
out in [271] . 

^''We shall not consider here the existence of closed timelike curves as a possible cause of a nonsingular 
universe. 

^^A complete analysis of the behavior of the energy conditions for different types of singularities has been 
presented in [25] , 
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As will be seen below, other ways to avoid the singularity are: 

1. Changing the way gravity couples to matter (from minimal to non- minimal coupling, 
see for instance the case of the scalar field in Sect 13]); 

2. Using a non-perfect fluid as a source, see SectJHl 

Finally, quantum gravitational effects also give the chance of a bounce (see Sect j9.2p 0. 
1.2 Extrema of a{t) and p{t) 

Let us study the relations imposed by EE between extrema of the scale factor, the energy 
density, and the energy conditions, in the case of one fluid. Let us recall that the sufficient 
conditions to have a bounce are o 0^ = and > 0, where 6 = 3a /a , and the subindex 
b denotes that the quantities are evaluated at the bounce. It follows from Raychaudhuri's 
equation for the FLRW model (Eqn.(l9])) with A = 0, 

+ j = -lip + Sp), (10) 

that at the bounce we must have (p + 3p)|^ < 0, independently of the value of e (as was also 
shown in the previous section). From the conservation equation, 

p = -{p + pW, 

we see that there may be extrema of p when 6e = (as in the case of a putative bounce) 
and/or when pe = —pe- The second derivative of the energy density is given by 

p = -{p + p)d~{p + p)e. (11) 

Let us assume first that 6^ = with pe + Pe ^ 0, which implies that Pe = and 

Pe = -{Pe + Pe)Oe, Oe = ~{Pe + ^Pe) ■ 

The different possibilities, according to the sign of 9^, pe + pe, and p + 3p are displayed in 
the following table: 



Pe + 3pe 


0e 


Pe +Pe 


pe 


pe 






< 


> 


< 


> 


min. 


min. 






> 


< 


max. 






> 


< 


< 


< 


max. 


max. 






> 


> 


min. 











^^A definition of a nonsingular space using the so-called principle of quantum hyperbolicity has been given 
in [M] 

^^We are assuming that a ^ 0. 
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We see that there are two cases that agree with what may be termed "normal matter" (rows 
2 and 4), in the sense that maximum (minimum) compression leads to maximum (minimum) 
energy density. Notice however that the case in row 2 violates the strong energy condition 
(see Appendix) . The other cases are clearly unusual: minimum density with minimum scale 
factor (row 1), and the opposite (that is, maximum density with maximum scale factor, row 
3)1^. Notice that it is the null energy condition p + p > (see Appendix) and not the SEC 
that is violated at these unusual cases. In fact, if the requirement p + p > is not satisfied, 
then the equation of energy conservation for a perfect fluid, 

P = -e{p + v). (12) 

says that compression would entail a decreasing energy density, which is a rather unexpected 
behavior for a fluid [^. Examples of the four behaviors will be found along this review. 

When an EOS p = \p plus the condition p > are imposecj^. we see that the case in 
row 1 is permitted for A < —1, and that in row 2, for A G (—1, —1/3). The case in row 3 is 
not allowed for any A, while that in row for is permitted for A > —1/3. 

Notice that all the extrema in p in Table 11.21 are global, since the other possibility (given 
by Pe + Pe = 0) leads to an inflection point in p, assuming that p = Ap. 



1.3 Appendix: Energy conditions 

We shall give next the general expression of the energy conditions, and also their form for 
the particular case of the energy-momentum tensor given by 

= diag(p, -p, -p, -p). (13) 

• The null energy condition (NEC) states that for any null vector, 

NEC ^ T^^^'k"" > 0. (14) 

In terms of Eq. (fT3|) , 

NEC ^ p + p>0. (15) 

• The weak energy condition (WEC) asserts that 

WEC ^ T^uV^v'' > (16) 
for any timelike vector. In terms of Eqn. (fT3|) . 

p > 0, andp + p > 0. (17) 



^^The former is precisely the behavior that allows for a bounce in loop quantum gravity [340j (see 
Sect. (|9.2[) ). while the latter is what is found in the so-called big-rip [23 . 

^^Fluids that violate the NEC are called phantom or ghost fluids, and have been studied in |333j . 
^°Notice that some models do not satisfy this conditions, see for instance Ean. p84p . 
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The strong energy condition (SEC) is the assertion that, for any timehke vector, 

SEC ^ (t^, - vf^v" > 0. (18) 

In terms of Eqn.( fT3|) . 

p + p>0, and p + 3p > 0. (19) 



Each of these three conditions has an averaged counterpart There is yet another con- 
dition: 

• The dominant energy condition (DEC) says that for any timehke vector 

DEC ^ Tf^^v'^v" > and Tf.^v" is not spacehke. (20) 

The different energy conditions are not independent. The following relations are valid: 

WEC =^ NEC, (21) 

SEC NEC, (22) 

DEC WEC. (23) 

Notice that if NEC is violated then all the other pointwise energy conditions would be 
violated pLOj. 
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2 Higher-order gravitational theories 



Higher-order terms in the action for gravity (such as -R^, R^^R^'^, etc.) typically appear due 
to quantum effects, either in the case of quantized matter in a fixed gravitational background 
[31] , or in the gravitational effective action as corrections from quantum gravity [32] or string 
theory W\ [33] . These terms are expected to be important in situations of high curvature, 
when the scale factor is small E^. The models that are engineered to work in the intermediate 
regime, where quantized matter fields evolve on a given classical geometry (the so-called 
semiclassical approximation) mirror the path taken in the early days of quantum field theory, 
in which quantum matter was in interaction with a classical electromagnetic background 
field. In the case of gravity, it is generally agreed that this approach may be valid for 
distances above ipi, although this statement can only be verified by a complete quantum 
theory of gravitation, not yet available. As we shall see in ChlHl some models go below ipi, 
incorporating effects expected to be present in the complete theory, but for the time being 
the quest of the "correct theory" at this energy level seems far from being settled. 

2.1 Quantized matter on a fixed background 

Let us start by considering the corrections coming from quantum matter in a given back- 
ground. As shown for instance in ^3], in the models based on the semiclassical approxima- 
tion the mean value of the stress-energy tensor T^i, of a set of quantized fields interacting 
with a classical geometry is plagued with infinities. These divergencies can be removed by 
a suitable modification of EE that follows from a renormalization procedure. In order to 
render the mean value of T^^ finite, the cosmological constant A and Einstein's constant k 
are renormalized, and a counterterm of the form 

AL = ^ {a R^ + p R^" R^,) (24) 

must be introduced in the Lagrangianl^. The corrections arise from the ultraviolet behavior 
of the field modes, which only probe the local geometry, hence the appearance of geometric 
quantities. After the elimination of the divergences and with a convenient choice of a and 
j3, EE with < > as a source preserve their form [M] : 

G,, + A('-'=") g^, = -«:(^^") < T^r^ > • (25) 

Note that such renormalization does not affect the conservation of the energy-momentum 
tensor, that is 

< TllU > 0. (26) 

^""^ Since in this case the non-Unear terms are always coupled to one or more scalar fields we shall consider 
it in Sect. [5X11 

^^As opposed to Lagrangians that are negative powers of R, which are currently being considered as 
candidates to explain the acceleration of the universe [73| . 

■^•^The relevance of this type of series development was discussed also by Sakharov [35] . 
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Since the constants introduced by the counterterm are to be determined by experiment, 
instead of fixing their values so as to ehminate the quadratic contribution to EE (as was 
done in [^), we can shift them as a ^ a + t] and /3 — > /3 + 777. The new equations are 

G,. + vix,. + 7 Z,.) + A^'-^") g,, = < Tj^^ >, (27) 



where 



and 



- X^lv = R{R^lu - ^ Rgfiu) + R ^fi-u - Ofiu^R, (28) 

Z^u = R - OR^u -liOR + R^p R''^)g^, + 2i?"'3 R^^p,. (29) 

Cosmological solutions of Eqn. (1271) in the case of the FLRW metric were studied in [36] . For 
a fiat universe, the equations take the form 

P + 3(^^)(p + p) = 0, (31) 
where tc = l/(c/i), and /i = —2ri{'y + 3). For the case of radiation (p = pca'^/a^), we get 

^.,,..p_^M^,^f«_^Al ^/M', ,32) 






3 V a 

where if = d/a, and pc = p(^c)- If we impose the existence of a bounce by the conditions 
af, > 0,ab = 0, and db > 0, it follows from this equation that p > 0. It as also shown in [36] 
tc < 3.33 X 10~^ sec. in order that the theory does not confiict with the three classical tests 
of GR. 

Vacuum solutions of Eqn. fl271) in the FLRW geometry were studied in [38]. Notice that 
taking the trace of Eqn. (1271) in the absence of matter we obtain 

R + hR + aR=0, 

where a = 1/{2t](1 + 7)), h = d[ln(—gy^'^]/dt. This equation is analogous to that of a 
damped harmonic oscillator. Depending on the sign of the parameter a, there may be 
damped oscillations for R around i? = 0, or exponentially decaying or growing solutions 



Corrections coming from one-loop contributions of conformally-invariant matter fields 
on a FLRW background were studied in |39j (see also [40|). They allow for nonsingular 
solutions that are not of the bouncing type since they describe a universe starting from a 
deSitter state. A thorough analysis of this setting was given in [69], where the back- reaction 
problem for conformally invariant free quantum fields in FLRW spacetimes with radiation 
was studied, for both zero [69] and non-zero [7DJ curvature and/or A. It was found that 
depending on the values of the regularization parameters, there are some bouncing solutions 
that approach FLRW at late times. 
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2.2 Lagrangians depending on the Ricci scalar 

On approaching the singularity, powers of the curvature may be expected to play an impor- 
tant dynamical role, hence other possible nonlinear Lagrangians are those belonging to the 
class defined by 

S = [ V^fiR)d\ (33) 



where f{R) is an arbitrary function of the curvature scalar, encompassing polynomials as 
a particular case. The problem of the singularity using this type of Lagrangians has been 
repeatedly discussed in the literature (see for instance [HI EI])- The EOM that follows from 
this action is 

/' R,, -\fg^.- Dfg^, + f\^.^^ = 0, (34) 
where /' = df /dR. This equation can be expressed in / and its derivatives as 

/' R^^u-^f g,u + /" (i?,M;. -^Rg^u) + /"' {R,^.R,u - R,xR'' g,.) = 0, (35) 

or, using the trace, 

/' (^R^, - ^Rg^^u^ + f" (^R,^;u - I g^^uOR^ + /"' (^R,^.R,. - \RxR'^ = 0. (36) 
The particular example given by 

f{R) = R + aR^ (37) 

was studied by many authors |l6l HTl HHJ H^. In principle a term Rap-ysR"'^'^^ should be 
included in the action, but the existence of a topological invariant yields 

5 j {RaMsR''^^^ - ^RapR""^ + R^) d'^X = 0, 

in such a way that the Riemann- squared term can be omitted. The equations of motion for 
the Lagrangian introduced in Eqn. (l37|l in the presence of matter are 

(1 + 2aR) R^, -\{R + «^') 9^.u + 2«(i?,^;, - URg^,) = -T^,. (38) 



If we restrict ultra-relativistic matter, the — component of this equation yields 

p = 1 + 51 _ 2a.» (() + ? + 1^ + + 2a« fi, (39) 

3 V 3 / a* 



where R = 29 + 4^^/3 + 6e/a^. At the point where the bounce occurs, 9^ = and 9b > 0, 
and Eqn.( l39il reduces to 

• 2 3e / 6ae 4a\ , , , 

Pb = - 2a9t + - 1 + — + ^ • 40 
al \ al 3 ; 
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Let us take as an example the case in which the section is Euchdean. If we want to have a 
minimum with positive energy density, it follows from Eqn. fl40l) that a < 0. As shown in [47j, 
such a choice for the action of the gravitational field admits solutions in the FLRW framework 
that allow a regular transition from a contracting to an expanding phase. Although negative 
values of a remove the initial singularity, it was shown in [171 HH] that the solutions with 
a < do not go to the corresponding FLRW solution (a oc t^/^) for large t. 
A theory that generalizes that defined by Eqn. fl?r|) . namely 

f{R) = R^aK' 

was studied in [H]. It was found that the FLRW solution for n = 4/3 and p = p/3 is regular 
for all values of t, and tends to the radiation solution for large values of t. Later, solutions 
of this theory with dust as a source were found to have similar properties in [42] . 
Another type of corrections, given by the Lagrangian 

£ = i? + A + 5i?2 + Ci?2ln|i?|, (41) 

were studied in [5D] (with B and C constants). The quadratic and logarithmic terms are 
consequences of vacuum polarization |^52j. Although this form of the Lagrangian does not 
eliminate the singularity in the FLRW solutions, addition of particle creation effects through 
a viscosity term does (see ChjS]) 0. 

The stability analysis of the FLRW solution in theories with L = f{R) was performed in 
[51] , along with necessary and sufficient conditions for the existence of singularities. Eqn.f l35l) 
in the case of a FLRW geometry in the presence of matter reduces to [57j 

/"d(a2 a +aad - 2a^ - 2de) + -fa^d + —fa^ + — a^Too = 0. (42) 

6 36 18 

The argument of the function / is given by 

R= ^(^ad + a^ + e). (43) 
ct 

Assuming that near the bounce the scale factor can be developed in a power series as 

a{t) = ao + ^ aif + 1 a2t^ + (44) 
2 



a necessary condition for the bounce was given 

/oao + 6ai/^<0, (45) 

where /o = f{Ro), and Rq = — 6aQ^(aoai + e), and it was assumed that Tqo > 0. In the 
quadratic case given by Eqn. fl37l) . this condition takes the form 

/ \ 

Qae — a^e — Qaa^a^ < 0. (46) 



'''*A Bianchi I solution of this theory with and without self-consistent particle production was considered 
in [53] . It was shown that particle production quickly isotropizes the model. 
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When e = 0, the condition a > is regained, but there are other possibihties when e = 
1,-1 [5T]. In the same vein, but without using a series development, conditions for a 
bounce in f{R) theories were studied in [380j The basic equations are, that follow from 
Raychaudhuri's equation and the Gauss-Codazzi equation are 

(ib _ _Pb fb 
\ab alJ 

These equations were used in [380J to analyze a possible bounce in the theories given by 
= BJ", f2{R) = R + aR"", fsiR) = exp{XR). Bounces for e = ±1 are possible in the 
case of /i. This case can describe an "almost-FRLW" phase folowed by an accelerated phase 
if n > 1 and n is odd for e = — 1 and R> 0. The same happens with n even and n < with 
-R>OorO<n<l with R < 0, where in the second case n can be only rational. For /2, 
closed bounces are allowed for every integer value of m (often together with open bounces). 
For m rational, closed bounces are not allowed in general for < m < 1. For m rational 
with even denominator there is no closed bounce for (m > 1, a < 0) and no bounce at all for 
negative m and a. In the case of /s, one of the following two conditions must be satisfied in 
order to have a bounce: A > and Rb > lia{2pb)/X, or A < and Rb < ln(2p;,)/A. 

Some exact solutions have been recently found in |275] for the theory defined by f{R) = 
R^-^^. For the vacuum case with e = 0, there is bouncing (entirely due to the dynamics of 
the theory), for < 5 < 1/4. There are vacuum solutions for 6 = 1/2 and e 7^ 0, are given 
by 

\1 — er^ ) 

This solution exhibits a bounce for k > 0. Bouncing solutions were also obtained for a 
perfect fluid with p = (7 — l)p in the case 5 = 1/(87 ~ 1) S 

We would like to close this section by pointing out that Eqn|l2] illustrates the fact that 
the flatness problem is not a priori a problem in theories other than GR (no definite behavior 
of |fi — 1| with time follows from H2|) . 



2.2.1 Saturation 



An interesting idea was proposed in [57j to limit the curvature by adding terms in the 
Lagrangian, following the lines that Born and Infeld 0^ devised to avoid singularities in 
electromagnetism. The Born-Infeld Lagrangian, given by 



BI 



^2 _ ^2 



1 - 



(3^ 



(47) 



■^^ Bounce solutions were also shown to exist in orthogonal spatially homogeneous Bianchi cosmologies in 
f{R) = i?" in [38T] . 

■^^Cyclic solutions were obtained in the case 5 = (87 — 4)/(2(7 — 87)) for a convenient choice of the 
integration constants. 
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is such that the invariant — S''^ cannot take values higher than The fact that it takes 
more and more energy to increment the field when it takes values near is a phenomenon 
called saturation 0. A similar cutoff may be postulated for the curvature tensor when 
quantum gravitational fluctuations become non- negligible, that is (presumably), when 

In [57j, non-polynomial Lagrangians f{R) were considered such that they reduce to R when 
R << ipf, and required that f{R) — >■ constant for — > cxd. This condition is of course 
not enough to determine the Lagrangian, but a quahtative guess can be made. A typical 
Lagrangian that fulflUs the above given conditions is 

An approximate solution of the EOM fl42p for fl48p by a development as a power series of t 
for e = was built in |221] . the solution being non-singular though strongly dependent on 
the non-linearities of the chosen Lagrangian. 

The idea of saturation was subsequently explored in [30j, where an explicit nonsingular 
solution given by 

^4^2x1/4 



a(t) = a[l + ^j , (49) 

was inserted in Eqn. fl42p . where cr is a small parameter. This expression tends to the 
radiation-dominated scale factor for (3H'^/a^ » 1. With this a{t) and using that R = 
— 3/3^(j^/a*, Eqn. fl42p can be rewritten as an ordinary linear second-order differential equa- 
tion for /(-R). This equation was integrated for all the values of the 3-curvature. The 
dependence of the resulting f{R) on the chosen form of a(t) was tested in the case e = 
with that obtained from a^{t) = 1 + 2(1 + a)t^ + 1^, which has the same asymptotic limit of 
Eqn. fl49p . The result in this second case is not distinguishable from the flrst. 

A related analysis was carried out in [273], where it was asked that the theory deflned 
by f{R) be asymptotically free (implying that gravity becomes weak at short distances, 
in such a way that pressure may counteract the gravitational attraction, thus avoiding the 
singularity), and also ghost-free (so that the bounce is not caused by negative-energy-density 
matter) The actions studied in [273] that satisfy these requirements were specifled by [fj 

oo 

f{R) = R + J2 CnRD'^R, (50) 

n=0 

and can be rewritten in terms of a higher-derivative scalar-tensor action: 
S = j d'^x^/^ {^R + ^ qDV - {ip{^ - 1) - co^^ 



■^^This is analogous to the fact that it takes an infinite amount of energy to accelerate a mass moving with 
V ~ c'm special relativity. 

■^^For the relation between f{R) theories and ghosts, see |288| . 

■^^It was shown in [273] that polynomial actions in R do not satisfy these requirements. 
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from which it follows that ip = R (from the EOM of $). After a conformal transformation 
and linearization it follows that the EOM for the scalar fields are [273] 

/ oo N 

^ = 300, = 2 I^qDV + co^ 

with $ = e*^. From these we get 

oo 

-6 5^Qn^+M0 = r(n) 



1 



and the scalar propagator is 

r(-p ) 

It is precisely the function V that controls the absence of ghosts and the asymptotic properties 
of the theory, which was parameterized in [273J as r(— p^) = e^^~^ \ with 7 analytic. To 
actually show the existence of bouncing solutions with the properties mentioned above, the 
scale factor 

a{t) = ao cosh ^ 

was imposed in the equation for Cqo written in terms of F and its derivatives, and compared 
with the r.h.s. composed of radiation and cosmological constant, thus yielding the following 
constraints on F: 

2ujV\uj) - 1 > 

(the latter coming from demanding that the bounce be caused by the nonlinearities, and 
not by the radiation energy density). The authors go on to show that the kinetic operator 
defined by 

7(0;) = kiuj — k2Uj'^ + k^uj'^, 

where fcj are constants, satisfies the constraints and has the correct Newtonian limit. So 
a bouncing solution that is ghost and asymptotically free exists for the theory defined by 
Eqn. flSOj) . although the Lagrangian in the original variable R was not exhibited. 



2.3 The limiting curvature hypothesis (LCH) 

A different proposal to deal with the singularity problem in the higher-order-curvature sce- 
nario is to adopt the limiting curvature hypothesis, introduced by M. Markov pU] as the 
limiting density hypothesis lj. The LHC postulates the existence of a maximum value for 
the curvature, in such a way that 

R^ < Ipf, R^^R^" < ipl W^p.sW'^^''' < ipl 

"^^For boucing solutions that implement this hypothesis through modifications of the EOS, see |376| . 
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etc, and that any geometry must approach a definite nonsingular solution (typically the de 
Sitter solution) when the limiting curvature is reached. This automatically guarantees that 
all curvature invariants are finite |291j . A nonsingular higher order theory was constructed 
in [61j in which every contracting and spatially fiat, isotropic universe avoids the big crunch 
by ending up in a deSitter state enforced by the LCH, for all initial conditions and general 
matter content 0- The action used in [61] was the linear action plus a non-linear term I2 
with the property that 

I2i9^^.) = ^ 9,.u = gj^J^ , (51) 

and enforced that I2 ^ for large curvatures using an auxiliary field (see below). In a 
subsequent paper [52], the method was applied to an isotropic, homogeneous universe, both 
in vacuum and in the presence of matter. The solutions corresponding to e = 1 display 
a deSitter bounce. In the case in which matter is present, it is shown that its coupling 
to gravity is asymptotically free. Later, the model was generalized to include a dilaton 
field [63j, in which case it admits fiat bouncing solutions. The starting point is the dilaton 
gravity action with an added non-linear term (J2) times a Lagrange multiplier ip subject to 
a potential V{ip): 

^ = -^1 d'^^ {^'l (^'^)' + ^ ^"^^2 + Vii^)^ . (52) 



The potential is to be tailored from the EOM and the constraint equations in such a way 
that I2, given by 

goes to zero for large curvatures. Notice that this form of I2 satisfies condition ( ISTl) . so all 
the curvature invariants are automatically bounded. Restricting to an FLRW metric with 
fc = 0, the EOM are 

1 n . 



H = -V'{ij), (54) 
X = -3Hx, (55) 
with X = 0, and a prime denotes derivative wrt ip. An example was given in [63j, where 

= ^TT^- (56) 
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was chosen. This potential yields the dilaton gravity action at low curvatures, enforces that 
I2 go to zero at large curvatures, and enables a bounce. By means of a phase space analysis 
of Eqns. (l53p - (|33]) . it was shown [53] that all the solutions are non-singular, and that some of 
them display a bounce either with or without the dilaton. In particular, the fiat bouncing 
solutions with a non-zero dilaton interpolate between a contracting dilaton-dominated phase 



'^Note that the LFH furnishes in this case a nonsingular universe without bounce. 
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and an expanding FLRW epoch, thus avoiding the graceful exit problem of pre-big-bang 
cosmology (see below). 

One obvious drawback of the LCH is that the non-linear terms are not dictated by first 
principles: they are chosen in such a way as to render the theory finite. 



2.4 Appendix: f{R) and scalar-tensor theories 

Higher-order Lagrangians can be related to scalar-tensor gravity (see for instance ^58j). Let 
us start with the function f{R) is given by 

f{R) = R + aR\ (57) 

The EOM that follow from this Lagrangian is 



2aR,^, - (1 - 2aR)R^, + g^, Qai?^ + ii? - 2aUR^ = 0, 



(58) 



the trace of this equation being 

□i?-^ = 0. (59) 
It was shown in |5H] that this theory is equivalent to the one given by the action 

S = j ^/^d^x [{l + 2a^)R-a^'^]. (60) 

Varying independently g^y and (p in the action given in Eqn. fl60l) . one obtains 

1 a 

(1 + 2a) {R^y - -Rgf,u) + - V^g^,u - '2a{^,f,;u - 9t,u) = 0, (61) 

and 

2a{R-ip) = 0. (62) 
In turn, as shown in [65] the conformal transformation 

9fj,,y = i^ + 2a0) gf,u, (63) 

takes this theory to Einstein gravity with a massive scalar field. 

Except in the case in which a vanishes (which is precisely the case in general relativity) the 
second equation yields that the scalar field is nothing but the scalar of curvature. Inserting 
this result into Eqn. (l6ip one arrives precisely at Eqn. fISS]) . The equivalence can be generalized 
to functions f{R) (see [66]) @. Based on the equivalence, the singularity problem in fourth 
order theories was analyzed in [59j for homogeneous cosmological models with a diagonal 
metric. 



^^It was later proved that all higher order, scalar-tensor and string actions are conformally equivalent to 
general relativity with additional scalar fields which have particular (different in each case) self-interaction 
potentials [324] . 
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3 Theories with a scalar field 



3.1 Scalar field in the presence of a potential 

Violations to some of the energy conditions are produced even at the classical level by some 
scalar field theories. From the singularity theorems discussed in Ch. 1, we can expect the 
existence of bouncing solutions in this scenario @. We shall see next examples of avoidance 
of the singularity in scalar field models that violate some of the energy conditions, as well 
as theories with nonminimal coupling. 

A universe filled with radiation and pressureless matter coupled to a classical conformal 
massless scalar field was studied in [76]. The coupling was provided by the action 

^ = -\j (^."^'" + ~ j^^^ ^^^^ 

where n is the mass of the particle, and 

f^i) [ i-gy^/^^ix" - a;^(r))rfr 



-/ j ijdr = -fj d'x 

(this interaction was suggested in f76] as a classical analog of the pion-nucleon coupling). 
Assuming that we have a FLRW universe filled with a uniform distribution of identical fi 
particles, in the continuum approximation, the field equation for ip takes the form 

+ kF = -fN, (65) 

where F = aip and = na^ =constant. The calculation of the trace of the total stress-energy 
tensor from Eg. (1641) yields 

T^" = -/in, 

so we get for the trace of EE 

Att 

a" + ea = —N^. (66) 

3 

Finally the Friedmann equation is given by 

a'2 + = — {F'^ + eF^ + 2Na^ + 2NfF + 2B), (67) 
o 



where i? is a constant that gives the amount of radiation. The system composed of Eqns. m5l 
1671) was solved in [76] for all values of e, and it was shown that a bounce is possible for the 
three cases when some relations between the integration constants are fulfilled. However, 
physical requirements show that only the e = +1 solution can bounce provided N'^f'^ > 2B. 
A nice feature of this solution is that it satisfies the weak energy condition. 
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The role of scalar fields in Cosmology has been examined for instance in |289j . 
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Another non-singular universe based on a scalar field was presented in [77]. A closed 
FLRW model was considered, with a conformally coupled scalar field as matter content, 
which can be thought as a perfect fiuid with comoving velocity defined by 



In this case, the energy density and the pressure are given by 
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7-77 P = 0' 



with p = (7 — l)p. From these equations we get 



a 
a 



2 ' 



d2 + l 



0. 



Introducing the conformal time through dt = a{ri)dri, and with the changes of variables 
u = a' /a, and u = w'/ (cw), with c = 87/2 — 1, the solution for 0(77) is [77| 

a{v) = ao[cos{ci] + d)]^^", 

where and d are integration constants, which were fixed resorting to the limiting curvature 
hypothesis (see Sect. 12.31) . The potential V was then reconstructed in terms of the scale 
factor (assuming that the EOS changes inthe different eras of the universe) and from 
7 = 1 + j)/p, and the evolution of (p was obtained by numerical integration. 
More general models, given by solutions of the theory 

S = [ d^x^/^{F{(l))R - (9„0(9"0 - 2\/(0)}, 



in which (f) is nonminimalluy coupled to gravity through F, were studied in |108] . where it 
was shown that there are bouncing solutions, which were later proved to be unstable |109j . 
A phase-space analysis of the models given by -F(0) = showed the existence of bouncing 
solutions, under certain restrictions on the constants of the potential V^(0) = otcf^ + /50^ + A 

m- 
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Nonsingular solutions for a scalar field in the presence of a potential were also studied in 
for theories defined by 

where u is determined by dU/du = ^cj)"^. The existence of a bounce was shown for a tailored 
potential given by 

Uiu) = A (u-' + i^a;"/(i-") - - 
\ a 

where A is a constant with dimensions of energy density, and a is a number parameterising 
the classes of theories 0- The bounce exists for a < 1/3, and e = +1. Later, this approach 
was generalized to Bianchi I cosmologies in [370j . 

So far we have examined a classical scalar field on a given background. A quantum scalar 
field in a classical geometry was studied in [TH [75] where, inspired by the features of 
the mechanism of spontaneous symmetry breaking, the authors seek a solution in which the 
fundamental state of is given by 

where rj is the conformal time of an open Friedmann geometry given by 

ds"^ = a^{T]) [dr]"^ - dx^ - sinh^ x {dO^ + sin^ ^^0^)] ■ (69) 
For a massless field the equation of motion for the scale factor reduces to 

a" , , 

- = 1. 70 
a 

From the Lagrangian 

we obtain the equation of the scalar field 0, given by 

0" + 20'- + 2aaV' = 0. (71) 
a 

Compatibility of these two equations with the assumption in Eqn.f l68l) yields the relation 

(T2) 

For the scale factor as function of the Gaussian time t we obtain 



i{t) = Vt^ - L^, (73) 



*This potential interpolates between p = p for p << A, and p < for high densities. 
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where L is a constant and 

r-f + f = 0. (74) 

It was shown in [71] that the solution / = is unstable, while the solutions = 1 are 
stable. From the equation for g^,^ and specializing for fi = u = we obtain the value of the 
constant L in Eqn. fl73|l : 

L^ = — (75) 
24cT ^ ' 

which represents the minimum allowable value of the scale factor. From standard quantum 

field theory in curved spacetime, 

it follows that i?|o) = < 0, which shows explicitly the expected violation of the weak 

energy condition that causes the absence of a singularity in this model. Note that the 
gravitational constant in the vacuum state is renormalized: 

1 1 02 l2at^-n/2 



It follows that Kren < for < 2^ and K^en > for t"^ > thus showing that a change in 
the sign of the gravitational constant can be induced by the non-minimal coupling of scalar 
field with gravity, yielding repulsive gravity. 

The phenomenon of repulsive gravity can also be generated at a classical level by means 
of a non- minimally coupled complex scalar field [368] . The Lagrangian is given by 

C = d^<j)dy* - cr(0*0)2 - ^/?(0*0) + K'^R + Cm, 



where a is the constant that measures the auto-interaction of 0, and C is the matter La- 
grangian. The EOM following from this Lagrangian are 

□ + 2(70*02 ^ ^ 

6 

where 

R = k(^1- ^0*0) , (76) 

e,u = \ {d,<P*d,<P + d,<P*d,<P - g,,{dp<P*d^<P - a(0*0)2) + ^g^M<P*<P) - \{<P*^\,^ , 

and Tfj^u is the energy-momentum tensor associated to matter. From Eqn. fl76l) we see that 
the gravitational constant is renormalized at the classical level by the scalar field. In fact, 
as shown in [368 J. for the open FLRW metric @ the scalar field has three vacuum solutions: 



"^^This scenario does not work for the closed case. 
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= 0, and = ±7/a(t), where 7 is a constant. Only the nonzero solutions are stable, and 
they are also more favorable from the point of view of energy |368] . Since they are inversely 
proportional to a, it may be argued that the scalar field was in a nonzero vacuum in the 
early universe. Hence, 

-1 

where Oc = (/t/12(T)^/^ signals the change of sign of the gravitational interaction. Nonsingular 
solutions were obtained in |368] for matter given by radiation (p = t/a^): 

a{t) = —= cosht, 
v2 

where w"^ = ~ |Ke. This case reduces to the case without matter for e = 0. 



1-^ 



3.2 Dynamical origin of the geometry 

We shall see in this section that a cosmological scenario displaying a bounce arises in an 
extension of Riemannian geometry called Weyl Integrable Space-Time (WIST) [204]. 

Let us begin by recalling that one of the central hypothesis of General Relativity is that 
gravitational processes occur in a Riemannian space-time structure. This means that there 
exists a metric tensor g^^ and a symmetric connection V^^u related by 

In other words, the connection is metric and can be written in terms of the metric tensor 
as follows 

rv = {U^ la'-^igp,,. + gp.,, - 9,u,p]. (78) 

A direct method to deduce such metricity condition is given by the first order Palatini 
variation (in which the variation of the metric tensor and of the connection are independent). 
The starting point is the Hilbert action: 

S = j ^/^Rd^x. (79) 

In a local Euclidean coordinate system, 

^R^lu = ^^'^^J.a■,u — ^lu■,a■, (80) 

where the covariant derivative represented by a semicolon must be taken in the non- 
perturbed background geometry. From this equation it follows that 

5C = {R,, - ]^Rg^,)^5g^- + ^g^'SR,,. (81) 
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Correspondingly 



5S 



-9{R 



Hence, 



and we obtain 



-99 
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2' 



-99n;uSl--{ 



0, 



-99 



0. 



^2) 



(83) 



^4) 



After some algebra it can be shown that space-time has a Riemannian structure, that is, it 
obeys the metricity condition, 

9^.e■,a = 0. (85) 

The other equation that follows from the variational principle yields Einstein's equations. 
The lesson we learn from this calculation is that the structure of the manifold associated to 
space-time is not given a priori, but may depend on the dynamics. Surely, we should check 
whether the addition of matter alters this feature. The answer is not unique: it depends 
crucially on the way matter couples to gravity. There will be no modification to the precedent 
structure if we adopt the minimal coupling (that is, if the strong equivalence principle is 
valid). However, when the interaction is non-minimal, the geometrical structure obtained by 
the Palatini variation is not Riemannian in general. The simplest way to show this is with 
an example. Let us take the Lagrangian which describes the non-minimal interaction of a 
scalar field with gravity in the form: 



Lint = y^Rfif). 

Following the procedure sketched above we get: 



^6) 



int 



+ 



and it follows that 



-9 f (Rf^u - 
sf^9f9''] 



^R9,.)59''' 



'9f9n;Jl 



i7) 



{V^f{v)9'''}.^, = 0. 
This equation shows that the covariant derivative of the metric tensor is not zero but 

where Qf^ux = —(}^f),\9^iu- Taking the cyclic permutation of Eqn. (13.21) yields 



X 

M a 



Q\, - Q 



19) 



(90) 
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The equation 

9f^u;a = -(^T^f),\9fj,,y (91) 

shows that the structure generated by the Lagrangian fl86|) using the Palatini variation is 
not Riemannian but, as we shall see in the next section, a special case of Weyl geometry. 

3.2.1 WIST (Weyl Integrable Space Time) 

A Weyl geometry is defined by the relation |228j 

gtMu;a = WaQt^y. (92) 

This equation implies that there is a variation of the length of any vector under parallel 
transport, given by 

M = ioWf.Ax". (93) 

This property has the undesirable consequence that the measure of length depends on the 
previous history of the measurement apparatus, as pointed out by Einstein in the beginning 
of the past century in a criticism against Weyl's proposal for the geometrization of the 
electromagnetic field |106j . Einsteins remark led to the abandonment of this type of geometry. 
However, there is just one particular case in which this problem disappears: the so-called 
Weyl integrable spacetime (WIST). By definition, a WIST is a particular Weyl spacetime in 
which the vector Wfj_ is irrotational: 

It follows that in a closed trajectory 

Ai = 0, (94) 

which solves the critic raised by Einstein. From the definition given in Eqn. fl^^ it follows 
that the associated connection is given by 

c",u = {l^^}-l (^M^" + ^M^" - ^''a.u) . (95) 

Using this equation we can write the contracted curvature tensor R^^^ in terms of the tensor 
R^i, of the associated Riemann space constructed with the Christoffel symbols | | • We 
obtain ^ ^ ^ 

where the covariant derivatives are taken in the associated Riemannian geometry and □ is 
the d'Alembertian in the Riemannian geometry. Thus, for the curvature scalar, 

= /?-3n(^ + ^(^,A<^'^ (97) 

in which R is the curvature scalar of the associated Riemannian spacetime. 

The expressions in Eqns. (1961) and (1971) are very similar to those obtained by a conformal 
mapping of a Riemannian geometry as shown in Sec J3.4[ 
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3.2.2 WIST duality: the Weyl map 

A Weyl integral spacetime is determined by both a metric tensor and a scalar field. In 
[228] . Weyl introduced a generalization of the conformal mapping, which he called a gauge 
transformation, given by 

9t^u 9f^,y = e^gf,u, (f ^ (p = (p + (98) 

in which x is an arbitrary function. Under such transformations the affine connection and 
the curvature and Ricci tensors are invariant: 

f>{W) ^ p(iy) 

Note however that this is not the case for the scalar of curvature, which changes as 

This property has been used to construct gauge-invariant theories, as we shall see next. 

3.2.3 Invariant Action in WIST 

From the behavior of the geometric quantities under a Weyl map, it is not difficult to write 
an action that is invariant under the transformation given by Eqns.( l98|) : 



Sw = j v^e-^i?(^). (99) 
This Lagrangian can be rewritten in terms of the associated Riemannian quantities as follows: 

Sw = j ^e-^ (^R-?>nip+^- ip^^ip^^^ . (100) 
After some algebra, we arrive (up a total divergence) to the result 



Sw= I v^e-^(i?-^^,^^'^). (101) 



Note that the kinematical term of the scalar field for the scalar field appears with the "wrong" 
sign. This can be interpreted as a ghost field term hidden in the WIST structure. 

3.2.4 A particular case of WIST Duality 

Let us go one step further and add to the above Lagrangian a kinematical term: 



Sk = I V^e-^ip,,ip'^. (102) 
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If we restrict to the case in which x (given in Eqn. (p8|) ) is a functional of (p, it follows that 
the complete action 

S = JV^ e-'^(i?(^) +/3^,^^'n (103) 
is invariant under the restricted map 



(104) 



which is a special case of the general transformation fl^Sl) . In terms of Riemann variables, 



S 



-9 e 



(105) 



There are three invariants of dimension (length)^ that can be constructed with the inde- 
pendent quantities of a WIST geometry: R^^\ Lp" (pa, and ip^-a, where (pa = 'P,a- Now, since 
the covariant derivative ";" in the WIST spacetime can be written in terms of the Riemann 
covariant derivative (denoted by "||") as 

■,a = V \\a - <^,a, 

the three invariants reduce to two. The most general action can then be written as 



S 



Independent variation of the metric tensor and the WIST field (p yields 
(the operator □ is calculated in the Riemannian spacetime) and 



(106) 



(107) 



R^'^^u - \ i?^"^^ 9^.. + ^,^.■,u - 2(e - + (e - \)9,u^,c.r'' = 0. (108) 

This equation can be rewritten exclusively in terms of the associated Riemannian structure 



R 



1 



where 



-Rg^y-\ ip,^ y?,^ + — y?'"* g^u = 0, 



A^ = ^(4^-3). 



(109) 
(110) 
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3.2.5 A nonsingular cosmological model in WIST 

Let us now show how a nonsingular cosmological scenario in the WIST framework can be 
constructed, following |204] . We shall work with the standard form of the FLRW metric: 



dr"^ 



1 — er" 



(111) 



As in the case of a standard scalar field, the WIST configuration can be represented by a 
perfect fluid, so that Eg. (11 09 p becomes Einsteins equation for a perfect fluid with v'^ = 5q, 
energy density and pressure p^, given by 



2 -2 



(112) 



In this interpretation, the WIST structure is equivalent to a Riemannian geometry, satisfying 
the equations of General Relativity with a perfect fluid having negative energy density as a 
source. The gauge vector (px for this geometry becomes 



(113) 



where the dot denotes differentiation with respect to the time variable. Use of Eq. (11071) 
yields a first integral for the function ip{t): 



where 7 = constant. In turn, EE (I109p for the Friedman scale factor a{t) are 



+ e 



6 



0, 



2a d + + e — — {'paf' 







(114) 



(115) 



;il6) 



where e is the 3-curvature parameter of the FLRW geometry. From Eqns. (11140 and (I115P we 
see that (3 — 4^) = — 2A^ < 0, and an open Universe is obtained (i.e., e = —1). Combining 
Eqns. ([112]) and (11151) we get the fundamental dynamical equation 



Oo 

a 



1 4 



(117) 



with ao = [7^A^/6]^/"^. Before entering into the details of the solution of the system of 
structural and dynamical equations (I115P and (11140 . let us comment some of the consequences 
of this cosmological model and list some interesting results. 

Features of the model 

An immediate consequence of Eq. (lll7p is that the scale factor a{t) cannot attain values 
smaller than ao. Let us consider a time reversal operation and run backwards into the 
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past of the cosmic evolution. As the cosmic radius a{t) decreases, the temperature of the 
material medium grows. In Hot Big Bang models such increment is unlimited; in the present 
theory, on the other hand, there is an epoch of greatest condensation in the vicinity of the 
minimum radius oq. Close to this period, there occurs a continuous "phase transition" in the 
geometrical background: a Weyl structure is activated, according to Eg. (11 141) : the Universe 
attains the minimum radius at {t = 0), and consequently an unbounded growth of the 
temperature is inhibited. Notice that since the Universe had this infinite collapsing era to 
become homogeneous, in the present scenario the horizon problem of standard cosmology 
does not arise. 

For very large times, the scale factor behaves as a ~ t. Thus, asymptotically, the geo- 
metrical configuration assumes a Riemannian character (since — t) in the form of a flat 
Minkowski space (in Milne's coordinate system). Consequently, in the present model the 
evolution of the universe may be started by a primordial instability of Minkowski spacetime 
at the remote past, due to Weyl perturbations of the Riemann structure through Eq. flQlj) . 
In order to prescribe the behavior of these perturbations, knowledge of the time dependence 
of the gauge vector ipx is required. Since the WIST function ip has a maximum at t = 0, 
the largest deviation of the Riemannian configuration corresponds to the epoch of greatest 
contraction near to the value Oq. 

Stability of the solution 

Among the difficult questions concerning bouncing Universes, one may count the problem of 
their survival with respect to eventual metric perturbations (see Sect.f lTT]) ). We shall show 
that during the stage of greatest condensation the WIST model of the Universe is stable. 
Applying the perturbative scheme 



to Eqs. (lll4p and (IllSp . one obtains 



a + 6a, 



5ip = Sa, 
a^ 



~ 2^ 6a. 

a^d 



Hence, 



(50 = — ^ aSd, 

6d A 4l — i 

— ~ a [a - agj 2 
Far from oq (i.e., for large t) we have a » oq; then, 

6d 



6a 
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da ~ dt 



so with {6a) i being the initial spectrum of perturbations, one obtains 



Sa ~ {Sa)i exp[a 



The solutions of the system Eqs. (4.12) and (4.15) are therefore stable against metric per- 
turbations in the course of the infinite collapsing phase. 

The exact solution 

No closed solution can be obtained in terms of the cosmological time, so it is convenient to 
move to conformal time t], in which case the solution is easily shown to be 



where rjo is an integration constant. The following qualitative plot shows the difference 
between this bouncing solution and the radiation-dominated model in standard cosmology. 
The scale factor has a minimum for a = oq, which corresponds to rj = rjo. Thus the Universe 



Figure 1: The qualitative plot shows (in conformal time) the scale factor for the bouncing 
model given by Eqn flllSl and the scale factor for radiation in the SCM, a oc rj. 

had a collapsing era for rj < rjo, attained its minimum dimension at rj = tjq, and thereafter 
initiated an expanding era. Both the collapse and the expansion run adiabatically, i.e., at a 
very slow pace. 

The correlate behavior of the Hubble expansion parameter H = {a /a) helps to understand 
the model (Fig. [2]). Indeed, the Hubble parameter H is a smooth function of the conformal 
time 7] and does not diverge at the origin of the expanding era; quite on the contrary, it 
vanishes at rj = rjQ. The corresponding evolution of the Cosmos may be outlined as follows: 
the Universe stays for a long period in a phase of slow adiabatic contraction, until H attains 
its minimum value. Then an abrupt transition occurs: a fast compression turns into a fast 



a 



iv) = aoVcosh2(r7 - ?7o). 
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Figure 2: Plot of the Hubble parameter in conformal time for oq = 1. 



expansion up to the maximum of H, and afterwards the expansion proceeds in an adiabatic 
slow pace again. While this image supplies a picture of the behavior of an Universe driven 
by ip{t), it is however incomplete, due to the fact that the production of large amounts of 
matter and entropy has been neglected. This topic will be discussed in Sect J3.2.9|) . 

3.2.6 The WIST function (^(t): structural transitions 

According to the basic conception of the scenario presented above, the WIST function (p(t) 
governs the cosmic evolution. Taking into account the solution Eg. (11321) for the scale factor 
a{t), the first integral equation f lll4p yields for (p{t) the expression 

(119) 

The behavior of (p{t) is qualitatively portrayed in Fig. [31 along with 0. Note that when 
a —>■ ±oo (i.e., for large times), ip —>■ ±77r/4aQ =constant, which is consistent with the 
assumption that the Universe originated from a Minkowskian "nothing" state. The behavior 
of the time derivative (p = which appears in Eq. (11121) of the energy density p^, of 

the "stiff matter" state associated to the WIST field is also shown in Fig|3l Since this 
function has a strong peak in the neighborhood of the minimum radius Oq, the greatest 
deviation from the Riemannian configuration happens at this point. In this sense, a sort 
of "structural phase transition" takes place when the Universe approaches its maximally 
condensed state. The increase of the (negative) energy of the WIST "fluid" precludes the 
collapse to a singularity, reversing the cosmic evolution into an expansion. Note that the 
"kinky" aspect of the behavior of the WIST function ip{t) in Fig 13] suggests a similarity 
between the Weyl structural transition described above and the propagation of instantons 
in Euclideanized models of quantum creation (see Eqn. (lll7p ). 
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Figure 3: Plot of (full line) and (dotted line) in conformal time for oq = 1. 



3.2.7 WISTons and anti-WISTons: On the geometrization of instantons 



In the derivation of the solution of the WIST structural function (p(t) (given by Eq. flll9p ). 
no attention was paid to the sign of the constant 7. Since the only information we have 
about 7 is that 7^ = Goq/A^, according to Eqs. (11141) and (11171) . 7 can be either positive or 
negative: 



7 



Al 



Hence, Eqns. (lll4p and (I119P actually yield two equations, as follows: 



(±) 



(±) 

(fQ arccos 



Oq 

a 



0W 



7 



(120) 



(121) 



in which ip\^ = 'j^^^ /2al = ±1^3/2 |7|^"'^. Thus the amplitude of the solutions ip^^^ depends 
exclusively on the dimensionless parameter ^ (see Eqn. flllOp ). The plot of the WIST functions 
(t) and ip^^^ (t) is given by the mirror image of Figl3] with respect to the horizontal axis. 
Note, however, that the energy density p^p of the "stiff matter" state associated with the 
WIST field (f{t) is the same in both cases, since from Eqns. (I112p and (I114p we have 



(122) 



Thus, in spite of the fact that the pairs of WIST functions (v^*-^-*, 0''^'') have 
different characteristics, they induce the same type of cosmological evolution. Their only 
distinction, in fact, is connected to length variations, since according to Eq. (l93p one now has 



39 



It is interesting to observe that the system is invariant with respect to the time re- 
versal operation t — (— t) if ip^^^ is concurrently mapped into ip^~^ and reciprocally. In 
this sense, the WIST instanton-like functions ip^^^ and ip^~^ may be called "WISTon" and 
"anti-WISTon" solutions, respectively, since an anti- WISTon may be described as a WISTon 
running backwards in time. According to Eq J107[ WISTons are defined up to an additive 
constant. 

A closer inspection of the equations governing the behavior of (p{t) reveals an instanton- 
like behavior typical of nonlinear theories of self-interacting scalar fields. Of course, the 
root of such nonlinearity is the fact that (p{t) is taken as the actual source of the curvature 
of the metric structure, which in turn modifies the D'Alembertian operator □ due to the 
introduction of <y9-dependent terms. A direct way to clarify this issue is to make explicit, 
by means of a change of variables, the hidden nonlinearity of the system of equations of 
motion involving the scale factor a{t) and the WIST function ip{t). Define the new variable 
s{t) = Using Eqns. ffTOTll and ffTTij) . we have 

I ^t^^""f=n ' (123) 

Taking s{t) to represent a generalized coordinate associated with a one-particle dynamical 
system yields the conservation equation 

+ V{s) = , (124) 
in which the associate potential V{s) is given by 

4 4 ^ S" 

-0 S -7^53 

(125) 

with = 6A~^7^/^. Thus the evolution of field s is equivalent to a unit mass particle moving 
in a potential with vanishing total energy. Due to the nonlinear character of this potential, 
the instanton-like aspect of functions (y9'^^^(t) is not surprising. Figure H] shows the behavior 
of V{s). The potential vanishes at s = and at s^'' = 'y^^^aQ^ its extrema are at s = 0, and 
at s^r^^ = (2/3) ^/'^7(^)ao^ (which are minima). However, the system cannot remain at the 
stable states V{sm '') = (— |) 7^cto^5 since in this case s ^ 0; this in turn implies, of course, a 
nontrivial, evolving cosmic configuration. This nonlinear scheme provides a succinct picture 
of the evolution of the Universe: its development is initiated at s = (which corresponds 
to Minkowski space time at t — > —oo), attains its minimum radius a{t = 0) = at either 
s^'' or s^'' and returns back to s = (which now corresponds to a Minkowski spacetime at 
t +oo). According to whether the system proceeds along the right or the left branches 
(i.e., from s = to s^'* or s^'') of the figure, the cosmic evolution is driven by a WISTon or 
an anti-WISTon, respectively. 
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Figure 4: Qualitative plot of V{s). 



The appearance of instanton-like configurations is a direct consequence of the fundamen- 
tal dynamical equation (11171) . in combination with the "structural" equation (11141) which 
prescribes the degree of "Weylization" of space time. 

3.2.8 Weylization 

We shall see next that the "structural transitions" discussed above are equivalent to a quan- 
tum tunnelling process in models of quantum creation from "nothing" . Consider a generic 
Einstein equation for a Friedman scale factor, 



It was shown in |414] that a semiclassical description of a quantum tunnelling process is given 
by the bounce solutions of Euclideanized field equations, i.e, of field equations in which the 
time parameter t is changed into {—it). Applying such an Euclideanization procedure to 
Eq. fll26p . one obtains 



In the case of an e = +1 universe driven by a (positive) cosmological constant A = 3^^^ this 
approach was used in [415] to obtain, instead of the classical de Sitter solution, namely 



— e + - pa 



(126) 
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a = +e pa . 



(127) 
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cosh(^t). 



the solution 
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corresponding to a de Sitter instanton - a "kink" configuration- propagating with negative 
classical energy, which bounces at the classical turning point a = oq = (l/<j) interpreted as 
representing the tunnelling to classical de Sitter space from "nothing." 

Now consider Eqn. (11261) in the case of a closed Universe driven by the energy density 
p = Sla^/a^]. The euclideanized version of Eq. (ll27p gives 



But this is precisely the fundamental dynamical equation ( 1117P of the WIST cosmological 
scenario. In this way, an equivalence is established between the Euclideanization of a closed 
Universe model driven by a positive energy density and a "structural transition" to a Weyl 
configuration which results in an open Universe model driven by a "stiff matter" state of 
negative energy. Just as in models of quantum creation the propagation of an instanton 
is seen to represent the tunneling of the Universe from a primordial quantum "nothing" 
state, in the present scenario the propagation of a WISTon (i.e., a deviation of the Rieman- 
nian structure) is tantamount to the development of the Universe from a primordial empty 
Minkowski space. 

It has been argued that solutions obtained through Euclideanization are in fact non- 
realistic, since they are to be interpreted as instantons, field configurations which tunnel 
across a classically forbidden region. Other authors endorse the view that such solutions cor- 
respond to an actual primordial phase of the cosmic evolution in which the basic Lorentzian 
nature of spacetime is changed into an Euclidean one. According to the present model, 
a different interpretation may be ascribed to these solutions, since an enlargement of the 
spacetime structure to a Weyl configuration - in which the geometry is characterized by the 
pair {gf^u, V\) of fundamental variables - supplies, at least in a particular case, the same basic 
behavior. It then becomes possible to conciliate the opposing interpretations of an "abstract 
soliton configuration" [260j and of a truly observable Euclidean cosmic phase |261j . The 
WIST solution is observable, whereas its basic nature is always Lorentzian. It is the Rie- 
mannian character of spacetime structure that results altered; allegorically, the choice is no 
longer Euclid or Lorentz, but rather Riemann or Weyl. 

3.2.9 Solution with matter generation 

We have mentioned above that the the model must be improved by taking into account 
matter creation. A non-singular solution in WIST that incorporates the effect of the creation 
of matter on the geometry was studied in [327j . Friedman equation in conformal time is given 




by 



/2 2 ^ 

a — a = 

6 



(129) 



while the second EE is 




(130) 
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The conservation of the stress-energy tensor in the case of ultra-relativistic matter is 

(aVmy + 4(«%)' = 0. (131) 

A particular solution to these equations that describes creation of relativistic matter only 
around the bounce, and enters a radiation phase with a constant scalar field in a short time 
is given by the expression [327j 

a{r]) = f3^/cosh{2r]) + ko sinh(2r7) - 2A;o(tanhr7 + 1), (132) 

with /3 = ao/Vl ~ ^o, and < ko < 1/7. The dependence of ip with r] can be obtained 
from Eqns. (11291) and (I130p . An asymmetry is to be expected both in the scale factor and 
in (p, since the evolution of this universe starts from the vacuum and enters a radiation 
dominated epoch. This is pictured in FiglS Notice that since the scalar field tends rapidly 




Figure 5: Plot of a and ip for ko = 1/7 and = 0.93, values chosen by imposing that the 
solution in Eqn. (11321) enters the radiation era for t ^ 10^^ seg. 

to a constant value, the production of matter (controlled by ip', see Eqn. (11311) ) stops soon, 
and the model enters a radiation phase without the need of a potential. In this sense, this 
solution describes a hot bounce, as opposed to cold bouncing solutions, which do not enter 
the radiation era unless they are heated up |328j . Another nice feature of this solution is 
that the scalar field (formally equivalent to the dilaton of string theory) goes automatically 
to a constant value for ^ oo, in such a way that the solution could be taken as the leading 
order of a perturbative development (as is the case in string theory). Again, no potential 
was needed in order to display this feature. 

3.3 Scalar-tensor theories 

Scalar-tensor theories are a generalization of the Brans-Dicke Lagrangian [78] , in which the 
constant appearing in the kinetic term of the scalar field becomes a function of 0. Among 
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the possible Lagrangians to describe these theories, one possibihty is [79] 



L = -f\<P)R + ^0,^0-'^ + 167r/:^,tter, (133) 



where the scalar field cj) couples non- minimally with the curvature through /(0). With the 
redefinition ip = f{4>), the Lagrangian becomes 

C = -ipR + '^<f,f,ip'^ + 167r£^atter, (134) 

with uj{(p) = \ f I f"^ and = df /dip. Brans-Dicke theory is a special case of this Lagrangian, 
/(0) oc (jP' which entails uo =const. This Lagrangian also describes the gravity-dilaton sector 
of low-energy string theory for = — 1 [33]. The differences between the two Lagrangians 
have been analyzed in [80] . Following the results of the discussion presented there, we shall 
use Eqn. 01341) as the definition of scalar-tensor theories. 
The equations of motion corresponding to Eqn. (11341) are 



Rnv = - ^df^uT) - '^^^^(p,f,p,u - ^^,f^;iy - ^9tiuO¥>, (135) 



[3 + 2u;{p:>)pp = T -LJ^p^^p''^. (136) 
Eqn. (ll35p suggests that it may be possible to find solutions in which matter satisfies SEC, 



but the whole r.h.s. is such that R^yv'^v^ > This implies, via the singularity theorem 
given in Sect.f ll.lT ) that nonsingular solutions may exist in scalar-tensor theories. Using 
Eqn. (1135]) . the inequality Rf^uV^v" > translates for the fiat FLRW case and EOS p = \p 
to 

1/1 oW^"'"^ f ^ ^' \ V 



4 



--(l + 3A)p-— --^ J -->0- (137) 

Solutions satisfying this constraint, and hence exhibiting a bounce, have been presented in 
J95j. for e = in the cases of vacuum and radiation (for which T = 0, see r.h.s. of Eqn. ( 11361) ) 
|. With these restrictions, Eq. (11361) written in conformal time takes the form 

a 3 + 207 

which integrates to 

^'a^ = ^==, (139) 

where A is a constant. Introducing the variable y = pa? and using Eq. (11391) . the Friedmann 
equation takes the form 

y'^ = ATy + A\ (140) 



"^^The same happens in some wormhole configurations in Brans-Dicke theory. See j87j . 
■^^A shadow of doubt has been cast on these results in [52], where it was shown that gravitons would still 
see a singularity, even if the rest of matter does not. 
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(r > is a constant coming from energy conservation) yielding for y{T]), 

y{i]) = A{r] + r]o) (141) 

in the case of vacuum, and 

y{v) = T{v + vo)' - ^ (142) 
in the case of radiation. Dividing now Eg. (11391) hj y = ipa? we obtain 



fV^WTl ^^^^r^ (143) 

J ^ J y{v) 

If this equation is invertible, we could obtain from it (p = ip{r]), and then y = Lpa^ yields 
a{j]). To integrate Eq. (ll43p . we need to specify the function ujiLp). The choice in [95] was 

2u{ip) + 3 = 2/5 (^1 - , (144) 

where a, /? > and are constants. With this choice of u;, Eq. (11431) can be solved for ^[j]) 
in the cases a = (which corresponds to Brans-Dicke theory), a = 1 and (5 = —\ (which 
defines a theory introduced by Barker [SS|), and a = 2. The latter was studied in [55]. The 
solutions for the vacuum case are given by 

, ^ A{v + vo){l + {v + Vo)') 



<^c(^ + ^o) ' 

i + {'n + voY 

with A = a/3/2/3. These solutions were shown to be nonsingular for (3 < 3/2. Hence the 
radiation solutions (which approach those for the vacuum for r/ — > [95]) are also nonsingular. 
All the solutions for a = 2 approach the FLRW radiation regime at late times because cp 
tends to a constant, and then uj{ip) oo, but in order to be in agreement with solar system 
experiments, a must be greater than 1/2 [95j. 

The case of stiff matter (defined by p = p) sourcing the scalar field was studied in |84j . 
Since the density of a barotropic fluid {p = (7 — l)p) evolves as p cx a"^''', this kind of matter 
is expected to dominate at early times, and the associated solutions give information about 
the early evolution of the universe. One of the results in [M] is that a necessary condition 
for d = when spatial curvature is negligible is = —6Mip/A, where A and M are positive 
constants, yielding a negative kinetic term for ip (see Eqn. (ll34|) . A thorough qualitative 
study of the case in which u}{ip) is a monotonic but otherwise arbitrary function of (p was 
presented in [53] , where the existence of nonsingular solutions in theories which agree with 
GR in the weak field limit was proved. 

The first term on the left hand side of Eqn. (11350 suggests that the gravitational constant 
is not actually a constant but varies with tp~^. Based on this idea, a generalization of 
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scalar-tensor theories (the so-called hyper-extended scalar-tensor) was advanced in [85]. The 
Lagrangian associated to these theories is given by 



C 



+ IGvriZmatter; 



(147) 



which reduces to Eqn. (11341) when G{ip) = 1/ip. Sufficient conditions on G{ip), uj{ip), and their 
derivatives in order to have bouncing cosmological solutions were given in [HS], generalizing 
the work of [81] for the case of ST theories. 

Another descendant of the original ST theory are the multiscalar-tensor (MST) theories 
|416] . which are the generic product of a compactification process of a higher-dimensional 
theory. The scalar content of a given MST theory depends on the details of the internal 
manifold that results from the compactification (usually gauge fields are set to zero in cos- 
mological applications). Typically, one or more fields are associated to the size of the extra 
dimensions. In string theory, the coupling constants depend on the expectation value of 
massive scalar fields (called moduli fields) also associated with the size and shape of the 
extra dimensions, the most popular example of them being the dilaton. The moduli are an 
inescapable ingredient of string theory, hence several problematic issues raised by them must 
be confronted, such as stabilization, overcritical density, and violations of the Equivalence 
Principle. Cosmological solutions of low-energy string theories have been extensively studied 
(see |90] for a review). Needless to say, the results depend on the field content, which in turn 
depends on the given string theory under scrutiny. 

A possible way to parameterize an action of a MST theory is |1U7] 



This Lagrangian represents pure multidimensional theories when =constant, x = constant, 
and uj = {1 — d)/d, where d is the number of compactified dimensions (assuming that they 
have the topology of a torus). The same case but with ip ^ constant and n = —2/d + 1 
corresponds to a two-form gauge field in higher dimensions. If this field is conformal, it is 
associated to a (rf + 4)/2-form, leading to n = —2/d. In the case of string theory, uj = — 1, 
and the field ip is associated to a three-form field Hfj_^x, leading to n = —1. The scalar x 
is related to another three-form field coming from the R-R sector of type IIB superstring 
theory. 

The existence of bouncing solutions for this Lagrangian in vacuum and in the presence 
of radiation for the FLRW geometry for all values of the three-curvature and for arbitrary 
values of u and n has been studied in |107j . The results show that generically there is bounce 
for n < 1 and u < 0. 

3.3.1 Corrections coming from String Theory 

Superstring theory is a candidate for a unified theory of the fundamental interactions, includ- 
ing gravity [310\ . Since the fundamental objects in this theory are at least one-dimensional. 




matter • 



(148) 



46 



geodesies of point particles are replaced by world- volumes. It is a valid question then to ask 
whether string theory has anything to say about the singularity problem. In this regard, 
it must be noted that in string theory, the gravitational excitations are defined on a fixed 
metric background. Since singularities in general relativity are boundaries of space-time, 
which are a consequence of the dynamics governing its structure, a fixed manifold is cer- 
tainly a restriction. Yet another difficulty is the breakdown of string perturbation theory 
in the regime of interest [301] . However, we have seen in the previous section that the in- 
corporation of the massless degrees of freedom (correspongind to the lowest order EOM), 
which applies on scales below the string scale and above those where the string symmetries 
are broken, may smooth out the singularity. One could go further and include higher-order 
corrections in the action of string theory. There are two types of corrections. First, there are 
the classical corrections arising from the finite size of the strings, when the fields vary over 
the string length scale, given by = ^/ol' . These terms are important in the regime of large 
curvature, and lead to a series in a' (the inverse of the tension of the string). Then there 
are the loop (quantum) corrections. The loop expansion is parameterized by powers of the 
string coupling parameter e'^ = S'string; which is a time-dependent quantity in cosmological 
models. In the so-called strong coupling regime, the dilaton becomes large and quantum 
corrections are important. 

The effective action at the one- loop level is given by (see for instance [96j ) 

S = J d'x^g || + \{S/<pf + \{S/a f + l[Ae^ - 5e(a)]i?^^} , (149) 

where is the dilaton, a is a modulus field, and X = {g is the string coupling), 6 is 
proportional to the 4-d trace anomaly, and ^(a) = ln(2e'^?7^(ie°")), where rj is the Dedekind 
function. The correction to the gravitational term is given in terms of the Gauss-Bonnet 
invariant, 

^GB — -K/x!/kA-K — ^rc^^n + n . 

The EOM that follow from this action in the case of a FLRW fiat spacetime with the metric 
g^, = diag(l, -e2-5,,) are ^96] @ 

3cj2 _ - 1^2 ^ 24/cI;^ = 0, (150) 

2u + 3ilj^ + + + Wfu^ + Sfdj'^ + 16fuu = 0, (151) 

a + 3ua + 6^dj'^(dj^ + u) = 0, (152) 
oa 

4) + 3u^-3\e'^u^{u^ + u) =0, (153) 

where / = ]^(Ae'^ — 6C,{a)). These equations are not linearly independent due to the conser- 
vation of Tfj^^. 

*^See [53] for the case of nonzero spatial curvature. 
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It was shown in [96] that there are solutions with bounce for S < 0, which interpolate 
between an asymptotically flat and a slowly expanding universe with a period of rapid 
expansion. The bounce is essentially due to the violation of the strong energy condition 
by the modulus field (the dilaton playing an unimportant role). In a subsequent paper p7] 
it was shown that non-singular solutions can be obtained under the assumptions that ^ is 
a smooth function that has a minimum at some point ctq, and grows faster than o"^ for 
a ±cxD, and 6 > 0. However, these solutions were later shown to be generically unstable 
for tensor perturbations [102J. Less symmetric models (Bianchi I |1UU I| and Bianchi IX [IIJT]) 
were also studied for this action, confirming the findings of |102j . 

Another attempt to avoid the singularity is to consider the effect of matter terms to the 
action of string theory. In |104j an action including dilaton, axion and one modulus field 
was considered along with matter (radiation or a "stringy" gas) and higher-order dilaton 
corrections in a fiat FLRW background in d dimensions. In this case, the results of |1U4] 
show that the energy densities of matter, axion and modulus are strongly suppressed in the 
inflationary phase driven by the dilaton, and hence the higher-order corrections coming from 
this field take the system through the graceful exit. 

Yet another model inspired in string theory is the so-called ekpyrotic universe and its 
extension, the cyclic universe which will be discussed in Sect 110231 

3.3.2 String Pre-Big Bang 

A very-well developed example of the string cosmology approach is the so-called "pre-big 
bang" |91j, which we shall call "string pre-big bang" (SPBB), to differentiate it from similar 
models not coming from string theory (see Sect J3. 21) . There are two properties of string 
theory that can be expected to play an important role in cosmology [S2]- First, in the short- 
distance regime, a fundamental length is expected to arise, thus introducing an ultraviolet 
cut-off and bounding physical quantities such as and a. Hence a bounce may be expected. 
Second, as we discussed before, at lower energies, the action of string theory is not Einstein's 
but a (multi)scalar-tensor theory, where one of the scalar fields is the dilaton, which controls 
the coupling constants. If these are really constant today (see ^3]), the dilaton must be 
seated at the bottom of its potential, but it may have evolved in cosmological times. The 
idea of the SPBB is that during the cosmological evolution, the kinetic term of the dilaton 
drove a period of deflation (or inflation, depending on whether we consider the Einstein 
frame or the string frame) "before the big bang" (that is, in the contracting phase j^, which 
can solve the horizon and flatness problems |397j . In this approach, the universe starts from 
a perturbative state, passes through a high- curvature and high-coupling stage, and then 
(hopefully) enters the radiation-dominated FLRW evolution. Duality symmetries present in 
the low-energy action of string theory are invoked to support this line of reasoning p3]: in 
the isotropic case, the gravidilaton EOM in the FLRW setting are invariant under a time 
inversion, 

-t^H ^ -H, 

^^This idea was also suggested in [204] . 
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and under the duality transformation 



a 



a = a 



-1 



(j)-^(j) = (j)~Q\na. 



(compare with the Weyl transformation, Eqn. fl98|) ). These transformations relate four branches 
of the solution (PBB, and post-big-bang expansion and contraction). In particular, to any 
expanding solution with decreasing curvature (such as those in the standard cosmological 
model), duality associates an accelerated contracting solution (see Fig |3.3.2|) . It is this pair- 
ing (which is possible only in the presence of the dilaton) that supports the whole idea of the 
SPBB. One of the issues of this idea is the joining of the two phases through the putative 



Figure 6: The four branches of the low-energy string cosmology backgrounds. Taken from 



singularity (the graceful exit problem). It has been proved in [313j that the graceful exit 
transition from the initial phase of inflation to the subsequent standard radiation dominated 
evolution must take place during a "string phase" of high curvature or strong coupling is 
actually required. The corrections to the lowest-order lagrangian can be parameterized as 



H 






where 




(154) 
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and Lq designates the quantum loop corrections. Several forms of Lq were studied in |311j . 
The existence of a bounce in the Einstein frame, yielding a solution to the graceful exit 
problem, was shown by numerical integration of the EOM in [3TT] for the case Lq = — 2(V0)^, 
Cq = — 2(V0)^ + -R^/3, and for the two-loop correction Cg = 2e'^i?|.^, in all cases by choosing 
the appropriate sign for the correction. 

An even more general form of the corrections was studied in [312j . where Cc was given 

by 

and 4:13 + 20 + d = —4a {G^'^ is the Einstein tensor). The quantum corrections were included 
by adding a suitable power of the string coupling, so the total effective Lagrangian is given 
by 

C = R+ {d^cpf + Cc + Ae'^Cc + Be^'^'Cc, 

and the parameters A and B set the scale for the loop corrections. Solutions with graceful 
exit were found in |312j for a large range of parameters, but it is very hard to obtain the 
transition in the weak coupling regime, whilst keeping the loop corrections small. 

A problem that remains to be solved is the stabilization of the dilaton to a constant 
value (otherwise there would be violations to the Equivalence Principle and to the observed 
"constancy of the coupling constants"). This was achieved in the previously mentioned 
articles in a number of ways: 1) by introducing by hand a friction term in the equation 
of motion of the dilaton, and then coupling it to radiation in such a way as to preserve 
overall conservation, 2) by "turning off" by hand the quantum Lagrangian by means of a 
step function, and 3) by the manipulation of the sign and size of the higher-loop corrections. 

3.4 Appendix: Conformal Transformation 

Consider the map 

g^^i, (x) = (x) g^,y (x). (155) 
Then, for the contravariant components: 

-g^""" {x) = n'\x)g^'' {x). (156) 

The conformal transformation of the connection is provided by 

f = r". + ^ {^,, K + - 9'" 9,u) , (157) 

and for the curvature tensor: 

where 
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Contracting Eqn. fllSSp we get 



i?", = n-2i?%--Q",--Q5^, (159) 

and contracting again, 

R = n-^[R + Qn-^nn]. (leo) 

A direct comparison of this conformal scalar of curvature and the Weyl scalar equation (l97j) 
shows that they coincide (up to a multiphcative factor) if we set 



and Eqn. fll60p takes the form 



which is exactly the transformed of the Ricci scalar for the WIST: 

R = e^R^^\ 



4 Maxwellian and Non-Maxwellian Vector Fields 
4.1 Introduction 

The model described by the FLRW geometry with Maxwell's electrodynamics as its source 
displays a cosmological singularity at a finite time in the past |210] . However, this is not an 
intrinsic property of the combined electromagnetic and gravitational fields. Indeed, mod- 
ifications of Maxwell electrodynamics (or, generically, massless vector field dynamics) can 
generate non-singular spatially homogeneous and isotropic (SHI) solutions of classical GR. 
We shall examine here two modifications that are relevant to the singularity problem: 

• The non-minimal coupling of the EM field with gravity, and 

• the self-interaction of the EM field. 

These modifications will be introduced by means of Lagrangians which depend nonlinearly 
on the field invariants or on the space-time curvature. In both cases, the singularity theorems 
(see ChJT]) are circumvented by the appearance of a large, but nevertheless finite, negative 
pressure in an early phase of the SHI geometry. 
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4.2 Einstein-Maxwell Singular Universe 



The fact that Maxwell electrodynamics minimally coupled to gravity leads to singular models 
for the universe in the FLRW framework is a direct consequence of the singularity theorems 
(see Ch{T]) . Essentially, this can be understood from the examination of the energy conser- 
vation law and Raychaudhuri equation, as follows. To be consistent with the symmetries of 
the SHI metric, an averaging procedure must be performed if electromagnetic fields are to 
be taken as a source for the EE [B]. As a consequence, the components of the electric £"1 and 
magnetic Mi fields must satisfy the following relations: 



^. = 0, ^ = 0, ^,^ = 0, (161) 

^ = -\'^^9^3. (162) 

M,M, = -^M'^g,,. (163) 

The symmetric energy-momentum tensor associated with Maxwell Lagrangian is given by 

E^, = F^^F'', + -^Fg^,, (164) 

in which F = F^i, F^'^ = 2{Jf'^ — S""^). Using the above average values it follows that the T^y 
reduces to a perfect fluid configuration with energy density and pressure given by 



E^,u = (p^ + p^) tv '"u - Qf^u, (165) 

where 

P^ = Sp., = -{<ff^ + M^). (166) 

The fact that both the energy density and the pressure in this case are positive definite for 
all values of t implies the singular nature of FLRW universes. In fact, the solution of EE for 
the above energy-momentum configuration gives for the scale factor the singular form |225j 



a{t) = ^alt - et2, (167) 

where Qq is an arbitrary constant. We conclude that the space-time singularity in the 
Einstein-Maxwell system is unavoidable. 



4.3 Non-minimal interaction 

Most of the articles concerning the interaction of Electrodynamics with Gravitation assume 
the principle of minimal coupling, which is a direct application of the strong form of the 
Equivalence Principle. In the absence of stringent limits from observation, ideally we should 
keep an open mind and consider other possibilities. Non-minimal coupling of the EM field 
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with gravity has recently been apphed in cosmology, following the trend initiated by scalar 
field theories interacting conformally with gravitation. These opened the way to the exam 
of more general theories, such as those in which curvature is directly coupled with the fields. 

There are seven possible Lagrangians for the interaction of the EM field with Gravity 
which can be constructed as linear functionals of the curvature tensor. They are divided in 
two classes. Class I is given by: 

A = RA^A^, 
C2 = R^uA^Ar 

These two Lagrangians are gauge dependent but no dimensional constant must be added 
since they already have the right dimensionality. As shown in |177j the EOM obtained from 
C2 in Einstein's gravity with the addition of a kinetic term for A'^ do not admit a FLRW 
solution. Thus, in the following we shall limit our analysis to £1. 
In Class II, there are five Lagrangians : 

C3 = R F^yF^^ , 
= R F^ F"" 

Cr = Wap,u (168) 
where W""^^^ is the Weyl tensor and the star in the Weyl tensor means 

Wo.,,. = W^^,^ = W,^^^ = \ri:;W,^,.. 

These Lagrangians are gauge independent but they all need the introduction of a length £0 
in order to have the correct dimensionality. 

Another Lagrangian sometimes studied in the literature that is not explicitly contained 
in this list is 

£3 = ^F^^F^"'. 

However, is not independent of {Ci,...Cj). Indeed, the double dual -R*^^* satisfies the 



identity 



or, equivalent ly. 



^aPiit ~ ^a|3^J.u - 'iWa/Sfiu ~ 2 ^ 9al3fiu, (169) 



- ^ Rga^3^^u■ (170) 
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Thus, 



Hence, £« = -Ca - | £3 - 2^5 . 



Raugpfi — Rf3fj.gau) F'^'^F^^ . 

I £3 — 2£5. 



4.4 An example of a non singular universe 

The first example of a nonsingular universe driven by the nonminimal couphng of EM and 
gravity was presented in |177j . using the £1 of the previous section: 

C= R-^F^''F^, + f3RA^A^^. (171) 

As mentioned in Sect J4.2l in order to obtain a SHI geometry in the realm of General Relativity 
having a vector field as a source, an average procedure is needed. In the present non-minimal 
case there is another possibility, which we shall now explore. Since this theory is not gauge- 
invariant, it is possible to find a non-trivial solution for such that F^'^ vanishes. 
The equations of motion that follow from the Lagrangian (11711) are: 

{1 + 13 A^) {R^, - ^ Rg^,) - (3 UA^ g^, + {3 {A^),^,, + (3 RA^A, = -E^,-, (172) 

Fl'J = -2I3RA^'. (173) 

From the trace of (I172p it follows 

R=-'il3uA^, 

which when inserted in the equation of evolution of the electromagnetic field yields a non- 
linear equation: 

F^''.^-Q[3'^{nA^)A'' = Q. (174) 

The non-linearity induced by the non-minimal coupling with gravity is a generic feature for 
any field. To obtain a solution in which the geometry is nonsingular for a SHI geometry 
without imposing an average on the fields |177j we can consider the case in which F^y is 
zero. This is possible due to the explicit dependence of the dynamical equations on the 
vector A^. We take the vector field A^ of the form 

A, = A{t) 5°. (175) 

Defining the quantity Q by 

n{t) = l + f3A'^, (176) 
the set of equations (11 72111 73 p in a FLRW geometry reduces to the following: 

3^ = -§. (177) 
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□ fi = 0. (179) 

The last equation implies that a^dQ/dt is a constant. Thus we set dfl/dt = ba~^. A 
particular solution of this set of equations for e = — 1 is given by |177] 

AHt) = 1-^ (180) 



a{t) = sjf' + al (181) 

where ao is a constant that measures the minimum possible value of the scale factor. When 
tto = the system reduces to empty Minkowski space-time in Milne coordinates. For ao 7^ 
this model represents an eternal universe without singularity and with a bounce The 
system (1177111791) can be written as a planar autonomous system, which was examined in 
[334j . . . . Notice that in recent years theories with negative energies have been examined in 
a cosmological context |325] . One way to achieve this goal is by introducing an ad-hoc term 
in the Lagrangian with the wrong sign. In the case of a scalar field this is given as 

1 
2 




S= I ^{R--d^^d^^). (182) 



A fluid with this odd feature can also be obtained by the non-minimal interaction of a vector 
field with gravity. Indeed, the solution presented in the precedent section can be interpreted 
as a perfect fluid with negative energy. The equations of motion presented in |177j can be 
re-written in the form: ^ 

were f2, given by Eqn. (11761) . depends only in time . The structure of the corresponding 
system of equations is equivalent to the equations of General Relativity in the SHI geometry 
having as its source the energy-momentum tensor of a perfect fluid with negative energy 
density and pressure given by 

J>=^P = -4 (184) 

In this way, fluids with the "wrong" sign in Einstein's equation can be interpreted as vector 
fields with non-minimal interaction with gravity. 



4.5 Nonlinear electrodynamics 

As pointed out in the introduction of this Chapter, linear electromagnetism unavoidably 
leads to a singularity. This situation changes drastically in the case of non-minimal cou- 
pling. In this section, we shall deal with another type of theories, in which it is the non- 
linearity of the self-interaction of the EM field that provides the necessary conditions for 

^°This form of the scale factor is similar to Melnikov-Orlov geometry [71], the difference being in the 
interpretation of the minimum radius oo and the source of the curvature. 
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a cosmological bounce to occur. The theories that will be examined are described by La- 
grangian which are arbitrary functions of the invariants F and G that is £ = C{F, G), where 
F = F^jjF^^^^G = \Tiap^ivF°'^ F^^^ . Their corresponding energy momentum tensor, computed 
from Eqn.Q yields 

T^, = -4 Cf F^'^F^, + {GCg - C) g^,, (185) 
where Ca = dC/dA, with A = F,G. It follows that 

p = -C + GCg - ^Lp<g'^, (186) 

V = H-GLG--{^.y^'^ - ^'^)i:.F- (187) 

3 

We shall start our analysis by studying a toy model generalization of Maxwell's elec- 
trodynamics generated by a Lagrangian quadratic in the field invariants as in [201] . that 
is: 

L = + + (188) 
where a and /5 are dimensionfull constants [^. 



4.5.1 Magnetic universe 

In the early universe, matter behaves to a good approximation as a primordial plasma 
|212[ 1213] . Hence, it is natural to limit our considerations to the case in which only the 
average of the squared magnetic field M''^ survives |211l 1212] . This is formally equivalent 
to put ^ = in (11621) . and physically means to neglect bulk viscosity terms in the electric 
conductivity of the primordial plasma. 

The Lagrangian (I188P requires some spatial averages over large scales, such as the one 
given by equations fll6ip - (ll63l) . If one intends to make similar calculations on smaller scales 
then either more involved Lagrangians should be used, or some additional magnetohydro- 
dynamical effect [229j should be devised in order to achieve correlation [230j at the desired 
scale. Since the average procedure is independent of the equations of the electromagnetic 
field we can use the above formulae (I16ip - fll63l) to arrive at a counterpart of expression fll65p 
for the non-Maxwellian case. The average energy-momentum tensor is identical to that of a 
perfect fluid (11650 with modified expressions for the energy density p and pressure p, given 
by 

p = ^Jf^l-8aJf^), (189) 

p = -J^^l-40ajr^). (190) 
6 

Inserting expressions (1189p -( fT90l) in the conservation equation ([2]) yields 

^ = ^, (191) 

^"'^If we consider that the origin of these corrections come from quantum fluctuations then the value of the 
constants a and f3 are fixed by the calculations made by Heisenberg and Euler. 
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a'^-^ 1-°-^ (192) 



where J^q is a constant. With this result, equation (jl]) leads to 

Since the right-hand side of equation (11921) must not be negative it follows that, for a > 
the scale factor a{t) cannot be arbitrarily small regardless of the value of e. The solution of 
Eqn. (ll92p is implicitly given as 

dz 

t = ± , . (193) 



where a(0) = ao- The linear case described by Eqn. (ll67p can be regained from Eq. (11931) by 
setting a = 0. For the Euclidean section, expression (11931) can be solved asE§ 



a2 = ^oy^(t2 + 12a). (194) 

From Eqn. (ll9ip . the average strength of the magnetic field ^ evolves in time as 

^2=3 1 _ ^ 

2^2 + 12 a ^ ' 

Expression (I194p is singular for a < 0, as there exist a time t = ^ —Via for which a(t) is 
arbitrarily small. Otherwise, for a > at t = the radius of the universe attains a minimum 
value (see FigjT]) ao, given by 

ag = J^V8a, (196) 

which depends on M^. The energy density given by Eqn. (11890 reaches its maximum value 

pmax = 1/Q4a at the instant t = tc, where 

= yi2a. (197) 

For smaller values of t the energy density decreases, vanishing at t = 0, while the pressure 
becomes negative (see FiglH], left panel). Notice that we have a minimum of a{t) along with 
a minimum of the energy density, entailing a violation of the NEC condition, in accordance 
with the first row of Tabl dl.2l 

Only for times t ^ v^4a the non-linear effects are relevant for the normalized scale-factor, 
as shown in Figure [H left panel. Indeed, the solution (11940 yields the standard expression 
(11671) of the Maxwell case at the limit of large times. Notice that the energy-momentum 
tensor (11851) is not trace- free for a 7^ 0. Thus, the equation of state = p-yip-y) is no longer 
that of Maxwell's; it has instead a term proportional to the constant a, that is 

1 1 6 

P = -p-^aJi-\ (198) 



^^Nonsingular solutions in Bianchi universes with nonlinear electrodynamics as a source were studied in 
[223] . 
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Figure 7: Plot of the scale factor as a function of t for different values of e and a. Taken 
from [20T] . 



This scenario has been generalized in several ways in [224J . First, the general expression for 
the scale factor was shown to be 

a{t) = ao(4a^t^ + AaoPot + (199) 

where 



Eqn. fll94p follows as a particular case from Eqn. fll99p . which describes a bounce with 

Omin = «o(8co'^^)^/'', tmin = " /?o/ (2ao) , ^nin = r— , Pmin = 0. 

2v 2« 

Solutions of this model with the addition of a cosmological constant A were also discussed 
in |224] . It was shown that nonsingular solutions are possible both for a constant A, and for 
certain choices of A = A(t). 

4.5.2 Born-Infeld electrodynamics 

A widely studied EM theory is that proposed by Born and Infeld, with Lagrangian 

Lbi = f 1 - Vx) (200) 
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Figure 8: Left panel: time dependence of p and p for e = and a > 0. Right panel: scale 
factor for e = and a > (full line), and for Mawxell's case (dashed line). Taken from [201] . 



where 

Note that, following Born-Infeld's original work, a constant term has been added in the 
Lagrangian in order to eliminate a cosmological constant and to set the value of the Coulomb- 
like field to be zero at the infinity. Using equation (11861) for the energy density we obtain 

P=4J^ (202) 



and for the pressure 

A straightforward calculation of p + Bp shows that this theory cannot yield a nonsingular 



P=^\^X-(3' + -----^] (203) 



universe. 



4.5.3 Bouncing in the Magnetic Universe 

The "magnetic universe" displays a very interesting property due to the nonlinear dynamics: 
its energy density can be can be interpreted as composed of k non-interacting fluids, in the 
case in which the dynamics is provided by the polynomial 

C = Y,CkF\ (204) 

k 

where A; G Z. The conservation of the energy-momentum tensor projected in the direction 
of the co-moving velocity f ^ = 5 q yields 

p + {p + p)e = Q. (205) 
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From the expression for the energy density and pressure given in Eqns. fll86p and (11871) with 
<f = we get that p = J2kPk^ and p = YjkPk where 



Pk = — Cfc2''J^^^ 
Pk = Ck2'j^^''(l 




(206) 



in such a way that we can associate to each power of k an independent fluid characterized 
by Pk and pk, with an EOS 



Inserting the total energy density and pressure (from the sum of pk and pk in Eqns.f l206l) 
and (12061) ) in the conservation equation (12051) we obtain 



The important result that this equation shows is that each fc-fluid is separately conserved, 
since the dependence of the conservation equation on the specific form of the Lagrangian 
factors out, in such a way that evolves with the scale factor as 



for any C of the form given in Eqn. (12041) . 
4.5.4 Two-fluid description 

It follows from equations (I189p . (11901) and (11911) that in the case of the nonlinear Lagrangian 
given by Eqn. (11881 it is not possible to write an equation of state relating the pressure to 
the energy density. This is a drawback if we want to use a fluid description of the averaged 
electromagnetic fleld. In order to circumvent such difficulty a two-fluid description can be 
adopted, because of the remarkable fact that there exists a separate law of conservation for 
each component of the fluid, as we saw above. The fact that the dynamical equation for 

factors (see Eqn. (12071) means that the fluids are conserved independently: the energy- 
momentum tensor can be separated into two pieces, each representing a perfect fluid which 
is conserved independently. In other words, there is no interaction between fluids 1 and 2. 
We shall see in Section 111.21 that the analysis of the stability of the non-singular universe 
described in this section is more transparent when using the two-fluid description. This 
case can be generalized to a multi-component fluid, but we shall restrict here to the 2-fluid 
application for a pure magnetic fleld. 

In order to get a better understanding of the properties of the cosmic geometry controlled 
by the magnetic fleld let us analyze the case in which the spatial section is closed (e = 1). 





(207) 




(208) 
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The crucial equations for such analysis are the conservation law, the Raychaudhuri equation 
for the expansion and the Friedman equation, that is: 

p + (p + p)^ = 0, (209) 

e + \e^ = -\{p + ?>p). (210) 



3 2 

(211) 



In the magnetic universe we have 

A necessary condition for the existence of a bounce is given by the vanishing of the expansion 
factor for a given value of t. This leads to an algebraic equation of third order in x = a^: 

^3_:^^2^4^^^4^^^ (213) 
6 3 

Using the fact that a is a very small parameter, it can be shown that this equation has three 
real solutions. Two of them are positives and the third is negative. Thus we retain only the 
positive solutions which will be called Xi and X2. The important quantity for our analysis 
is contained in the expression 

^2 

Pb + 3pfe = -f (x2-24«jro'). (214) 

Thus, at one of the points, say Xi there is a local maximum for the scale factor; and at the 
other, X2 there is a minimum for < 24a^^. Note that at the bounce (where = 0), 
there is an extremum of the total energy: = 0. The analysis of the second derivative in 
the bounce depends on the location of X2 through the equations: 

p, = \'^ {x" - 16a^o') {^^ - 24a^o') • (215) 

At X = Xi the density is a minimum. For x = X2 the extremum depends on the location 
of the bounce respect to the point in which the quantity p + p changes sign. For the case 
in which VoaM^' < X2 < 2AaJffo^, it follows that the density has a maximum at Xi. On 
the other hand if X2 < IQaJ^o^ it is a minimum. To understand completely the behavior 
of the energy density the existence of other critical points for p must be addressed. This is 
controlled by equation (12091) . Thus, the extra extremum (which are not bounce or turning 
points) occur at x such that 



2 



P + P=-^{x'-lQaje^) = 0, (216) 

that is, at points in which the scale factor takes the value \/l6a^^^. Direct inspection 
shows that these are points of maximum density. 

Another consequence of nonlinear electromagnetism in cosmology is the occurrence of 
cyclic universe, as will be discussed in Sect J10.2l2l 
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4.6 Appendix 
4.6.1 Repulsive gravity 

A peculiar result which provides a framework to generate cosmological scenarios without 
singularity comes from the nonminimal interaction of EM with gravity, rendering gravity 
repulsive. The theory is defined by 

L = |i? - 1 F^.F'^^ + (3R A^A^^ , (217) 

where /3 is a dimensionless constant. This Lagrangian is not gauge-invariant and can be 
interpreted in terms of a photon having a mass which depends on the curvature of the 
geometry. 

Variation of gfj,u and yield the equations of motion: 

^+PA^^G,, = Pg,,nA^-[3A\^.^^-l3RA^A,-E^,, (218) 
FZ = -2/3RA'', (219) 



where E^^i, is Maxwell's energy-momentum tensor given by equation fll64p . As will be shown 
next, this set of equations allows a renormalization of the gravitational constant. Consider 
for instance the case in which A^A'^ = Z = constantT^ 0. Then 

(^i + pZ^ G^, = -pR A^A, - E^,. (220) 

Taking the trace of this equation we obtain i? = 0, and inserting this result back into 
Eqn.([220]) we get 

where the renormalized constant k is given by 

4 = i + 3Z. 

K K 

Thus, Eqns.( l2T8l) and ( 1219^ can be written as 

R^, = -kE^,, F'^^ = 0, (221) 

which are nothing but Maxwell's electrodynamics minimally coupled to gravity with a re- 
normalized gravitational coupling plus the condition A^A^ = constant = Z. 

The addition of other forms of neutral matter, such that the corresponding energy- 
momentum tensor is traceless, takes the Lagrangian to 



i^ = v^|^i?-^ F.^F^-' + /3RA,A^^ + L^'")! 



(222) 
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where L*^™^ represents the Lagrangian for all other kinds of matter such that Tj^f^ g'^'^ = 
T^™) = 0. The equations of motion in this case are given by 

+ P A^^ G^, = pnA'g^, - (3Al.^^ - PRA^A, - E^, - T"^^ (223) 

F'^^ = -2PRA''. (224) 

Taking again the case A^A^ = constant, yields R = 0. Then Eqns. (1223112241) take the reduced 
form 

F^.; = 0, 

where k was given above. Thus, the renormalization of the gravitational constant by the 
non-minimal coupling represented by the presence of the term RA^A^^ in the Lagrangian in 
the state where A^A^ is constant is still valid in the presence of matter with null trace. 

4.6.2 Global Dual invariance 

While observation must be the ultimate judge of the choice among the possible couplings, if 
it scarce or not available, we can resort to criteria coming from theoretical considerations. 
One of them is related to the invariance of the Lagrangian under a given transformation, 
such as the dual rotation. A dual map is a transformation on the set of the bi-tensors F^^ 
such that 

F^.u ^ Fl^ = cos OF^, + sin 9F;^. (225) 

Classical Maxwell's electrodynamics is invariant under such transformation only if the angle 9 
is constant. In a Minkowskian background it is not possible to implement such invariance for 
a local map 9 = 9{x). However, this can be achieved in the case of a non-minimal coupling of 
the electromagnetic field with the metric of a non-fiat geometry. In fact, using the identities 

F 

c 

it can be shown that the combined Lagrangian: 

1 f ^ 1 



^DI ~ -^5 — ^ 



Li = [r,u -\r 9,}j F^^F^r (226) 



is invariant under local dual rotations: Ldi = Ld\. This is a remarkable property which has 
no counterpart in the fiat space limit. 
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5 Viscosity 



A full knowledge of the global properties of the universe cannot be achieved without giving 
a description of the thermodynamics of the cosmic fluid. In the last decades, this task 
was addressed in three distinct periods. In the first period the universe was treated as a 
system in equilibrium in which all global processes were described by classical reversible 
thermodynamics, in such a way that total entropy was conserved. The salient feature of 
this phase was the development of the standard cosmological model, which comprises the 
homogeneous and isotropic FLRW geometry, and the characterization of the matter content 
of the universe as a one-component perfect fluid in equilibrium. In order to solve the EE, 
the energy density p and the pressure p were considered functions of the cosmological time 
only, and they were related by a linear EOS p = Xp. The FLRW models generated in this 
way share the common property of having an initial singularity (with A > —1/3). 

Later, it was realized [314] that the validity of thermal equilibrium near the initial sin- 
gularity is perhaps too strong an assumption. A second phase then started, in which the 
description of the cosmic fluid was improved by allowing viscous processes. Some of the 
motivations for this alteration are the following: 

• The examination of the possible role of viscosity in the dissipation of eventual primor- 
dial anisotropics (chaotic cosmologies), 

• The effect on the existence and/or the form of the singularity, 

• The application in cosmology of results obtained from non-equilibrium thermodynam- 
ics. 

In 1973 a FLRW cosmological model without singularity was presented [5lj (see also [55]). 
using a viscous fluid as a source. The energy-momentum tensor was given by 

T^,u = (p + p) v^v^-pg^^, 

in which p = pth — C^', where pth is the thermodynamical pressure, ^ is a viscous coefficient 
and 6 is the Hubble parameter, which is exactly the case of the energy-momentum tensor 
representing particle creation |315j . The SEC in this case is given by the inequalities 

P + Pth> 

and 

p + 3pth > 0, 

which are weaker than the correspondent ones in the case of perfect fluid, hence allowing for 
the absence of singularity. The solution found in |5l] is nonsingular, and past-eternal. 

More general forms for the dependence of viscous quantities have been investigated for 
arbitrary Stokesian regimes in which the fluid parameters become more general (for instance 
nonlinear) functions of the expansion. With these modiflcations, there are non-singular 
cosmological solutions, but they may suffer from a possibly worse disease than the initial 
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singularity: they are unstable and display non-causal propagation. In fact, the instability of 
the model in [51] was proven by the analysis made in [315j . It was also proved in |315j that 
the avoidance of the singularity is not generic. In other words, the singularity is not avoided 
for any type of viscosity (that is, for any dependence of the coefficients of viscosity on the 
expansion factor). 



In this second phase, local equilibrium [56] is still imposed, in such a way that the 

thermodynamical variables are described as if the dissipative fluxes - e.g. heat flux - do not 
influence local variables like for instance the entropy, although as a whole the system is not 
in equilibrium. As another example, a fluid in the regime 



was analyzed in |246j . both for a = (3 = constant, and a = 0, /5 = Mp"^, with M and m 
constants. In the second case, nonsingular solutions were found using tools from dynamical 
systems analysis. 

Let us remark that in general, the imposition of local equilibrium leads to causal difficul- 
ties, allowing dissipative signals to travel with infinite velocity of propagation. These causal 
problems were the focus of the third phase, where extended irreversible thermodynamics was 
used |247j . In this theory, the basic quantities become dependent not only on local variables 
of classical thermodynamics but also on the dissipative fluxes. This has very important 
consequences, the most important one being the preservation of causal connections for the 
whole system. In |248] . a FLRW universe was studied in this context, the net consequence 
of the assumption of extended irreversible thermodynamics being to provide an additional 
equation of motion for the non-equilibrium pressure tt, with p = pth + tt, given by 



(where tq is the relaxation time) which preserves the causal structure. Thus, contrary to the 
previous case in which the viscous term is assumed to be a polynomial on 6, here it must 
obey Eqn. fl227p . The other quantities relevant to thermodynamics (that is, the entropy flux 
and the particle flux per unit of proper volume n) are determined by 



Assuming an EOS given by pth = Ap, the cases ^ = constant, and ^ = /3p, (with /3 = 
constant) were analyzed in [248] . always with tq = constant, and nonsingular solutions were 
discovered in both cases, for A = and A = 1/3. The relevant equations of this system can 
be put in the form of an autonomous planar system: 




p = pj^a9 + pe"^ 



tqti + tt = -^6. 



(227) 
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where A is the cosmological constant. The set of integral curves of this system was studied 
in 12491. where it was shown that the solution found in 12481 is stable. 



Bifurcations in the early cosmos 

Quadratic dissipative processes were analyzed from a new perspective in j317j . where it 
was shown that dissipative processes may lead to the appearance of bifurcations. This is a 
consequence of the application of a theorem due to Bendixson |318] to the system of EE that 
describes a universe with curvature controlled by a dissipative fluid. Indeed, let us consider 
a planar autonomous system that contains a parameter, say a, of the form 

X = F{x,y;a) 

y = Gix,y;a), (229) 

where the functions F and G are non-linear and the parameter a has a domain V. Applying 
methods of qualitative analysis to this system and restricting to the two-dimensional plane 
r of all integrals of this system, one arrives to the notion of "elliptical" and "hyperbolic" 
sectors, that characterize, as the names indicates, the behavior of the integral curves in the 
neighborhood of a multiple equilibrium point (that is, an isolated points that is a zero of 
both F and G). Let us call S and Ti. the number of elliptical and hyperbolic sectors of a 
given equilibrium point M = {xq, y^) of F, respectively. Then the Poincare index is defined 
by the formula 

£-n , 

Ip = h 1. 

2 

This is a measure of the topological properties of the integral curves in the phase plane F. 
If above a certain value of V the topological properties of the system (12291) change, then 
there is an abrupt change of behavior of the physical system in the vicinity of the unstable 
equilibrium point. The crucial consequence of the above-given theorem is the appearance 
of indeterministic features. In j317] this theorem was applied to spatially homogeneous and 
isotropic cosmological models, whose dynamics is described by a planar autonomous system, 
given by 

p = -^pe + ae^ + (3e\ 

where a (referred to in the theorem) can be either a, (3 or 7, and the energy-momentum 
tensor is 

T^,y = (p + p) v„ - pg^^, 

where 

p = Pth + aO + (3e'^, 

with pth = (7 - 1)P- 

The viscous terms (parameterized by a and /3) can be a phenomenological description 
of particle creation in a nonstationary gravitational field as proposed in [31 9] and [320J. 
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Applying the methods of quahtative analysis to the system given in Eqn. fl230p it was shown 
in [317] that for 7 — 3/? < 0, the Poincare index Ip{B) = —1 (saddle point); for 7 — 3/5 > 0, 
Ip{B) = 1 (two-tangent node). This situation characterizes a bifurcation in the singular 
point, when p = 9 = 00. This bifurcation, caused by dissipative processes involving quadratic 
viscous terms generates a high degree of indeterminacy in the development of the solution 
of EE, which enshrouds the past of this model of the universe. In this case, nothing can be 
stated about the existence of the initial cosmological singularity. 

6 Bounces in the braneworld 

Theoretical developments coming from string theory have revived the idea that our universe 
may have more than 4 dimensions (first considerated by Kaluza in the context of unification 
of gravity and electromagnetism) . Among the multidimensional models, those with one or 
more branes that live in a bulk space have been thoroughly studied recently (see for instance 
|330j ). In these models, the matter fields are typically confined to a 3-brane in 1 + 3 + c? 
dimensions, while the gravitational field can propagate also in the d extra dimensions, which 
need not be small, or even finite, as shown in one of the models introduced by Randall and 
Sundrum where for d = 1, gravity can be localized on a single 3-brane even when the 

fifth dimension is infinite. The Friedmann equation on the brane is modified by high-energy 
matter terms and also by a term which incorporates the nonlocal effects of the bulk onto the 



where e is the 3-curvature, H = d/a, p is the energy density of the matter on the brane, 
= Sir /M M PI is the fundamental 5-dimensional Planck mass, = Sn/Mpi, and 



where A is the tension of the brane, and A is the 5-dimensional cosmological constant. Finally, 
Uq is the constant corresponding to the non-local energy conservation equation. This term 
comes from the projection of the Weyl tensor of the bulk on the brane |33U] . From Eqn. (12311) 
we see that a necessary condition to have a bounce with p > in the e = 0, — 1 cases is that 
either A < or W < 0, or both. 

The case that includes matter in the bulk, without cosmological constant for a flat FLRW 
d + 1-dimensional was studied in [113] . A neceessary condition in order to have a bounce is 
that dH/dt > 0, with 



brane [TT2l[33n] : 



^ =3 + y^-^ + 36^ nuJ "HJ ' 



(231) 
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where (in a notation slightly different from that used in Eqn. fl23ip ) k is the bulk gravitational 
coupling, Gn the effective Newton constant on the {d + 1) -dimensional brane, E is the 
projection of the bulk brane Weyl tensor on the brane, and = (— -R, P) is the projection 
of the bulk energy-momentum tensor on the brane. It follows from this equation that a 
necessary condition to have a bounce without resorting to exotic forms of matter (that is, 
matter that violates p > or pp > 0) is a negative Eq |113j . This is precisely the approach 
taken in [2761 1114j . where a brane evolving in a charged AdS black hole background was 
studied. Bouncing solutions were found for both critical (A = 0) and non-critical (A ^ 0) 
branes, the bounce generically depending on the parameters of the black hole, and on the 
matter content of the brane PI . 

The abovementioned necessary condition was explicitly checked in the case of the dilaton- 
gravity braneworld [113] , and bouncing solutions were obtained for a a fiat FLRW brane in 
a static spherically symmetric bulk0. This solution describes (in the string frame) a pre-big 
bang model where the transition between the branches is realized at low curvature and weak 
coupling, thus providing an example of succesful graceful exit without resorting to quantum 
or "stringy" corrections. 

Notice that the extra dimension(s) could be spacelike or timelike. The latter case was 
analyzed in [115j . The usual incantations [330] for the case of an extra timelike dimension 
and an homogeneous and isotropic brane lead to [115] 

e A SvrGp p^ C , , 

*^^=3+V^-4^^' ''''' 

where G and A are the effective gravitational and cosmological constant, respectively, and 
M is the 5-dimensional Planck mass. Notice that the minus sign in front of p^ may lead 
to a bounce instead of a singularity, since this term grows faster than the others, leading 
to if = 0, this feature being independent of the equation of state and also of the spatial 
curvature of the universe. The simplest of these bouncing universes, described by 



"=^'-k ''''' 



will be discussed in Sec J10.2.Tl since it may lead to a cyclic universe. 

The case with an extra timelike dimension in this scenario was also extended to Bianchi I 
universes [115j . which exhibit an anisotropic bounce as long as the shear scalar aapcr"'^ does 
not grow faster than as a goes to zero at the end of the contraction phase. All these 
results were obtained by neglecting the induced curvature on the brane, which can trigger 
the formation of a singularity at the beginning or at the end of the evolution [115] . 

Another model along these lines was introduced in [392] . where a "test brane" {i.e. one 
that does not modify the ambient geometry) moves in a higher- dimensional gravitational 
background. Using the thin-shell formalism, in which the field equations are re-written as 



^^The bounce in the model presented in [114j was analyzed from the point of view of the causal entropy 
bound in |417j . and its stability was put in doubt in [418j . 

^^Bouncing solutions for a domain wall in the presence of a Liouvillc potential were found in [329] . 
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junction conditions relating the discontinuity in the brane extrinsic curvature to its vacuum 
energy, the motion of domain walls in de Sitter and anti-de Sitter (AdS) time-dependent 
bulks was discussed. This motion induces a dynamical law for the brane scale factor, and 
it was shown in |392j that in the case of a clean brane the scale factor may describe a 
non-singular universes. In order to build the class of geometries of interest, two copies of 
(d+l)-dimensional dS (AdS) spaces Ml and M2 undergoing expansion were considered. 
From each of them, one identical d-dimensional region fij (z = 1,2) was removed, yielding 
two geodesically incomplete manifolds with boundaries given by the hypersurfaces dQi and 
8^2- Finally, the boundaries were identified up to an homeomorphism h : dQi dQ2- 
Hence, the resulting manifold that is defined by the connected sum M1#M2 is geodesically 
complete. The starting point is the action 

S=^ d'+'x^{R-2A) + ^ d^x^lC + a d^x ^, 
16vr Jm Stt Jq^ Jq^ 

where the first term is the usual Einstein-Hilbert action with a cosmological constant A, the 
second term is the Gibbons-Hawking boundary term, K,mn is the extrinsic curvature, and 
a is the intrinsic tension of the ci-dimensional brane. The spatial coordinates on dVL can 
be taken to be the angular variables 0j, which for a spherically symmetric configuration are 
always well defined up to an overall rotation. Generically, the line element of each patch can 
be written as 

ds^ = -de + A\t)[r^dnl^^^) + (1 - kr^)-^dr% 

where e takes the values 1 (-1) for dS (AdS), is the corresponding metric on the 

unit d — 1-dimensional sphere, and t is the proper time of a clock measured in the higher- 
dimensional spacetime. In order to analyze the dynamics of the system, the brane is allowed 
to move radially. Let the position of the brane be described by a;^(r, 0i) = (t(r), a(r), 0j), 
with r the proper time (as measured by co-moving observers on the brane) that parameterizes 
the motion, and the velocity of a piece of stress-energy at the brane satisfying u^um = — 1- 
With these assumptions the brane will have an effective scale factor A'^{t) = a'^{t)A'^{t). The 
constraint 

dr^ f^^Mj^ 



dt V 1 — (-0? 

along with the result of the integration of EE across the boundary (done with the junction 
conditions) |392j yields two differential equations for A and a. For the case of a background 



composed by two patches of dS undergoing expansion, A{t) = i cosh.{t / i) , and e = 1, where 
£2 — d{d — 1)/|A| is the dS radius. In this case the EOM for the brane is 



Att ±dsinh(t/£) + [ai cosh{t / - a 

o — 



2a 



4'-"V - 1) (1 - a2 - [£dcosh(t/£)]2)i/2 

Nonsingular analytical solutions of this equation for cr = can be obtained, while for cr 7^ 0, 
numerical methods must be used. This latter case also yields bouncing solutions (see Fig 16]). 



69 




Figure 9: Effective scale factor A{t) plotted for a = A'ira/{Lp \d — 1)) = 10 ^. The insert 
displays the dependence of A with a. Taken from [392] . 



The motion of a test brane in a background produced by a collection of branes was 
discussed in |331] (the so-called mirage cosmology). Adopting spherically-symmetric back- 
grounds, it was shown that although there is a singularity in the evolution of the 4-d brane, 
the higher- dimensional geometry is regular. The origin of the singularity on the brane is ac- 
tually the embedding of the brane in the bulk, in such a way that the singularity is smoothed 
out when the solution is lifted to higher dimensions. 

The effect of inflation on a bouncing brane was used in |408j to set limits on the parameters 
of the braneworld. Specifically, the model consists of closed a FLRW metric embedded in 
a 5-d conformally fiat bulk with one extra timelike dimension, containing a conformally 
coupled scalar field (the infiaton field) and a radiation fluid, evolving on the brane with 
corrections due to the bulk. The non-singular bouncing solutions considered were oscillatory 
and bounded, or initially bounded. They are in principle stable and would never enter an 
inflationary phase with an exponential growth of the scale factor since they correspond to 
periodic orbits of the integrable dynamics in the gravitational sector. The introduction of a 
massive scalar field, even in the form of small fluctuations, turns non-integrable the dynamics 
of the system |408j . As a consequence, non-linear resonance phenomena are present in the 
phase space dynamics for certain domains of the parameter space of the models, and the 
associated dynamical configurations become metastable, allowing the orbits escape to the de 
Sitter infinity in a finite time. From the conditions for these orbits to happen, limits on the 
parameters (a, m, Eq) are set, where a is the brane tension, m is the mass of the scalar field, 
and Eq is a constant proportional to the total energy of the fluid. 
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Yet another turn in the mirage model was introduced in |251j . where the brane moves 
in an open orbit around a non-trivial spherically-symmetric background. In this model, the 
brane is moving on a Calabi-Yau manifold generated by a heap of D3-branes, and the mirage 
effects dominate the evolution of the Universe only at early time, i.e. when the brane moves 
in the throat of the background manifold. The new feature is the influence of the angular 
momentum of the test brane on its motion in the higher-dimensional space. In fact, the 
effective 4-d metric has two parameters: the energy U and the angular momentum L of the 
4-d brane, which determine the form of the orbit. In particular, to have an open orbit in an 
asymptotically Minkowskian background, 

L^_4(f/ + 2)f/3 > 0. 

As discussed in |251j . the effective metric corresponding to orbits satisfying this constraint 
display cosmological contraction during the ingoing part of the orbit, expansion during the 
outgoing part, and a bounce at the turning point l3- 

Another model based on the brane scenario is the ekpyrotic universe |105j . the cyclic 
version of which shall be considered in SectlTOl 



7 Variable cosmological constant 

General Relativity allows for the introduction of only one arbitrary constant, the so-called 
cosmological constant A. At least two attitudes can be taken regarding A [292]. The first 
one is to consider it as a derived quantity, that emerges from vacuum fluctuations (see for 
instance [293j). One way out of the huge disagreement between theory and observation in 
this case |138] is to assume that A is actually time-dependent. The second attitude that 
can be adopted is that A is, along with G, a fundamental parameter of the theory, to be 
determined by observation PI |294j . In fact, from a gravitational point of view what matters is 
the "effective" cosmological constant, since the matter Lagrangian can sometimes contribute 
with a A-like term, as in the case of the scalar field in the presence of a potential with a 
minimum: 

Aeff = A + 

where A is the "bare" cosmological constant. Any change in 0min during the evolution 
leads to changes in the value of Aefr. In fact, the effect of the evolution of the universe on 
the ground state is to add a temperature dependence, which can be translated into a time 
dependence |253j . A model along these lines based on a gauge field (instead of a scalar 
field) was presented in [262J 0. This is another motivation to consider a variable A, that 
is not a constant but a function of spacetime coordinates, in such a way that its value is 



^^Further effects of the angular momentum on the motion of the brane, including cyclic universes, were 
studied in [382]. 

^^Notice that this second attitude is somewhat different from Einstein's original ideas leading to GR, since 
there would be curvature even in the absence of matter, caused by A. 

^^In fact, any classical nonlinear field theory (such as nonlinear electromagnetism) admits a fundamental 
state that generates a cosmological constant [326 . 
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determined by the dynamics of the theory under scrutiny (following the line of reasoning of 
other "variable constant" theories, see Section (TDJ). In fact, a time-dependent cosmological 
constant has also been called upon to explain the current accelerated expansion and the fact 
that this phase started in the recent past. 

In the case of A = A{t), EE for the FLRW metric take the form 

1 A(t) e 
d A{t) 1 



a 3 6 
and the continuity equation is given by 



(p + 3p), (236) 



p + 3- (p + p) = -A. (237) 
a 

As seen from Eqn.f l237|l . A can supply or absorb energy from ordinary matter and radiation. 
In fact, it follows from this equation that 

TdS = -VdA. (238) 

Hence, A is a source of entropy. Requiring that dS/dt > implies dA/da < through cosmic 
expansion. 

Assuming that only radiation is present, Eqn. fl237p gives 

dp dA 4p _ ^ 
da da a ' 

which can be integrated to 

where p = po when a = oq, and the subindex denotes quantities evaluated at t = 0. Notice 
that the model is completely determined in this case by providing the function A = A(a), 
since Eqn. fl239p then yields p = p(a), and a = a{t) follows from Eqn. fl235p . A cosmological 
model based on this scenario was discussed in [252j . where the dependence of A on a was 
fixed by imposing that p = Pc for all values of t, where Pc = 3if^ is the critical density. It 
follows from Eqn. (12351) that 

A = ^. (240) 

The conditions A > and d > give e > 0, hence e = 1. In the model presented in [252J, at 
t = the universe had only a nonzero cosmological constant. With po = 0, Eqn. fl239f240p 
give 

P(a) = ^ fl - 4V (241) 



a^ V a^ 
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Note that po = implies that Qq 7^ 0, in such a way that the singularity at t = is absent. 
An estimation of Oq was made in j252j by assuming that the maximum temperature reached 
is Tmax ~ Mpi, which gives 

ao ~ 4i X 10-'°(GeV)-\ 



where N = N{T) is the effective number of degrees of freedom at temperature T. 

The fact that this model does not display a horizon problem was also shown in [252] . In 
fact, the time at which global causality is established is given by 

tc = ttQ sinh — ~ 2.3ao, 

which indicates that global causal connection was established at a very early time. The 
model is also free of the monopole problem, but it is worth noting that there is an inflationary 
period. From Eqn. fl235p we get 

a'^ = al + t'^. (242) 

A peculiarity of this model is that a — cx) for t — >• oo, even though e = 1. Needless to say, 
other choices of A would give a different asymptotic behavior. 
The same form of A, namely 

m = (243) 

where 7 is a constant to be determined by observations, was studied in [238], but without 
the assumption that p = pc- The conservation equation (123 7p can be solved for dust and 
radiation. Inserting the solution in Eqns. (12351) and (12361) we get 



- + - = 

a? 3 



a 6 

where T = e — 27/3 for radiation, and T = e — 7 for dust, and p'-*-' is the energy density 
of dust or radiation for the case A = 0. These equations show that the effect of assuming 
that A oc is to shift the curvature parameter e by a constant value. A nonsingular 
cosmological model based on the model presented in [238j has been analyzed in |241] . Notice 
that Eqn. (l243p along with condition dA/da < require that 7 be positive. A positive A for 
all t implies, through Eqn. (12361) that there may be a zero in a, and hence the possibility of a 
bounce. For this to happen we need that d be zero at the putative bounce. Supposing there 
is a bounce, it follows from Eqn. (I235p evaluated at the bounce that 

«~Voao = 6-7. 

Hence, po > implies that e > 7 > 0, and so e = 1. Introducing the Ansatz (I243p in the 
Friedmann equation, we get 

aV = (27-l)(a2-a^), (244) 
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so it follows that 7 > 1/2. Hence, 1/2 < 7 < 1. This equation can be integrated to get 

= (27 - l)^ + al 



which leads to bounded-from-above densities and temperature Ifj. 
Yet another form for the dependence of A, given by 

A = Al + A2 a-*", 

where Ai, A2 and m are constants (with A2 > 0), was studied in |295j . The analysis of the 
dynamics was carried out using the analog of the one-dimensional problem of the particle 
under the influence of the potential V{a) given by 

V{a) = -AiS^- - A2S — + &a-", 

a + 2 a — m + 2 

where a = 1 + 3A, 5 = 1 + A, 6 is a positive integration constant, and p = Xp. Denoting by 
r the maximum of the potential, cyclic solutions are obtained for the cases e = 1 with Ai, 
A2 > 0, and r > —1, and for Ai < 0, A2 > 0, and m < 2, regardless of the sign of e. 
The proposal in Eqn. (12431) was later generalized in |239j to 

A = 3(3H^ + ^, (245) 

where (3 and 7 are dimensionless numbers, and H = a/a. following |254j . With this Ansatz, 
the Friedmann equation for a radiation-dominated phase can be rewritten as 

a' = ^ + Aoa-'^'^, (246) 

which allows a bouncing solution at t = for < 0, /? < 1/2, e = 1 (with po > 0). The 
value 7 = 1 was chosen in [239] so that dS/da is always greater than zero, thus solving the 
entropy problem. In this case, the model gives < 1 for all t. 

A thorough review of variable-A models has been presented in [296] . The models analyzed 
were power-laws of the different relevant parameters, namely 



Al = Ar^, A2 = Ba-"", A3 = Ci/", A4 = Vq 



where A, B, C, V, i, m, n, and r are constants. Let us state from [296] the relevant 
results for this review: (1) no bouncing models were found for Ai with k = and £ = 
1,2,3,4, irrespectively of the sign of A. (2) For A2, it was shown (numerically) that there 
are nonsingular models for dust, e = 1, with m = 1, Qq = 0.34, and 0.68 < Qqa < 0.72, and 
also with higher values of m and Qq. (3) For A3, the value n = 2 admit analytical solution. 
For this n, there are bouncing solutions for 7 > 2/3 and e = 1 with C > 3(37/2 — l)Qo, and 
also for 7 > 2/3 and e = —1, for C < 3(37/2 — l)Qo. (4) Only the value r = 1 was explored 
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The evolution of perturbations in this model was studied in [240] . 



74 



for A4. Defining Aq = —'DqQ/3, there are closed bouncing solutions for Aq > —^0, and open 
bouncing solutions for Aq < — ^^o- 

The examples given above show that varying-A scenarios are worth examining because 
they address a number of pressing problems in cosmology (horizon problem, entropy, initial 
singularity)!^. Furthermore, many of them are simple enough to draw definite conclusions 
about their viability. One of the drawbacks is perhaps the lack of strong motivation for choos- 
ing any given form of A. In this regard, let us remember that many of the varying-A models 
can be reverse-engineered to scalar-field models with a potential. Unfortunately, in most 
cases the corresponding models lack predictive power or clear particle physics motivation 
[292]. 

8 Past-eternal universes 

In this section, we shall examine some models which are nonsingular but do not exhibit 
a bounce. Historically, perhaps the most important example of these is the Steady-State 
model |264j As mentioned in Sect 12. ![ nonsingular solutions that start from a deSitter 
state were discussed in [39| HOj. Another example is that discused in |6T] in which every 
contracting and spatially flat, isotropic universe avoids the big crunch by ending up in a 
deSitter state enforced by the hmiting curvature hypothesis. 

8.1 Variable cosmological constant 

As noted in [255] . in all the articles mentioned in Sectd the dependence of A on a and a 
was set either from "first principles" (for instance quantum gravity, as in [238] ). or by ex- 
trapolating backwards current cosmological data, including the current value of A. However, 
another view can be taken. Since A can be considered as a remnant of a period of infiation, 
a complete model should also describe the era of infiationary expansion. This is precisely 
the proposal in |255] , where A was taken as 

A(H) = 3PH' + 3(1-/3)^, (247) 

where Hi is the timescale of infiation, and /5 is a parameter. Note that when H = H^, 
A = 3H/, as required by infiation, while A ~ 3(3H^ for large cosmological times. In the case 
of e = 0, and for 

P = (7 - 

an equation for the Hubble parameter follows [255] : 



^^Nonsingular cosmological solutions for the case in which the cosmological constant is replaced by a 
second-rank tensor A'^ were studied in [242j . 
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For an updated version, see Sect ll0.2T5l 
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whose solution is 




where C is a 7-dependent integration constant |^. This equation can be integrated to yield 



where a* is an arbitrary value of the scale factor. It follows from this equation that the 
evolution of the universe starts from a deSitter stage a ~ c^** for Ca^'^^^"'^^/^ << 1, and 
evolves towards a FLRW phase, a ~ t2/37(i-/3) for Ca^^'^^'^^/^ » 1. 

8.2 Fundamental state for f{R) theories 

A novelty in some theories described by Lagrangians that depend only on R is the possibility 
of the emergence of an intrinsic cosmological constant. This is not the case, however, in 
theories generated by Lagrangians that are a linear combination of and R^^R^^ as can be 
seen by a direct inspection of the EOM (|35l) . The proof of this assertion follows from the fact 
that the tensors Xfj.v ^^id Z^^, appearing in the EOM [27] are traceless in the case of a constant 
curvature scalar. {R^u = ^dfiu)- However, restricting to the f{R) case, Lagrangians that 
are not linear in R^ can bypass such prohibition. The existence of a deSitter solution in the 
absence of matter occurs when the function obeys the condition 



It follows straightforwardly from Eqn.(reffr) that R^^, = hg^^ is a possible state of the 
system. 

8.3 The emergent universe 

Another example of past eternal universe was given in [336]. This model uses general rela- 
tivity plus a scalar field with a potential, and matter. The relevant equations are 



— = constant. 



(248) 



A typical example is provided by the exponential Lagrangian 




+ 3if0 + = 0, 



a 



a 



Here the value e = was chosen, but this restriction was lifted in [263] . 
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From these, it follows that to have a minimum of the scale factor we need to impose the 
conditions 

^(l + 3u;)p + 02< 

and 

-0i +Vi + Pi = 



where the subindex i means that the quantities are evaluated at ti, the time at which a is 
minimum. Assuming positive potentials and energy density, it follows that only e = +1 is 
allowed. It follows that 

1 2 

-(1 - uji)pi + Vi = —, 



where = Aj, and 



2 



so a model can be constructed with pi = and constant 0^. This can be achieved in the 
limit t oo with the potential |336] 



V{(P) = Vf + {V, - Vf) 



exp 



a 



where 0/ is the value of the field for which V is minimum, and a is a constant energy scale. 
In order to achieve the Einstein universe state in the far past, some fine-tunning on and (pi 
is needed, which is not necessarily a hindrance }336j^ . The choice of such a highly-symmetric 
state as the initial state is supported by various arguments: it is stable agains some types 
of inhomogeneous linear perturbations, it has no horizon problem, it maximizes the entropy 
within the family of FLRW radiation models, and it is the unique highest symmetry non- 
empty FLRW model (with a 7 dimensional group of isometries). The model was elaborated 
further in [337] . where it was shown that an explicit form for the potential can be found 
such that the model leaves the inflationary stage and enters a reheating phase, followed by 
standard evolution. 



9 Quantum Cosmology 

As discussed in SectlH there are reasons to suppose that at very high energies some of the 
hypotheses of the singularity theorems are rendered invalid: if the universe ever attains this 
regime, an important role is to be played by quantum gravitational effects, in such a way 
that a quantum theory of gravitation is needed to have a proper description. 

^^In particular, the initial scale factor could be chosen in such a way to avoid the quantum gravity regime. 
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Although there is yet no complete realization of quantum gravity, there are some attempts 
to tackle the singularity problem in a quantum framework. A standard method of quantizing 
General Relativity is canonical quantization [299] where the momentum and Hamiltonian 
constraint equations are interpreted as operators, and it is required that they annihilate the 
quantum state. The Hamiltonian constraint gives the Wheeler-DeWitt (WdW) equation 
|267j . which depends on the choice of the factor ordering in the products of generalized 
momenta and "velocities". For some choices of the ordering, the WdW equation turns it 
into a Klein- Gordon equation on an indefinite DeWitt metric in the infinite-dimensional 
superspace (space of three-metrics), with a potential term |267j . In addition to the WdW 
equation, initial conditions must be specified, the two most popular being the "no-boundary" 
[338], and the "tunnelling" condition [339j . 

In practice, the infinite degrees of freedom of the superspace are truncated to obtain a 
minisuperspace model, usually under the assumptions of isotropy and homogeneity. Once a 
solution to the WdW equation has been found, there is the question of how to interpret it 
and extract probabilities from it. 

Among other issues related to the WdW equation, there is the fact that a suitable initial 
condition must be chosen to get a solution. It would be desirable that the initial condition 
be somehow determined by the dynamical law (see for instance |159j ). In fact, the most well- 
accepted proposals mentioned above do not solve the singularity problem [164] . Moreover, 
in the quantization following the ADM procedure, time is fixed by a gauge choice, and the 
results are dependent of this choice |271] 

As we shall see below, there are other approaches to Quantum Cosmology which may 
yield a nonsingular universe in the regime where the WdW equation is valid. We shall discuss 
two possibilities: the Bohm-de Broglie interpretation of QM, and Loop Quantum Cosmology 
(LQC). 



9.1 The ontological (Bohm-de BrogUe) interpretation 

If the universality of quantum mechanics is assumed, the Universe must be describable 
by a wave function (furnished by a yet-to-be-discovered quantum theory of gravity and 
matter fields) in every step of its evolution. Moreover, this description must have a well- 
defined classical limit. The orthodox interpretation of Quantum Mechanics (the so-called 
Copenhagen interpretation) [16U] is ill-suited for the task of describing the universe, since it 
assumes the existence of a "classical apparatus" external to the system to solve the measure 
problem by forcing the collapse of the wave function. Clearly, there is no classical apparatus 
outside the universe. Therefore, the least we can say is that an alternative to the Copenhagen 
interpretation is needed. One such alternative that has received some attention recently 
is that of Bohm and de Broglie (BdB)[158j 0- In classical physics, the dynamics of a 
point in configuration space is determined by the principle of extremal action, yielding the 

^■^In this regard, it was shown in [300] that a Bianchi I universe, quantized following the ADM recipe with 
a particular choice of the time coordinate 161J in the presence of dust is nonsingular. 

^''Other possibilities (not free of problems, though) are the many- worlds interpretation [304], non-linear 
quantum mechanics |305j , and decoherence [306j . 
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classical EOM. According to the BdB interpretation, in quantum physics the evolution of 
the configuration variables is guided by a quantum wave which obeys Schrodinger's equation. 
The associated Hamilton- Jacobi equation displays a new term (of quantum origin, see below), 
that can be interpreted as part of the potential. It should be emphasized that the BdB 
interpretation furnishes a framework to make predictions based on the wave function of the 
system, which must be obtained by some means (for instance, through the WdW equation). 

Let us briefiy review first the quantum mechanics of a single particle in the BdB inter- 
pretation, and afterwards the results will be translated, mutatis mutandis, to the context of 
FLRW cosmology. The Schrodinger equation for a non-relativistic particle in a potential V 
is given by 

at \ 2m J 
With the replacement = Rexp(iS/h), this equation becomes 

^^vYfl^^Uo. (250) 



dt \ m 

This last equation suggests that VS/ m can be interpreted as a velocity field, leading to the 
identification p = VS, in such a way that Eqn. fl249p is the Hamilton- Jacobi equation for the 
particle in the classical potential V plus a "quantum potential" Q = —h^'V^R/2mR. The 
BdB interpretation argues that a quantum system is composed of a particle and a field, and 
that quantum particles follow trajectories x{t), independent on the existence of an ouside 
observer. These trajectories can be determined from 

or fromp = mx = VS*, after 5* and R are determined using Eqns. fl249p and fl250p . In practice, 
since S is the phase of the wave function, it can be read off from the explicit solution of 
Schrodinger's equation. 

Let us analyze an example developed in [272] . where the Lagrangian was given by 

C = V^ {R-C^ct>,,<P'^), 

where C^; = (cu + |). From the metric 

ds^ = -N'dt^ + "f^l' , (dr^ + r^dn') , 

and the definitions /J^ = Airfpi/^V, <p = (j)y/CjQ, we get 



n = N\ -p-'^ + pA-e 



2a ' " 2a3 2(3^ 
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with Pa = —aa/{P'^N), = a^(f)/{P'^N). Defining a = a/P, setting /3 = 1 and a = Ina, we 
get 



where 

Pa 



AT ' N 

Notice that p,^ = A; is a constant of the motion. We shall restrict to the case e = since it is 
analytically tractable. The classical solutions are given by 

a = 3kf^^, (f)=-\nt + C2, 
3 

where C2 is an integration constant. Depending on the sign of k, this solution contracts to 
or expands from a singularity. 

The Wheeler-DeWitt equation corresponding to the Hamiltonian given in Eqn. (12511) is 
given by |272j 

da'^ dcfy^ 

The solution, obtained by separation of variables, reads 

^(«,0) = j F{K)A^{a)B^{cf>)dK, 

where k is a separation constant, F{k,) is an arbitrary function of k, 

Af^{a) = ai exp(ma;) + 02 exp(— ma), 

(for e = 0), and 

Bk{(P) = 6iexp(z/t0) + 62exp(— 

A direct application of the formalism sketched for the case of a one-particle system, gener- 
alized to several degrees of freedom yields from the Hamiltonian (12511) [27'2j 

Q{a,(p 



2R V 

with the "guidance relations" 

as e3"d dS e^" 



da N ' dd) N 



A state is now needed to read off from it 5* and R. A Gaussian superposition was chosen in 
[272] , given by 

^(a,0) = j F^B,{(t))[A,{a) + A^,{a)] dn, 
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with 

F{k) = exp ' ^ ^ 



a2 



and 02 = 62 = 0. Performing the integration in k, we can extract from the result the phase 
S which, when inserted into the guidance relations (in the N = 1 gauge) furnishes a planar 
system: 

00"^ sin(2rfa) + 2(isinh(cr^Q;</)) (252) 
exp 3a (2(cos(2(ia) + cosh(cr^a0))) ' 

— acr^ sin(2(ia) + 2(icos(2(ia) + 2(icosh(cr^a0) 



exp 3a (2(cos(2(ia) + cosh((T^a0))) 

The plot of this system (see Fig lTOl) shows that there are bouncing trajectories for a > 0, 
and also oscillating universes near the centre points (white points in the plot). The BdB 




Figure 10: Field plot of the planar system (12521) - fl253p for a = d = 1. Taken from |272j . 



interpretation has been applied to mini-superspace models in Quantum Cosmology (see for 
instance [26911271] ). and non-singular solutions have been found for models with scalar fields 
or radiation [407] . The bounce is due to the action of the quantum potential, which generates 
a repulsive "quantum force", large enough to reverse the collapse. 

One of the advantages of this formulation is that, starting from WdW equation, it yields a 
dynamics that is invariant under time re-parameterizations. Notice however that the results 
are dependent on the state chosen to represent the system. 

9.2 Loop Quantum Gravity 

Loop Quantum Gravity is a background- independent, non-perturbative canonical quanti- 
zation of gravity in which the classical metric and the extrinsic curvature are turned into 
operators on a Hilbert space [165] . The classical description of space-time is replaced by 
a quantum counterpart, in such a way that quantum effects are important at very short 



81 



scales, for instance near putative singularities. In this scenario, the evolution of the universe 
is divided in three epochs. First there is a quantum epoch with high curvature and energy, 
described by difference equations for the wave function of the universe. These are a direct 
consequence of the discreteness of space and time, the step size being dictated by the lowest 
non-zero eigenvalue of the area operator (see [159] ). It is this discreteness that modifies the 
behavior near the singularity, leading to a theory that is not equivalent to the WdW descrip- 
tion (even in the isotropic case), which furnishes a continuous spectrum for the scale factor. 
A semiclassical epoch follows, with differential equations for matter and geometry modified 
by non-perturbative quantization effects. Finally, a classical phase is reached, described by 
the usual cosmological equations. 

Since difference equations are often difficult to analyze or to solve explicitly, and at such a 
fundamental level, the emergence of space-time in inhomogeneous models with many degrees 
of freedom from the underlying quantum state is hard to understand, a suitable strategy is 
to use special models allowing exact solutions. Care must be taken in the extension of results 
from particular examples to more general cases. In any case, it may be instructive to have 
a detailed understanding of how the singularity is resolved in some instances. 

Yet another convenient simplification is to work in an effective semiclassical theory, which 
takes into account only some quantum effects. This theory can be understood as governing 
the motion of a wave packet that solves the difference equation [132] , and can be obtained as 
an asymptotic series of correction terms to the equations of motion in the isotropic case [131j . 
For instance, in the case of a matter term generated by a scalar field under the influence of 
a potential, the effective Klein-Gordon equation is [163j 



D{q) = - q'/\l[{q + 1)"/^ - |g - Ij^/^ - llg[(g + 1)7/4} -\q- l\'/Sgn{q - 



with q = a? /a1 and a1 = 7£pj73, where 7 ~ 0.13, and j is a quantization parameter, 
which takes half-integer values. This equation represents an approximate expression for the 
eigenvalues of the inverse volume operator |152j . The function D varies as a^^ for a <^ a^, 
has a global maximum at a ~ a*, and falls monotonically to D = 1 for a > a^,. In turn, the 
effective Friedmann equation is given by 




(254) 



where 





(255) 



and the effective Raychaudhuri equation is 




(256) 
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These approximations are valid for Oj < a < a*, where Oj = ipi- Below a, the quantum 
nature of spacetime cannot be replaced by an effective theory, while above a* we recover 
classical cosmology. It was shown in |163j that a closed universe with a minimally coupled 
scalar field will bounce (avoiding the so-called big crunch) as soon as a ~ a* for any choice of 
the initial conditions. The bounce in this case is due to the change of sign of the "friction" 
term in Eqn. (l254p . which becomes frictional for a << a*, freezing the field in some constant 
value, and turning the effective EOS into a cosmological constant EOS [163] . Similar results 
were obtained in the case of anisotropic models [133j. 

The previous example incorporated quantum gravitational effects on the matter (repre- 
sented by a scalar field) Hamiltonian, but there may also be modifications of the gravitational 
Hamiltonian due to quantum geometry. Recently, some calculations illustrating the effects 
of quantum geometry on both the gravitational and matter Hamiltonians were carried out 
in the case of a spatially homogeneous, isotropic e = universe with a massless scalar field 
(a system which is singular both classically and according to the WdW formalism in the 
Copenhagen interpretation of QM). It was shown in [164] that the singularity is resolved in 
the sense that a complete set of Dirac observables on the physical Hilbert space remains well- 
defined throughout the evolution; the big-bang is replaced by a big-bounce in the quantum 
theory due to the quantum corrections to the geometry; there is a large classical universe 
on the "other side", and the evolution bridging the two classical branches is determinis- 
tic, thanks to the background independence and non-perturbative methods 0. Notice also 
that no boundary condition was imposed (it was asked instead that the quantum state be 
semiclassical at late times) 

Surely the major limitation in all the analysis of LQC is that, since a satisfactory quan- 
tum gravity theory which can serve as an unambiguous starting point is not available yet, 
the theory is not developed by a systematic truncation of full quantum gravity. Another 
limitation is the restriction to isotropy and homogeneity. 

9.3 Stochastic approach 

A different approach was introduced in |321j . which starts form the observation made in 
|322j that the universe could be enlarged through an "analytic extension". In [322], such 
an extension is achieved from the geometrical construction of a semiclosed universe, namely 
a closed Friedmann model extended by gluing a given geometry to the FLRW before the 
maximum expansion. This gluing can be done in different ways, through the junction condi- 
tions. In [322j an asymptotically fiat geometry was chosen. A collection of this configuration 
(called friedmon in |321j ) was considered in }321j , in such a way that each member of the 
collection perceives the remaining systems as a perturbative effect of random character, as 
in a stochastic process. Noting that in the case of an open universe, the Friedman equation 

^^In a subsequent paper the Hamiltonian was modified to forbid the bounce at low densities [302] . 
^^An analysis along the same lines was carried out in [303j for the case e = — 1, and it was shown that the 
singularity is avoided too. 
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takes the form of the energy conservation for a harmonic oscillator, namely 

d^ + -Aa^ = 1, 
3 

a Hamiltonian can be defined by setting q = a, p = a, and the quantum theory of the 
harmonic oscillator can be developed according to [323] . A straightforward calculation leads 
to the result 

E[a\t,W)]=al, + ^V3j, 

where E is the expectation value, aci is the classical value of a, and W is the white noise 
associated to the stochastic process. One arrives at the result that the net effect of the 
environment is to preclude the collapse of the model, the minimum of the radius being large 
if A is small. 



10 Cyclic universes 

Oscillating universes have been explored in several contexts in an attempt to solve some 
problems of the standard cosmological model. The first example of such universes was that 
presented in the seminal paper by Lemaitre |3], who stated that "The solutions where the 
universe successively expands and contracts, periodically reducing to an atomic system with 
the dimensions of the solar system, have an incontestable poetic charm, and bring to mind 
the Phoenix of the legend" |3j 0. Let us briefly recall some of the issues of the standard 
model and the solution that oscillating models can provide: 

• The flatness problem. The Friedmann equation can be written as 

\^tot{t)-l\ = ^, 

As already discussed in Section [H in a situation in which the universe is dominated by 
matter or radiation, the difference \Qtot{t) ~ 1| grows as a power of t. Since present 
data indicate that Qtot is very close to 1, it must have been incredibly close to one far 
in the past, if Qtot 7^ 1 initially. This is the so-called flatness problem. As we shall 
see below, in a cyclic universe f2tot starts deviating from 1 only when a approaches its 
maximum. Since the maximum grows with the number of cycles, in a sufficiently old 
cyclic universe it may take a long time for fitot to deviate from 1 |147j . 

• The horizon problem. In the SCM, light signals can propagate only a finite distance 
between the initial singularity and a given time t, provided the energy density changes 
faster than a~^. Hence, microphysics would not have enough time to take the universe 
to its high degree of homogeneity. In the cyclic model the age of the universe is given 
by the sum of the duration of all the previous cycles. This would solve the horizon 
problem, provided correlations safely traverse the bounce. 
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Note however that Lemaitre did not produce an exphcit solution for the cychc universe. 
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Some implementations of the cyclic model may also solve the so-called "coincidence problem" 
(why did the universe begin its accelerated expansion only recently?). The model in |155] 
has its parameters tuned in such a way that the fraction of time that the universe spends in 
the coincidence state is comparable to the period of the oscillating universe. 

Oscillating models have been also used to explain the observed values of the dimensionless 
constants of nature. In [231j . the value of these constants is randomly set after a bounce 
(see also |297j ). In order to see whether cosmological evolution establishes any trend in the 
behaviour of the "constants", cyclic models were studied in |232j as solutions of varying- 
constants theories, such as the varying a theory presented in |233j . the Brans-Dicke theory, 
and the variable-speed-of-light theory [234J @ The cyclic solutions were studied both for 
non- interacting and interacting scalar field (which models evolution of the "constant") plus 
radiation, and the bounce was caused by negative-energy scalar fields. In all three theories, 
the models showed monotonic changes in the constants from cycle to cycle (the scale factor 
qualitatively behaving as explained in jl49j ). 

10.1 Thermodinamical arguments 

The existence of oscillatory solutions in the FLRW model was shown by Tolman (see [6J 
and references therein). His argument can be understood from a purely mechanic point of 
view, by modelling the Friedmann equation as a one-particle system: examination of the 
effective potential for a closed universe shows that there are oscillatory solutions for some 
values of the parameters of the model (assuming that there is a mechanism to revert the 
contraction into expansion before the singularity). These solutions are permitted from a 
thermodynamical point of view, since the matter term in the FLRW model is a perfect fluid, 
whose entropy is constant. Hence the expansion is reversible, although at a finite rate. In 
more realistic models however, entropy generation is inevitable, arising from various sources 
(such as viscosity effects from particle creation). However notice that, as discussed in [6], the 
entropy of each element of the fluid need not attain a maximum, as would be the case in an 
isolated thermodynamical system, because the energy of the fluid element is not constant. 
In fact, each time a given element of fluid returns to the same volume, its energy density is 
higher than in the previous passage through the same volume, due to a lag behind equilibrium 
conditions. The increment in the entropy leads to non-reversibility, which forbids identical 
oscillations. As a consequence of the raising energy density, the maximum value for a grows 
in each cycle @ This can be easily seen from Friedmann equation, taking the case A = 0, 
e = 1 as an example: 

a + 1 = -pa . 

3 

After one cycle, the 3-volume goes back to a value it had before when a does. Since p grows 
with the number of cycles, this growth can only be attributed to an increment in a. Hence 
a sufficiently "old cycle is strongly peaked, and Qtot remains close to 1 until a is very near 

^^A word of caution regarding this latter type of theory was issued in [235j . 

^^Notice that in these considerations neither the mechanism that allows safe passage through the singularity 
nor the details of the entropy generation are given. 
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the maximum, thus yielding a solution to the flatness problem. 

Starting from the fact that the entropy of the universe today is finite, and making the 
reasonable hypothesis that the increment in the entropy through each bounce shares this 
property, Zeldovich and Novikov [143j among others (see \287\ ) have estimated the number 
of cycles back to an initial state (which should not be singular, to keep the idea of a cyclic 
universe attractive). 

To move from qualitative arguments to actual calculations, the key issue is the production 
of entropy. The irreversible energy transfer from the gravitational field to particle generation 
was the source of entropy considered in |150j . while it was suggested in [108[ 1151] that 
black hole evaporation could be responsible for the entropy growth. An analytical study 
that showed the correctness of Tolman's arguments was presented in pH9] , where closed 
Friedmann universes with A 7^ were scrutinized, including an ad-hoc mechanism of entropy 
generation, and assuming that there is a bounce, without entering in the details of its 
realization. The entropy growth was implemented by relating the constant coming from 
the conservation laws 

Pio" = const. = Ci, 

where i denotes radiation or dust, and a = 4 or 3 respectively, to the expression for the 
entropy in each case. Let us take the case of radiation, in which 

g 

Sr = constant = -vr^/JT^a'^, 
3 

so we can set T^a^ = const. = 7. From this equation and the conservation law it follows 
that 

Or — -7— Or, 

thus linking the increment in entropy to the change in the constant appearing in the solution. 
In the same way it is shown that Cm is related to Sm through a similar expression. In |14"9] 
it was assumed that the entropy is constant within a cycle, but increases at the beginning of 
each cycle through the increment in the constants Cr and Cm- The behaviour of models with 
different combinations of matter, radiation and cosmological constant were studied for posi- 
tive and negative A. The results show that for A > A^ (where A^ = Ac{Cr, Cm)) the universe 
stops its oscillations with increasing maximum and starts an ever-expanding phase (see Fig. 
[TT]) . In other words, when the oscillations become large enough the cosmological constant 
dominates over the matter and radiation terms, the oscillations cease, and the universe enters 
a deSitter regime. If A < A^, the oscillations are not interrupted. Oscillations in anisotropic 
models were also studied in [149] . paying attention to the question of isotropization after 
a large number of oscillations. As the entropy increases, the volume of Bianchi I universes 
with A < oscillates with growing maximum amplitude, while the shear anisotropy vanishes 

70| 



''"Axisymmetric Bianchi type IX, dust Kantowski-Sachs, Bianchi IX, and some features of inhomogeneous 
cyclic cosmological models were also studied in [149) . 
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Figure 11: The plot shows that transition of a A > 0, radiation-filled universe from the 
oscillating phase to the ever-expanding phase, due to the growth of entropy (given by the 
increment in A^i = iy4Cr.„A/3). Taken from [149j . 



A more sophisticated model was studied in [236] . where FLRW two-fiuid out of equilib- 
rium models were considered. Exact solutions were found for a particular cases of the energy 
exchange, conserving the total energy. In the case of nonzero spatial curvature, cyclic models 
were shown to exist. The energy exchange between the fiuids was modelled by a function s 
such that 

p + SH'yp = s, pi + 3HTpi = -s, 

where 7 — 1 and F — 1 are the EOS parameters of each fiuid . Solutions of these equations 
along with 

=p + pi--^ 

were found in |236j for different forms of s, for the cases radiation and dust, radiation and 
scalar field, and radiation and negative vacuum energy. In the second case, a new feature 
appears (as well as the "runaway stage" mentioned in |149] ): the increment in magnitude 
of the minima in the scale factor as time increases. This was interpreted by the authors as 
a consequence of the energy exchange: the scalar field reached negative energy values after 
transferring energy to radiation. Surely this behaviour depends on the specific form of the 
function s. The examples studied in |236] suggest that caution is needed when it is said that 
cyclic models can solve the fiatness problem, since in some of them the cycles cannot become 
indefinitely large and long-lived, while in others the minimum of the expansion increases. 



87 



10.2 Realizations of the cyclic universe 



We present in this section some concrete examples of theories that yield cyclic regular so- 
lutions {i.e. which actually bounce at the minimum of the expansion without presenting 
singularities), along with some of its successes and conundrums. 



10.2.1 Changes in the matter side of EE 

One way to generate a cyclic universe is to add matter that will certainly produce a bounce, 
and consider next what conditions are to be imposed on it to produce oscillations. A nec- 
essary condition that the extrema of the expansion factor must satisfy is given hj H = 0, 
with 



in 

This amounts to p — f{p) = 0, where the function /(p) is positive. A cyclic universe has been 
generated along this line in [128] . where "wall-like" and "string-like" matter (whose energy 
scales as and respectively) generate the required /(p)lj. These rather exotic sources 
can be also thought as originating from scalar fields under the influence of a potential, using 
the procedure presented in |129j . A modification of the Friedman equation coming from 
brane models was used to fix the form of f{p) in jl37] . where 

see Sectjni The dominant component in this model is the so-called "phantom" matter, which 
has an energy-conditions-violating equation of state characterized by 

uq = — < -1. 

PQ 

Since the energy density of matter with state parameter uj scales with the expansion as 

we see that p grows with the expansion. Surely before reaching an infinite energy density, 
quantum gravity effects will take over the evolution. The somewhat paradoxical situation 
arises in which very high-density effects must be incorporated in the description of the 
universe for both very small and very large values of the scale factor. The central idea in 
[137] is that the same physics causes then the bounce and the turnaround, both governed by 
Eqn. fl257p . After a bounce, the universe follows the standard evolution until the phantom 
energy dominates. This energy may erase every trace of structure [138] . and dominates the 
evolution until high-density effects are again important, producing the turnaround. As will 
be discussed in Sect J10.3l one of the problems to be faced in the collapsing phase is the 

'^"'^ Earlier attempts along these lines, imposing that p oc —a"", and p = p oc —a^^ are respectively given 
in |377j and |378j . For a somewhat different approach, see [379]. 
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merging of black holes into a "monster black hole. The energy density the universe must 
reach in order that black holes are torn apart was shown in [137] to be 



Pbr oc Mp 



Mr 



M J 327r \1 + 3uq\' 



This energy density must be reached before the turnaround, characterized by pta = 2|o"|. 
The value a ~ rriGUT is enough for all but Planck-mass black holes to be torn apart (some 
of them evaporate before the universe enters the phantom energy stage). These Planck- 
mass remnants may help in explaining the dark matter puzzle [137j Some problems still 
remain in this model. First, the generation of structures in the contracting phase needs to 
be addressed, to see that the black hole problem does not recur. Second, as stated before, 
entropy production would lead actually to quasi-cyclic evolution. 
A similar model has been studied in |157j . given by 

H' = ^p + up' + ^, (258) 

where z/ is a real constant. Analytical solutions of this equation have been found in the case 
of dust, and their generic feature seems to be the replacement of the initial singularity by a 
bounce, some solutions displaying also a cyclic behaviour (those for A < and u < 0). 

An interesting twist to the entropy problem in cyclic universes was introduced in [279], 
where a model described by Eqn. 02581) was studied, with the cosmological constant replaced 
by a dark energy component with EOS p = up and u < —1, matter and radiation as normal 
components, and u < 0. The model takes advantage of the Big Rip phenomenon, where 
bound systems become unbound and their constituents causally disconnected as a result 
of the increasing value of the dark energy density. As a consequence of the Big Rip, the 
universe would disintegrate in a huge number of disconnected patches. The new ingredient 
of the model is that the turnaround is placed an instant before the "total Big Rip" , when 
each patch would contain almost no matter at all, and only a small amount of radiation 
|280] and dark energy. Due to the Big Rip, the huge entropy of the universe is distributed 
between the enormous number of patches, hence leaving each patch with very low entropy. 
The subsequent contraction of each patch is free of "formation of structure" problems, and 
proceeds until a bounce occurs. After the bounce, a normal inflationary phase follows (vastly 
increasing the entropy), and the cycle starts again. 



10.2.2 Cyclic universes in nonlinear electrodynamics 

As discussed in Sect J4.5[ nonlinear electrodynamics can describe a nonsingular universe. Here 
it will be shown how a cyclic model arises from the theory given by the Lagrangian |298] 

C = -^F + aF'-^ (259) 

''^Details about the evolution of this model and its relation with the so-called coincidence problem can be 
found in [155] . 
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where a and 7 are constants, with the dependence of the magnetic field on the scale factor 
given by Jif = M^ja^ (see Eqn. (12081) . The time-evolution of the scale factor can be qualita- 



tively described by the effective potential, which arises from Friedmann equation written as 
a "one-particle" system. For the case at hand, the effective potential is given by 



A B 



Ca\ 



(260) 



The constants in V{a) are given by 



A 



B 



6 



^2 



C 



7 



2^4' 



and are all positive. The analysis of V{a) and its derivatives implies solving polynomial 
equations in a, which can be reduced to cubic equations through the substitution z = a^. 
The existence and features of the roots of such equations are discussed in |309] . A key point 
to the analysis is the sign of defined as follows. For a general cubic equation 



X + px = q, 



the discriminant D is given by 



D 



I)' 



We will denote by Dy the discriminant corresponding to the potential and Dy' that of the 
derivative of V. From the behaviour of the potential and its derivatives for a ^ and a 00 
we see that only one or three zeros of the potential are allowed. The case of interest here 
(given by Dy > 0, Dy = 0) is plotted in Fig. [12], which shows the qualitative behavior of 
the potential for typical values of the parameters. The model is nonsingular for any value of 




Figure 12: Qualitative plot of the effective potential for Dy > 0, Dy 
line corresponds to e = 1. 



0. The lower dotted 
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e, and a cyclic model is obtained for e = 1. 

This setting was generalized in |419j . where the Lagrangian 

= + (261) 

was considered, with a,/3 and fi constants. As shown in [419J, four distinct phases can 
be described with this Lagrangian: a bounce, a radiation era, an acceleration era and a 
turnaround. This unity of four stages, christened tetraktys in [419] . constitutes an eternal 
cyclic configuration. The cyclic behavior is a manifestation of the invariance under the dual 
map of the scale factor a{t) —>■ l/a{t), a consequence of the corresponding inverse symmetry 
of the Lagrangian (126 ip wrt the electromagnetic field [F — > 1/F, where F = F^'^F^j^). 
Restricting to a magnetic universe, as defined in Sect J4.5.Tl the Lagrangian Ct yields for the 
energy density and pressure given in equations fll86til87p : 

P = -c^F^ + \f+^-§,. (262) 

p = F^ + —F- — — + — -. 263 

^ 312 3F 3F2 ^ ' 

As we saw in Sect J4.5.3l for any Lagrangian that is a polynomial in F, 

= Jiff) a~^. 

As discussed in [419] . the combined system of equations of the FLRW metric and the mag- 
netic field described by General Relativity and NLED, are such that the negative energy 
density contributions coming from Ci and £4 never overcome the positive terms arising from 
C2 and £3. Before reaching undesirable negative energy density values, the universe bounces 
(for very large values of the field) and bounces back (in the other extreme, that is, for very 
small values) to precisely avoid this difficulty. These events occur at the values Pb = Pta = 0, 
which follow from Friedmann's equation in the case e = 0. Notice that this is not an extra 
condition imposed by hand but a direct consequence of the dynamics described by Ct- 

Let us now turn to the generic conditions needed for the universe to have a bounce and a 
phase of accelerated expansion. From Einstein's equations, the acceleration of the universe 
is related to its matter content by 

3- = -l(p + 3p). (264) 
a 2 

In order to have an accelerated universe, matter must satisfy the constraint (p + 3p) < 0, 
which translates into 

U > j^. (265) 

It follows that any nonlinear electromagnetic theory that satisfies this inequality yields ac- 
celerated expansion. In the present model, the terms C2 and £4 produce negative accel- 
eration and Ci and £3 yield infiationary regimes (a > 0). Raychaudhuri's equation shows 
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that : imposes further restrictions on a{t) at a bounce. Indeed, the existence of a mini- 
mum (or a maximum) for the scale factor imphes that at the bounce point B the inequahty 
{Pb + Sps) < (or, respectively, {ps + ^Pb) > 0) must be satisfied. Note that, as already 
mentioned, at any extremum (maximum or minimum) of the scale factor the energy density 
is zero. Four distinct periods can be identified according to the dominance of each term of 
the energy density. The early regime (driven by the F"^ term); the radiation era (where the 
equation of state p = l/3p controls the expansion); the third accelerated evolution (where 
the 1/F term is the most important one) and finally the last era where the dominates 
and in which the expansion stops, the universe bounces back and starts to collapse. The 
bounce (for an Euclidean section) was discussed in Sect J4. 5. 31 The standard, Maxwellian 
term dominates in the intermediate regime. Due to the dependence on of the field, this 
phase is defined by J^"^ >> yielding the approximation 



P 



2 

p f» (266) 

When the universe becomes larger, negative powers of F dominate and the energy density 
becomes typical of an accelerated universe, that is: 



P 



1 /i^ 



2 ^2 



In the regime between the radiation and the acceleration eras, the energy content is described 
by 

P=^ + 



2 



or, in terms of the scale factor. 



1 ^ /i ^ 



2 2^2 V ; 



For small a it is the ordinary radiation term that dominates. The 1/F term takes over 
only after a = \/ M^j /i, and grows without bound afterwards. Using this matter density in 
Eqn. fl264p gives 

a < 1 _3/ 

a 2 2e^2 

To get a regime of accelerated expansion, we must have 
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which imphes that the universe will accelerate for a > a^, with 

^0 



For very large values of the scale factor, the energy density can be approximated by 

(269) 

and the model goes from an accelerated regime to a phase in which the acceleration is 
negative. When the field attains the value Fta = IQa'^fi'^ the universe stops expanding and 
turns to a collapsing phase. The scale factor attains its maximum value 



a. 



4 _ '^^0 



jr. 



'2 



Analytic forms for the scale factor in each regime can be found in |419j . 



10.2.3 Cyclic universes in loop quantum gravity 

There are realizations of cyclic models in the effective equations for loop quantum gravity 
(some features of which have been presented in Sect j9.2p . As discussed in Sect J9.2[ the 
Klein-Gordon equation for a scalar field under the influence of a potential, the Friedmann 
and Raychaudhuri's equations in the semiclassical regime are modified due to quantum 
gravity effects (see Eqns.( l25Hl256l) ). It was shown in |134] that positively curved universes 
sourced by a massless scalar field can undergo repeated expansions and contractions due 
to the modifications described above. This was achieved by rewriting Eqns. fl254ti256|) in 
the form of the classical FLRW model with the addition of matter described by an effective 
equation of state, given by 

^ ^ = _1 + (i _ ^ d\nD \ 

Peff 02 + 2DV V 6 dlna J ' 

A violation of the null energy condition, leading to a bounce, is accomplished when < — 1, 
which amounts to rflnD/rflna > 6, or a < 0.914a^, [134j . with 

D{q) = (^^J g^/2{7[(g + 1)"/^ - \q - - llg[(g + 1)7/4} - |g - ir/%n(g - 1)]^, 

with q = a? /a1 and a1 = 7£p;j'/3, where 7 ~ 0.13, and j is a quantization parameter, which 
takes half-integer values. When V = 0, u is independent of the kinetic energy of the field, 
and an oscillatory behaviour follows. The addition of a potential leads to the interruption of 
the cycles as soon as the potential dominates the motion (in analogy to what was discussed 
in Sect ilO.l1 for the cosmological constant), and a period of inflation may follow [134i| . This 
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analysis was later extended to the case of spatially flat universes, with both negative an 
positive potentials |135[ 1332] . 

Yet another realization of a cyclic universe in this scenario is the so-called emergent 
universe from a loop |136] . As mentioned in Sect l8.3l the Einstein universe is unstable, so 
perturbations drive the universe away from this state. This situation partially changes when 
loop quantum gravity corrections are considered. Using a phase-space analysis, it was shown 
in |136j that a new static solution appears in the semiclassical regime {a < a^) for positive 
potentials (for < this is the only solution). This new solution (called loop static, LS) is 
stable, and the universe oscillates around it, for V < K, with = 39/(1367r/p;a^), while for 
V > K, the equilibrium point corresponding to LS merges with that of the Einstein universe. 
So in the model proposed in |136j . the universe is initially at, or in the neighbourhood of the 
static point LS, with in the plateau region of the potential with (pi > 0. After undergoing 
a series of non-singular oscillations in a (possibly) past-eternal phase, while the field evolves 
monotonically along the potential, the cycles are eventually broken as the magnitude of the 
potential increases, and the universe enters an inflationary epoch. For this model to work, 
the potential must be such that dV/dcj) for cf) —oo, and increase monotonically to 
exit the cycles An example of a suitable potential is given by 



V = a 



exp 



m/Vs) - 1 



2 



where a and (3 are parameters that may be constrained by the CMB spectrum. As in the 
case of the classical emergent universe discussed in Sect J8.3t there are some fine-tuning issues: 
the scalar field must start in the asymptotically low-energy region of V. 



10.2.4 The cyclic universe based on the ekpyrotic universe 

The starting point of the ekpyrotic scenario ^105j is five-dimensional heterotic M-theory [409], 
where the fifth dimension terminates at two boundary Z2 branes, one of which is identified 
with the visible universe. There are two different versions of the ekpyrotic scenario, the old 
[116] , where there is a bulk brane between the boundary branes and the new 1117], where 
only the boundary branes are present |410j . The initial state in both cases is supposed to 
be the vacuum state, where the branes are flat, parallel and empty. The branes are drawn 
together by the action of an attractive potential, and collide inelastically over cosmological 
times. Part of the kinetic energy is transferred to the branes and used to create matter and 
radiation. After the collision, the universe enters a "standard" big bang phase, until dark 
energy domination at the end of the matter era, which causes an accelerated expansion, 
diluting the content of the universe. The whole process can be described by a 4-d effective 
theory, with the action (in the Einstein frame) given by 

Se = J d'x Qi? - i(Vvp)^ - Vi^) 

Other constraints are imposed by succesful reheating. 
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plus higher-order corrections, where the conveniently-tailored potential V{(f) is responsible 
for the main features of the model. The potential is slightly positive for (f > 0, and goes to 
zero as </9 — i> — oo. For (f < 0, the potential has a minimum and is very steep and negative. 
The minimum corresponds to the close approach of the branes, which happens at such short 
distances that quantum gravity effects are relevant. The field if moves rapidly through the 
minimum, and the branes collide as ip —oo. Both the old and the new model were shown 
to have problems due to excessive fine-tuning [HI], so a cyclic version was introduced |119j . 

In the cyclic ekpyrotic model, it is assumed that the interbrane potential is the same 
before and after collision (instead of being zero, as in the non-cyclic model). After the branes 
bounce and fly apart, the interbrane potential ultimately causes them to draw together and 
collide again. To ensure cyclic behavior the potential must vary from negative to positive 
values |119j . The model may be adjusted in such a way that, at distances corresponding to 
the present-day separation between the branes, the inter-brane potential energy density is 
positive and corresponds to the currently observed dark energy, providing roughly 70% of the 
critical density today. The cosmic acceleration restores the Universe to a nearly vacuous state 
and as the brane separation decreases, the interbrane potential becomes negative. As the 
branes approach one another, the scale factor of the Universe, in the conventional Einstein 
description, changes from expansion to contraction. When the branes collide and bounce, 
matter and radiation are produced and there is a second reversal transforming contraction 
to expansion so a new cycle can begin [119j . Figure [TO. 2. 41 shows a plot of several forms of 
the potential that allow for a cyclic universe in this scenario |121j . A qualitative description 
of the model can be given in terms of this figure as follows. Currently, the field is in region 
(a), at the point indicated with a dark circle, where the potential is flat and drives cosmic 
acceleration. Eventually, the field rolls towards negative values of V (region b), where cosmic 
expansion stops and the universe (being nearly vacuous as a consequence of the acceleration 
phase) enters a phase of slow contraction, where the spectrum of density perturbations is 
generated from quantum fluctuations in ip. In region (c) the kinetic energy of if dominates 
the energy density. At the bounce, part of this kinetic energy is converted into matter 
and radiation, while the perturbations in ip are imprinted as density fluctuations in the 
matter/radiation fluid. Meanwhile the field quickly returns back to (a) where it comes to 
a stop, and the universe enters the radiation-dominated era, so commencing the next cycle. 
As recognized by its authors, the model presents two weak points (as is the case with many 
cyclic models): the passage through the would-be singular point, and the propagation of 
perturbations 0. It is difficult to achieve the bounce without passing from the semi-classical 
regime to the high-energy fully quantum regime, where our use of the effective 4-dimensional 
theory breaks down. The problem is that the kinetic energy and the Hubble rate typically 
reach Planckian scale as the branes approach. In fact, in the semi-classical regime where 
loop corrections can be applied, brane collision may be prevented 0. 

Recently, a "new ekpyrotic cosmology" was presented in [411] . where a NEC- violating 
ghost condensate was merged with an ekpyrotic phase to generate a non-singular bouncing 

''''This second problem will be discussed in SeclTTl 

''^Some other problems of the model were discussed by Linde [122] . 
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Figure 13: Plot of several allowed forms for V{(p) . 



cosmology. The authors claim to obtain a pre-bounce scale-invariant spectrum using the 
mechanism of entropy perturbation generation |412j . This is accomplished by having two 
ekpyrotic scalar fields rolling down their respective negative exponential potentials, and 
having its own higher-derivative kinetic function. Notice however that the results of this 
model have been challenged in |413] . 

10.2.5 Oscillatory universe from the Steady State model 

The Steady State model [264] was proposed as alternative to the Big Bang model, and has 
fallen into disfavor because the observations of the CMB. However, its authors have advanced 
a new scenario, called the quasi-steady state model (QSSC, see |144[ I243t 12441 [2^ ). In this 
model, the singularity is avoided by the action of a scalar field C{x), which creates matter in 
compliance with the Weyl postulate and the cosmological principle, and has negative energy 
and stresses. The cyclic solutions in the QSSC can be expected from physical grounds 
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as follows |243] . To create a particle, C{x) must have energy- momentum equal or larger 
than that of the particle. When C is above the threshold, it creates particles and fuels the 
spacetime expansion (since it has negative stresses). To this overall expansion an oscillation 
is superimposed. The creation of particles and the expansion set C below the threshold, 
slowing down the number of created particles, and the expansion. Here, the cosmological 
constant takes control and causes contraction. The contraction rises the background level 
of the C field, and the cycle starts again. As shown in |243j . a solution to the EOM of this 
theory in the FLRW setting that oscillates in this way is given by 

a{t) =e'/^ {l + r]cose{t)) , 

with 6{t) fti 27rt/Q, where P is the long term "steady state" time scale of expansion, Q is 
the period of a single oscillation (with P » Q), and rj is a parameter . 

10.2.6 Other models 

Due to the recently discovered dark energy component of the universe, several forms for the 
dependence of the EOS parameter with the redshift have been analyzed [420] . In fact, some 
data suggest that uj{z) evolved from a value larger that —1 to a value smaller that —1 at 
some recent redshift. One of the models that describes this crossing is the quintom model 
|139j , where to is parameterized as 

a;(lna) = c^o + t^^i cos[y41n(a/ac)], (270) 

with luq, uji, A, and ac to be fitted by observations 0. It was shown in [140j that for a certain 
choice of the parameters, a universe filled with quintom matter (that is, matter wit u given 
by Eqn. fl270p plus radiation and normal matter expands and contracts cyclically, yielding an 
inflationary period at the beginning of each cycle, and an acceleration period at the end. 

Perhaps it is convenient at this point to remind that a closed universe has not been 
discarded by observation yet (and in fact, cannot be discarded with certainty due to the 
errors inherent to any experiment), though theoretical prejudice and observation tend to 
favor Q = 1. As we saw in Chapter [5l a nonzero bulk viscosity ( modifies the fluid pressure 
according to 

P = Po- KH, 

where po is the equilibrium pressure. The asymmetry in the pressure depending on the 
sign of H causes the increment in energy and entropy, leading to ever-increasing cycles. It 
was shown in |146] that a similar asymmetry can be caused by scalar fields in a pure non- 
dissipative setting. Starting from a FLRW setting plus a scalar field under the influence 
of a potential which displays a minimum, an asymmetry in the pressure, given by p ~ —p 
for H > 0, and p ^ p for H < is generated by the oscillations of the field around the 
minimum |146j . By imposing appropriate conditions to force a bounce [a —>■ a, a —>■ —a, 

''^Constraints on this form of dark energy were studied in [142j . 
''''Similar ideas were studied in |141| . 
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^ 0, — s> 0), it was shown that there is an in increment in the maximum radius of 
expansion of the universe in each cycle, due to conversion of work, done during expansion, 
into expansion energy. The flatness problem is gradually ameliorated in this model, since 
the universe becomes considerably long-lived and more fiat after each expansion. 

To close this section, other models of a cyclic evolution for the universe are listed next: 

• String theory-inspired cyclic universes inspired, starting from the property that there 
exists a minimal length, ^pl. See |130j . 

• Classical spinor field under the influence of a quartic potential in a FLRW background 
was discussed in |154j . It was shown that V = \^ + mipip — X{ijjipY gives rise to 
oscillations in the scale factor, for certain choices of the parameters. 

• A cyclic scenario that takes into account matter and radiation evolution if the proton 
has a finite lifetime was studied in |287j . 



10.3 Issues of the cyclic models 

Cyclic universes are not free of problems. As was put forward in |123j . during a matter- 
dominated cycle, black holes with masses ranging from stellar to galactic will form. During 
the contracting phase they will coalesce into a "monster black hole" with mass equal to the 
mass of the universe. Its entropy can be estimated by 

S = ^A = 27iR^ = SvrM^ > 10^^\ 

where the mass within one Hubble volume (~ IO'^^Mq) was used. However, the entropy 
of the radiation in the present Hubble volume is ~ 10^^, in such a way that black hole 
formation in a previous cycle would lead to a huge excess of entropy generation. In this 
scenario, the excess must have somehow been eliminated by the bounce. But there are some 
ways out of this problem. Sikkema and Israel |124j have suggested that the inner horizon 
of the monster Kerr black hole absorbs strongly blue-shifted gravitational radiation emitted 
during the last moments of the collapse. This radiation increases the mass of the core of the 
black hole by a huge amount, rapidly reaching Planckian values, and correspondingly greatly 
reduces its specific entropy. If quantum effects produce a bounce, this process would allow 
the expansion to begin from a state of relatively low disorder Durrer and Laukenmann |147j 
have proposed another solution to Penrose's problem. They have remarked that black hole 
thermodynamics is valid only in asymptotically fiat exteriors (a fact which was also noted in 
[124] ). They also noted that the entropy in the radiation we observe today is actually due 
to the previous matter cycle, which may have had shorter duration than the current cycle, 
leading to less clumping and consequently less entropy production 

Another issue of cyclic models was raised in |156] , where the evolution of a cosmic string 
network was considered in a bouncing universe. It was shown that the string network displays 

''^Gravitational perturbations were also studied in [147j 
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an asymmetric behaviour between the contraction and expansion epochs. In particular, while 
during expansion a cosmic string network will quickly evolve towards a linear scaling regime, 
in a phase of collapse it would asymptotically behave like a radiation fluid. A cosmic string 
network will add a significant contribution, in the form of radiation, to the energy (and 
hence also entropy) budget of a contracting universe, which will become ever more important 
as the contraction proceeds. Hence it establishes the need for a suitable entropy dilution 
mechanism. This process will also operate, mutatis mutandis, for other stable topological 
defects. Conversely, if direct evidence is found for the presence of topological defects (with 
a given energy scale) in the early universe, their existence alone will impose constraints on 
the existence and characteristics of any previous phases 

11 Perturbations in bouncing universes 

As discussed in the Introduction, inflation can solve many of the shortcomings of the SCM, 
but it also has problems of its own. Bouncing models may provide an alternative (or perhaps 
a complement) to standard inflation, since in principle the problems of the SCM come from a 
"shortage of time" for things to happen early after the big bang [HI]. The arguments in SectH] 
show that an accelerated contraction has the necessary features to solve the problems of the 
SCM |341j . Let us recall that if in the contracting phase the Hubble radius decreases faster 
than the physical wavelength corresponding to fixed comoving scales, quantum fluctuations 
on microscopic scales can be stretched to scales which are cosmological at the present time, 
exactly as it happens in inflationary models (see for instance [91]). Figure [TT] shows a sketch 
of the structure of a space-time in which standard inflation starts at ti and ends at tn. During 
inflation, the Hubble radius H'^(t) is constant, and it grows linearly afterwards, while the 
physical length corresponding to a fixed co-moving scale increases exponentially during the 
period of inflation, and then grows less fast than H~^{t). The figure shows that for a given 
k, the fluctuation can be (causally) produced well inside the Hubble radius, "leave " H~^(t), 
and "re-enter" in an appropriate way to describe the structures we observe today. 

Figure [TT] shows a universe that undergoes a contracting phase, a bounce, and then enters 
an expanding epoch, assumed to be that of the SCM. In this case, the Hubble radius decreases 
relative to a fixed comoving scale during the contracting phase, and increases faster in the 
expanding phase. Fluctuations of cosmological interest today are generated sub- Hubble but 
propagate outside the Hubble radius for a long time interval. There is however, one main 
difference with respect to the standard inflationary scenario. In the latter the curvature scale 
R oc is (almost) constant, while in the former, it grows until it reaches a maximum and 
then decreases 0. This difference may lead to observational consequences 0, particularly 
regarding the generation of a primordial spectrum of inhomogeneities through parametric 
amplification of the quantum fluctuations of the background fields in their vacuum state [174j . 

''^This assertion is valid in models in which quantum effects intervene in such a way that Rmax oc A^^„, 
which is the case of loop quantum gravity for instance, where Xmin oc ipi. For the models in which H reaches 
a null value, H"^ can be replaced by ic = \/cL^/a", see Eqn. (|355l 

*°See [342j for a qualitative discussion of these consequences in the case of string pre-big-bang cosmology. 
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Figure 14: Behavior of the co moving scale k and of the Hubble radius if ^ as a function of 
time in inflation. Taken from [393]. 



These, when decomposed in Fourier modes, satisfy a canonical Schrodinger-hke equation, 
whose effective potential is determined by the so-called "pump fleld" , which depends in its 
turn on the background geometry. There are then two properties of the background in a 
bouncing universe that can affect the flnal form of the perturbation spectra |342j : (1) the 
growth of the curvature scale, and (2) the flelds which, together with the gravitational fleld, 
determine the background. Property (1) has two important consequences. The flrst, is that 
bouncing scenarios may lead to "blue" {i.e. growing with frequency) metric perturbation 
spectra, instead of being flat, or decreasing ("red"), as in standard inflation. A growing 
spectrum leads to the formation of relic backgrounds whose amplitude is higher at higher 
frequency, hence more easily detectable. A typical example is that of gravitational waves 
in SPPB |342j (see Eqn. ( 12791) ). The second is that the growth of the curvature may also 
force the comoving amplitude of perturbations to grow (instead of being frozen) outside the 
horizon (see [385J for this effect in the SPBB) 0. 

Regarding Property 2, one of the interesting consequences is the ampliflcation of the 
fluctuations of the EM fleld, due for instance to the non-minimal coupling with a scalar fleld 
(such as the dilaton, or the scalar fleld in WIST, see Sect jll.4|) . A relic background of scalar 
particles is also generated, which may be related to dark matter |343j . 

There is yet another salient feature of the perturbations in a bouncing universe. Since in 
the far past of this type of models the universe is assumed to be almost flat, one can impose 
vacuum initial conditions for the perturbations based on simple quantum fleld theory in flat 

Consequently, special attention must be taken in the application of linear perturbation theory, see [385j . 
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Figure 15: Behavior of the co moving scale k and of the Hubble radius if ^ as a function of 
time in a bouncing universe. Taken from [393]. 



space [366], instead of having to set initial conditions in a high-curvature regime. 

It must be remarked that solving for the perturbations in bouncing models is in principle 
a nontrivial task, since there are potential ambiguities that may arise at the bounce, not 
present in standard inflation!^. Two views can be taken to tackle the study of perturbations 
in such a scenario. The first one is to devise first a detailed model of the bounce, and then 
study the properties of the post-bounce perturbations. The problem in this case is that total 
control of the high-energy physics involved in the bounce is needed, which is not always 
achieved. It may also happen that the bouncing solution under scrutiny is quite artificial 
from the physical point of view, as for instance if it is not embedded in any fundamental 
theory. But in any case some lessons may be extracted from the examples, as we shall see 
in Sect irm 

A second attitude is to make some simplifying assumptions and try to work out predic- 
tions that are independent of the UV physics that most surely governs the bounce. This 
possibility has led to a great debate [345] . In particular, in order to avoid the specification of 
the details near the high-curvature regime, matching conditions are used, leading to ambi- 
guities. The dependence of the post-bounce spectrum on the matching conditions has been 
addressed by many authors, as will be discussed in Sect J11.3l 

At this point, it is perhaps necessary to say that there are at least two alternative pro- 

^^For instance, at the bounce the comoving Hubble scale diverges. Hence all scales are inside the Hubble 
scale, at least for an instant. However, there are some issues common to both scenarios, such as the 
transplanckian problem (see for instance |400j ) . 
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cedures to deal with gravitational perturbations in a relativistic setting. Since Lifshitz's 
original paper |166] , it has been a common practice to start the examination of the theory of 
perturbations of General Relativity by considering variations of non-observable quantities, 
such as 5g^y. The main drawback of this procedure is that it mixes true perturbations and ar- 
bitrary (infinitesimal) coordinate transformations, which are unphysical. As shown in [349j . 
|192j . |193] . |174j . this problem can be solved by adopting gauge-independent combinations 
of the perturbed quantities expressed in terms of the metric tensor and its derivatives. The 
dynamics of these gauge-independent variables is then provided by the EE. 

A second method exists, based on the quasi-Maxwellian (QM) formulation of Einstein's 
equations. The advantage of this method is that it is gauge- independent from the start, thus 
deahng only with observable quantities [173], [M], [ITOl [HH [172] [176] . We shall briefly 
review both methods in Sect lll.51 including a summary of the relation between them. 

In the next sections we shall discuss examples of the two approaches. From an ob- 
servational point of view, the crucial question is whether bouncing models can furnish a 
nearly-scale invariant spectrum of adiabatic scalar perturbations after the bounce, as de- 
manded by the measurements of the WMAP |394] . Sloan survey [395] . and 2df |396j . It is 
also of interest to see if bouncing solutions lead to observable consequences that are markedly 
different from those of inflation (see Sect. [T^ . 

11.1 Regular models 

In the previous chapters, we have seen that it is possible to generate bouncing models 
in a wide choice of scenarios, essentially by any of the mechanisms presented in Sect Jl.ll 
Obviously, the outcome is very dependent on the choice, but specific models can be sometimes 
useful in the hope of extracting tendencies of a more general behaviour. In this sense, scalar, 
vector, and tensor perturbations have been studied in many exact backgrounds displaying a 
bounce. An incomplete list includes the following: 

• General relativity with radiation and a free scalar field having negative energy [187j . 

• Two scalar fields [TO (19511555] . 

• A 5d Randall-Sundrum model with radiation, in which the extra dimension is timelike 

• Two perfect fluids |179] . 

• A nonlinear EM Lagrangian [34 7j . 

• A scalar field with higher-order corrections from string theory, with an exponential 
potential (this case covers the SPBB and the first version of the Ekpyrotic universe) 
[316] . 

• A non-canonical scalar field, with Lagrangian C = p{X,(f)), where X = 1 / 2g^^ d ^(pd^ip 
[360] . 
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• Bounce due to quantum cosmological effects using Bohmian solutions of the canonical 
Wheeler-de Witt equation [358] . 

• Non-local dilaton potential stemming from string theory [351] . 

We shall present next a short discussion of scalar, tensor, and vector perturbations in 
some of these scenarios. 

11.1.1 Scalar perturbations 

The evolution of scalar perturbations through a bounce has been a subject of intense debate 
(see references in |35Uj ). A consensus for the case of a two-component bouncing model in GR 
seems to have been reached after |350j . This model is described by a flat FLRW metric, and 
one of the components has negative energy density (to produce the bounce) and is important 
only near the bounce. The components interact only gravitationally, and the component that 
dominates away from the bounce has an intrinsic isocurvature mode, in order to describe 
scalar fields or perfect fluids. The result obtained in |350] is that the spectrum of the growing 
mode of the Bardeen potential in the pre-bounce is transferred to a decaying mode in the 
post-bounce @, @- 

Since the phenomenology associated to the decaying mode is known to differ from ob- 
servation |352j . we may ask what can be done to lift the negative result of |350] . One 
possibility is to allow the fluids to interact. Another one is to incorporate in the background 
solution the decay of the normal component to radiation 0. Yet another possibility is to 
consider higher-order corrections. This has been done in several string-inspired models 0, 
in the gravi-dilaton regime by exploring regular backgrounds (such as those presented in 
Sect J3.3.11) . as in [354J . [355] . [3^6J . [353J . The results presented in these articles show that 
although it may be possible to generate a nearly scale-invariant spectrum in the pre-bounce 
phase, it corresponds to the decaying mode in the expanding phase 0. An exception is the 
model presented in Sect lll.LTl Another exception may be the ekpyrotic model, where there 
are results indicating that a scale-invariant spectrum may be obtained in the post-bounce 
phase [356] @ 

Another set of models comes from the quantum evolution of the universe. As discussed 
in Sect J9.1l bouncing solutions are possible (without the need of a "phantom" field) in the 
context of the WdW equation, when the Bohm-de Broglie interpretation is used in the mini- 
superspace approach. A feature of this scenario is that a full quantum treatment of both 
background and perturbations is possible [357] I359j . The model analyzed in [358] is GR plus 

^■^These result is supported by the references cited in |350j and also by the results in [179j . 
^''Notice that mode- mixing is possible with e = 1, as for instance in |344| . 
^^See Sect 13 .2.91 and 327] for an exact solution that has this feature. 

^^The string pre-big-bang model without corrections furnishes a highly blue-tilted spectrum = 4 of 
scalar perturbations [385]. 

^^The SPBB model may yield the right spectrum when axion fluctuations are considered [361] . 

*®See ref. [362] for another model in which the growing mode in the contracting phase goes over into the 
dominant mode in the post-bounce phase. 
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a perfect fluid, in which the scalar perturbations can be described in terms of a single degree 
of freedom, related to the Bardeen potential $ (see Appendix). The Bohmian quantum 
trajectory for the scale factor is given by 



a(T) = oo 



3(1-1^) 



(271) 



with p = up. The normal modes of the scalar perturbation satisfy the equation 



uk^ - — ] Vk = 0, 



(272) 



where a prime means derivative wrt conformal time. Following the usual procedure of ex- 
panding the modes for large (negative and positive) values of T, matching the expansions, 
and then transforming to the Bardeen potential, the power spectrum defined by 



yields for the post-bounce phase |358j 



9P 

71^ 



12uj 



1 + 300 

An analogous calculation for the tensor modes gives 



l + 3uj' 



(273) 



(274) 



(275) 



Notice that a scale-invariant spectrum follows both for the scalar and the tensor perturbations 
for the case of dust {u = 0), which is the fluid supposed to dominate the evolution at the 
time of the matching of the solutions (not necessarily the same governing at the time of the 
bounce) |358j . An important lesson that follows from this example and the one presented 
in [179] (see Sect Jll.Ol) is that the spectral index is quite insensitive to the details of the 
bounce, being determined mostly by the dominant component. The example also shows that 
the bounce is important in the mixing of the modes, which is relevant for the amplitude of 
the modes in the post-bounce phase. 



11.1.2 Vector perturbations in a contracting background 

It is a well-known result of perturbation theory that vector perturbations (VPs) only exhibit 
decreasing solutions in the context of an expanding Universe (see for instance |174j ) 0. 
However, as shown in |189j . VPs can increase in a contracting fiat background, with a 



®^ Another interesting result is that the simplest models of inflation do not produce VPs, see for instance 
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perfect fluid as source. Hence, they might provide a signature of a bounce. As shown in the 
Appendix, the relevant equations are S\ = C\/a?^ where C is a constant, and 



n « (276) 

Note that increases for = 0, and stays constant for radiation, but S^y. always increases 
for decreasing a. As argued in |189] . VPs cannot be neglected in the SPBB scenario, in such 
a way that near a bounce, the metric perturbations may become too large for the use of 
linear theory (depending on the value of the C\) @. Related results were presented in [365j, 
where it was also shown that the growing vector mode matches with a decaying mode after 
the curvature bounce, in the context of a low-energy flat gravi-dilaton model lj. 

Since many bouncing models are generated by a scalar fleld, a relevant question is whether 
VP are important in this type of scenarios. One important point is that VPs are not 
supported by a scalar fleld at flrst order. At second order, the scalar, vector, and tensor 
modes couple, and VPs can be generated by scalar-scalar mode couplings [364] . Considering 
exponential potentials and power-law solutions, the ratio of the amplitudes of second order 
vector perturbations in contracting and expanding phases was studied in |364] . The relative 
magnitudes of the second order vector perturbations in the two phases depend on the scaling 
solutions chosen, but at least in one of the examples studied (dust-like collapse, [198] ). the 
observable differences between the collapsing models and the inflationary scenario could be 
large, assuming that the transition between the two phases does not signiflcantly alter the 
ratio. 

11.1.3 Tensor perturbations 

The spectrum of gravitational waves can be a very powerful tool to discriminate between 
different models of the universe, since gravitational waves decouple very early from matter 
and travel undisturbed, as opposed to EM waves. In particular, in the context of the SPBB 
scenario, the ampliflcation of tensor perturbations is greatly enhanced wrt the standard 
inflationary scenario for large comoving wavenumber h |383] . This result was conflrmed in 
j , with a gravi-dilaton background solution of the EOM 

g: = I (d,^d^^ - Ui:d^^\ , (277) 



given by 



2 V"''' " ' 2 / 

□(^ = 0, (278) 



a(r/) = (-r/)^/^ (^(r^) = ^ ln(-?7) + const., 

1 + V 3 



Quantum corrections to the evolution of vector modes were studied in the context of loop quantum 
gravity in [367j 
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This is not nececssarily so in multidimensional cosmological models, also analyzed in [365j . 
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the typical amplitude for the normalized vacuum tensor fluctuations outside of the horizon 
over a scale is given by 



(279) 



PI 



where Hi ^ (ai?7i)~^ is the final contraction scale o, while the result in the standard 
inflationary expansion does not have the In dependence (see for instance |398] ). The possible 
influence of the nonperturbative phase, where the curvature and the dilaton are very large, 
was studied by imposing a bouncing solution in [38^, and by taking into account higher- 
derivative a' and quantum corrections (see Sect jS.S.l]) |387j . [384] . The results in these 
papers show that that the low frequency modes, crossing the horizon in the low-curvature 
regime, are unaffected by higher-order corrections, and also that the shape of the spectrum 
of the relic graviton background, obtained in the context of the pre-Big Bang scenario, is 
strongly model-dependent. 

This analysis was continued in |399j , where cosmological perturbations in the low-energy 
string effective action with a dilaton coupling -F(0) were studied, with the addition of a 
Gauss-Bonnet term, a kinetic term of the type (V^)"^, and a potential V{(f)). Scale-invariant 
spectra in the string frame and a suppressed tensor-to-scalar ratio were obtained by imposing 
slow-roll inflation in the Einstein frame. The results show that it is practically impossible to 
obtain these conditions without the second-order corrections given by Eq. fll54p . both with 
and without the Gauss-Bonnet term. 

Analytic and numerical results for the tensor post-bounce spectrum have been obtained 
for a two-component model defined by p± = uj±p± |179] . The flat background is given by 



a(r) = flo 1 + ^ 



r2 



with 



^0 



« = 777 7> «0 = , Tq 



3(^-^+) V7+/ ^Pi 7i 

7+ and 7_ are constants, with 7_ < 0, to produce the bounce. The tensor spectrum, 
assuming that —1/3 < uo^ < 1, and that the potential that arises from Eqn. fl355p has only 
one extremum at r = 0, is given by [179j Ph oc fc"^, where 



12w^ 



l + 3cj+ 

Note that the spectral index does not depend on the EOS parameter of the "exotic" fluid 
(contrary to the case of the spectral index for the scalar perturbations). This was to be 



®^ Scalar perturbations of this model were also investigated in [385j . and present amplitudes and spectra 
similar to the tensor perturbations. 
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expected since large wavelengths are comparable to the curvature scale of the background 
at a time when the universe is still far from the bounce, so the behaviour obtained in this 
case can be taken as generic, i.e. independent of the details of the bounce. 

Yet another example of the calculation of a tensor spectrum in a bouncing model was 
presented in Sect Jll.LTl based on the quantum evolution (using the Bohmian quantum 
trajectory) of a universe dedcribed by GR plus a perfect fluid. The result is (see the comments 
after Eqn.([27SD) 

In fact, the tensor-to-scalar ratio in this model was estimated as T/S* ~ 5.2 x 10~^, and the 
the characteristic bounce length-scale Lq ~ 1500£p;, (assuming that < 1.01) which is a 
value in the range in which quantum effects are expected to be relevant, while at the same 
time the Wheeler-de Witt equation is valid (without corrections from stringy/loop effects). 

11.2 Scalar perturbations in exact models using the quasi-Maxwellian 
framework 

As mentioned in the introduction of this chapter, perturbations can also be studied using the 
quasi-Maxwellian (quasi-Maxwellian) method. In this section we apply it to two exact bounc- 
ing solutions. The first one is generated by the non-minimal coupling of the electromagnetic 
field with gravity (see Sect. (14. 4p ). As discussed in the Appendix, in the quasi-Maxwellian 
formalism the scalar perturbations are completely described by the variables E and S, which 
obey the equations (I376p - fl378p : 

E = -^-±^pE--eE, 

2 3 ' 



E. 



6A / P\ 1 

with p = Xp, and k is the wave number (the subindex k in E and S has been omited). 
Combining these, we obtain the equation for the time evolution of the electric part of the 
perturbed Weyl tensor: 

E + e(^ + \] e + EX = 0, (281) 



^3 

where X is a function of the background functions given by 

^ ^3e + e / 2\ 2 + 3A^2 

Defining a new function g{t) hj g = E a^'^, where a = — (4+3A)/2, we obtain from Eqn. (12811) 

9 + xit)9 = 0, (282) 
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where lj 2 

X{t)^a- + (-) +X. (283) 

In the case of the bouncing universe given by Eqn. fll8ip . we have 

{t^ + al)^ g+{at^ + (3al)g = 0, (284) 

where a = - 7/4 and /S = k'^/3 - 1/2. With the change of variable ;z = 1/2 - zt/(2ao), 
this equation takes the form 

dz 
where 

= _ + "frtr (286) 

^ ' 4^2 (^ - 1)2 4^2 _ ]^)2 

After a direct calculation, Eqn. (12851) can be transformed into a hypergeometric equation 

;z(l - z) ^ + [c - (a + 6 + l)z]^ - afecu = 0, (287) 

where 



2+/(z)^7 = 0, (285) 



1 + (288) 



6=l-yl-a, (289) 

c = ^. (290) 

The solution for gi^z) is given by 

, . c , , c— a— fa— 1 , . , , 

g{z) = z^{z-\) — ^ — uj{z) (291) 
ot, in terms of the hypergeometric function F(a, 6, c; z), 



(7(2;) = ^! (^ - 1)^ i F + - a, i - ^ ^ - a, ^; ^ I . (292) 
Finally, the solution for the electric part of the Weyl tensor, is given by 



E, = .(-4a^)-t (. - l)-i F ( i + ^2 - ^, i - y/2 - ^, I . I (293) 



^■^Notice that, as shown in the Appendix, this equation is actuaUy a consequence of a transformation 
that takes the variables {E, S) (which are not canonicaUy conjugated) into a new pair of variables that are 
canonically conjugated. 
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where s is a constant. Restricting to z G 3?, it follows that this solution is regular for z < 1, 
and can be analitically extended for all values of z. Hence, the perturbation is regular. 

Notice that the power spectrum of the perturbations can be obtained using (see Ap- 
pendix) 

Pk = k-'\Ek\^. (294) 

The second example we shall study in this section is the model presented in Sect J4.5.T| 
the perturbation of wihch was analyzed by the quasi-Maxwellian method in |347j . In this 
model, the singularity is avoided by the introduction of nonlinear corrections to Maxwell 
electrodynamics, given by 

L = ~F + aF^ + (295) 

where F = F^^F^'^, G = \rjaj3^uF°^^ F^^" , a and [3 are arbitrary constants. After an average 
procedure (see Sect J4.5.Tl) . the expression for the scale factor for the "magnetic unverse" with 
e = is: 



a{tf = JTo 



^(t2 + l2a) 



nl/2 



(296) 



The interpretation of the source as a one-component perfect fluid in an adiabatic regime 
leads to instabilities [388] . which are artificial, as will be seen next. The sound velocity of 
the fluid in this case is given by |216j 

^ = ^ = ^ (297) 

dp p 9ip + p)- ^ ' 

This expression, involving only background quantities, is not defined at the points where 
the energy density attains an extremum given by 6* = and p + p = 0. In terms of the 
cosmological time, these points are determined by t = and t = ±tc = 12a. Notice that 
they are well-behaved regular points of the geometry, indicating that the occurrence of a 
singularity is in fact caused by an inappropriate description of the source. This difficulty can 
be circumvented by splitting the part coming from Maxwell's dynamics from the additional 
non-linear a— dependent term in the Lagrangian. As a result, we get two noninteracting 
perfect fluids: 

T^, = T«+T^2)^ (298) 

where 

T!,'J = {pi+Pi)v^v,-pig^,, (299) 
T^u) = ip2 + P2) v^v^ - p2 gf,u, (300) 
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and 



Pi = 


^ ^2 
/ 


Pi = 


1 -^2 

6 


P2 = 


-4a 


P2 = 


3 



(301) 

(302) 
(303) 
(304) 

From this decomposition it follows that each of the components of the fluid satisfies the 
conservation equation, thus showing that the source can be described by two non-interacting 
perfect fluids with equation of states pi = 1/3 pi and p2 = 5/3 p2- This splitting should be 
understood only as a mathematical device to allow for a fluid description. 
From the considerations presented in Sect Jll.5^ we obtain [347j: 

2Ai(36 + P) 
+ X2)P2 

Ei + ^^Ei = -^(l + Ai)piSi, (307) 

E2 + = (1 + A2) P2S2, (308) 

where k is the wave number. As shown in |170j . the scalar perturbations can be expressed 
in terms of the two basic variables Ei and Ej, and the corresponding equations can be 
decoupled. The result in terms the E^ is 

E, , '-l^heE, . (f . A.) ,i(l + - (2£±^!)^} E. ^ 0. (309) 

Note that in this expression there is no summation in the indices, and j ^ i ■, and Aj = (|, |) • 
In the first case the equation for the variable Ei becomes 



+ 1 KE2, (306) 



E, + ^OE^ + 



-e'^ - Pi - P2 - — 

3 3a^ 



El = 0. (310) 



Let us analyze the behavior of the perturbations in the neighborhood of the points where the 
energy density attains an extremum [i.e. the bounce and the point in which p + p vanishes). 
The expansion of the equation of Ei in the neighborhood of the bounce (at t = 0) up to 
second order, is given by: 

Ei + atEi + {b + bit'^)Ei = 0, (311) 
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where the constants a and b are defined as follows 



-4 (312) 

^~~vKmc ^^^^^ 

bi = -X- A- (314) 



Defining a new function / as 



fit) = E,{t)e^^{\+\-'-^h-^\e), (315) 



and introducing the coordinate ^ by 



a 



2 



i = -ie^^w~^, (316) 

we obtain for / the confiuent hypergeometric equation |389] 

e/ + (l/2-O/ + e/ = 0, (317) 

where 



lib -a/2) 1 
4(6i-aV4)V2 ~ 2' 

The solution of this equation is given by 



(318) 



f{t) =AM\d, 1/2, -it^s^h - ^ j , (319) 

where A is an arbitrary constant and M{d^ 1/2,0 is the confiuent hypergeometric function, 
which is well-behaved in the neighborhood of the bounce. Hence the perturbation E\{t) is 
regular and given by 



Ei{t) = AM {d,l/2,-it'^\lbi-^ 



X e.p{{-^^ + '-^b,---]t^}. (320) 

After a similar procedure, the perturbation E2 obeys, in the same neighborhood, the following 
equation: 

E2 + atE2 + ib + bif)E2 = 0. (321) 
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This is the same equation we obtained for Ei, with different values of a, 6 and bi given in 
this case by 

a (322) 
3 

^^n~^7^mc ^^^^^ 

^ ^ _ A. (324) 

The solution is given by the real part of 



E2{t) = AM [d,l/2,-it\lbi-^ 



ail 0? 



X ^^v{-\\-"-^\\h^-^-\t^). (325) 

so the perturbation -E'2(t) is well-behaved. At the neighborhood of the other critical point, 
given by t = the equation for the perturbation Ex is given by 

Ex + aEx + (6 + 6it) Ex = 0, (326) 

with 

a=A. (327) 

By the following variable transformation: 

Ex{t)^exp(^-^w{t)y (330) 

the differential equation goes to 

w+ {b- {a/2f + bxt)w^0, (331) 

and the solution is 

A b-(a/2f + bxt\ 
w{t) = Wo Ai ^-^^ '-j . (332) 

The Airy function Ai is regular near t = tc, and so is Ex- Finally we look for the equation 
of E2 at the neighborhood of t = tc- 

E2 + aE2 + {b + bxt) E2 = 0, (333) 
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where 



9 



(334) 




b 




) 

:) 



(335) 



b 



5^3 / 1 m2 



(336) 



1 



2t2 \t, 6^0 



This equation differs from Eq. fl326p only by the numerical values of the parameters a, b, and 
bi so we obtain the same type of regular solution 



Hence, it was shown by a direct analysis of a specific nonsingular universe, that in the neigh- 
borhood of the special points in which a change of regime occurs, all independent perturbed 
quantities are well-behaved, and the model is stable with regard to scalar perturbations. 

A similar analysis has been carried out for the model described by Eqn. (12961) in the 
case of tensor perturbations in [401] . The result shows differences between gravitational 
waves generated near a singularity and those generated near the bounce. While in the first 
case the system exhibits a a node-focus transition in the {E, E) plane, independently of the 
perturbation wavelength A, in the bouncing model the trajectories may exhibit a focus-node- 
focus transition, or no transition at all, depending on the value of A. 



As mentioned in SectJTTl another approach to the description of perturbations in a bouncing 
universe uses the idea of matching a contracting with an expanding phase. The hope here 
again resides in the fact that some general features can be extracted from given examples, 
since the matching may be done in such a way as to avoid a very detailed specification of 
the high curvature phase. Inasmuch as the result depends on the matching conditions, this 
issue was the subject of a long debate |345j . We shall present next some examples of this 
technique. 

The case of a scalar field with an exponential potential (inspired in the string pre-big 
bang and the ekpyrotic model) was studied in [183] . A matching between a contracting, 
scalar field-dominated phase and an expanding, radiation-dominated phase (and also of the 
corresponding perturbations) was done using the Israel conditions [402] . It was assumed 
that the slice of spacetime in which high-energy physics takes control is very thin, and can 
be approximated by a spacelike surface, with a negative surface tension (to be specified by 
the underlying physics) required by the jump in the extrinsic curvature. Neglecting possible, 
but subdominant, anisotropic surface stresses 0, and depending on the chosen surface, it 

^''This restriction was lifted in [403] . 




(337) 



11.3 Matching 
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was found that a scale-invariant spectrum could be transferred from the contracting to the 
expanding phase. A similar model has been studied in |198j . where it was shown that the 
value p = 2/3 of the power law a(t) oc (— t)^ was adopted for the scale factor generates a 
scale-invariant spectrum of adiabatic curvature fluctuations in the collapsing phase. The 
chosen background corresponds to a contracting Universe dominated by cold matter with 
null pressure. As a result of the glueing, the spectrum is matched at the bounce to a scale- 
invariant spectrum during the expanding phase. This model was also shown to generate a 
scale- invariant spectrum of gravitational waves, as already realized in [200J. 

It is useful to assume that the physics of the bounce is encoded in the transfer matrix T, 
defined by 

[DA ^ /Tn T,2\ (D^ 
\S+) \T2i T22) \S- 

T gives the degree of mixing between the dominant (D) and sub-dominant (S) modes before 
and after the bounce for a fixed comoving wave number k. Several combinations are possible, 
such as one for which the spectrum is initially not scale invariant but is turned into it because 
of a nontrivial k dependence of the transition matrix. Due to the fact that the bounce lasts 
only a short time, it is conceivable that it does not exert any influence on the large scales 
that are of astrophysical interest today. This implies that T does not depend on k [183] . in 
such a way that a scale invariant pre-bounce spectrum is transmitted without change to the 
post-bounce phase. This hypothesis has been tested in [184J . It was shown by way of an 
example (a bouncing solution in general relativity, with positive curvature spatial section, 
with a scalar field as a source, by using an expansion of the bouncing scale factor around 
the e = 1 de Sitter-like bouncing solution) that T may depend on k, provided that the null 
energy condition (NEC) is very close to being violated at the bounce, hence affecting the 
large scale behaviour of the scalar perturbations (see however [196[ 1197] ). Note however that 
it was shown in |184j that the spectrum of gravitational waves is not affected by the bounce. 

The authors of |18Uj have obtained the most general form of the transfer matrix respecting 
local causality. In particular, they have shown that no local-causality-respecting matching 
condition can lead to a scale invariant spectrum for both the pre-big-bang and the ekpyrotic 
model, in agreement with the result of |178j . They also studied a non-local model based on 
string theory and showed that under certain conditions a post-bounce SIS is possible. 

A different line of attack was pursued in [185] with the central assumption that the 
bounce in a spatially flat universe is governed by just one physical scale (chosen as r^^, 
the cosmological time at which the bounce occurs). Working in GR and incorporating all 
the eventual new physics in the matter side of EE, the general solution to the problem of 
the propagation of perturbations through the bounce was presented in |185] . It was shown 
that the spectrum of the Bardeen potential in the expansion phase depends critically on the 
relation between the comoving pressure perturbation and the Bardeen potential in the new 
physics sector of the energy-momentum tensor. Only if the comoving pressure perturbation 
is directly proportional to the Bardeen potential (rather than its Laplacian, as for any known 
form of ordinary matter), the pre-bounce growing mode of the Bardeen potential persists in 
the post-bounce constant mode. This would open the door to models with a scale-invariant 



(338) 
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spectrum (hence in agreement with observations) for those cases in which there is very slow 
contraction in the pre-bounce. This result is supported by numerical analysis of a toy model 
in which 6p (x |185] . Examples of this type of behaviour for the perturbations are given 
by models with spatial curvature (which cannot be treated however with this approach) and 
also by models with modifications coming from extra dimensions (such as the one presented 
in [US]) [T85] . 



11.4 Creation of cosmological magnetic fields 

The origin, evolution, and structure of large-scale magnetic fields are amongst the most 
important issues in astrophysics and cosmology. The standard model for the generation of 
this fields is the dynamo, which amplifies a small seed field to the current observed values of 
1 — few fiG. There are several mechanisms to produce these seeds, but the prevalent view 
is that they have a primordial origin |373j . In particular, the vacuum fluctuations of the 
EM field may be "stretched" by the evolution of the background geometry to super-horizon 
scales, and they could appear today as large-scale EM fields. For this to happen, conformal 
invariance of the EM equations must be broken. This is the case in models such as dilaton 
electrodynamics [374] and Weyl integrable spacetime (see Sect J3.2l and |375] for a list of 
references on the subject). 

As a previous step in the details of the case of the EM field, let us discuss the creation of 
massive scalar particles in a bouncing universe with e = — 1, following [371j . The expansion 
factor is given by a(t) = + Qq, or a{rj) = Oq cosh 77 in conformal time, as in the examples 
studied in [HI HZTj. The EOM for the scalar field is 



With the mode decomposition 



□0 + ( + -^R ]4> = 0. 



0fc(x) = a{T]) ^/Vk(f)xfc(?7), 



where k = {k, J, M) and the Y]^{x) are given in terms of the spherical harmonics (see [372j ). 
the function Xk{v) satisfies the modified Mathieu equation: 

where A = — (fc^ + ^rri^al), and h = |mao. The number of created quanta in the (asymptot- 
ically flat) future can be calculated with the solutions of this equation that have the right 
asymptotic behaviour, and following standard techniques. In the limit h << 1 (i.e when 
the Compton wavelength of the particle is much greater than oq), the result is [371] 



2 smh TTK \ ^ / 
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where k is the index in the Mathieu functions M_^j.{ri^ h), and is a comphcated function of 
A and h, which in the hmit for small h reduces to 

A = -P- J"' +0(/i'), 
2(A;+1) ^ ^' 

and if is a. phase, independent of h. The expression for varies from to 4x exp{—2nk) 
for large k, and shows that for a given k, the particle number depends on the product mao- 
The creation of magnetic fields in a bouncing universe in models that break the conformal 
invariance with a coupling to a scalar field was studied in |327l 14041 I375j . In the latter, 
canonical quantization was applied to the model given by 



where u is the scalar field, and F^p an abelian field, with /(cu) = exp(— 2a;). The modes of 
the potential = e'"^ satisfy the equation 

<:\v) + {k'-V{ii))Atl (339) 

where a = +, — designates the base of travelling waves, a = 1,2 describes the two transverse 
degrees of freedom, and V{ri) = —oj" + u'"^. For the background described in Eqn. flllSp . 

where a = v^/A, where A^ is the coupling constant of the scalar field to gravity. The mode 
equation fl339p admits analytical solutions in terms of hypergeometric functions, in terms of 
which the Bogolubov coefficients, and the expression for the energy density of the magnetic 
field pm can be calculated |375j . The amplification factor with respect to the conformal 
vacuum peaks for the modes with momenta such that k ~ 1.31, and is given by 

oc exp I ^ I , (341) 



(Pm)c/ 

for ?7 >> 1. The conditions for the spectrum to be greatly amplified today are [375] 

ao << ctr, A << 1, 

where tr is the time at which the scalar field is negligible, in such a way that the EM field 
is free again. 

At a comoving scale of about 10 kpc, the strength of conformal vacuum fiuctuations 
is of the order of 10~^^ G. To reach the strength required to feed the galactic dynamo, 
Bseed oc 10^^° G, , which is a conservative estimate, we get from Eqn. (13411) that A ~ 0.1. 
Taking for the comoving scale the size of the universe (^ 4 x 10^ kpc), the amplification 
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Figure 16: Plot of the potential (see Eq. fl340D ). for A = 0.1,0.3,0.5 (solid, dotted, dashed 
line respectively). 



factor becames 10^^, and we need A ~ 0.07. So the strength needed in both cases can be 
achieved by a modest value of A, the coupling constant of uj to gravity. 

These results were obtained in a model that did not take into account the effect of the 
creation of matter by the decay of the scalar field. The solution presented in Sect. (13. 2. 9p . 
namely 

aij]) = f3^/cosh{2r]) + ko sinh(2r/) - 2A:o(tanhr/ + 1), (342) 

with P = ao/y/l — ko, and < ko < 1/7. incorporates this feature, and its influence on the 
creation of photons was discussed in |327j . The result, displayed in Fig. (111. 41) . shows that 
there is a substantial increment in the number of photons if we take into account the effect 
of matter creation. 

11.5 Appendix 

In this appendix, we give a short summary of two gauge-invariant methods that can be 
applied to study the perturbations in cosmological scenarios. 

11.5.1 Perturbations using Bardeen variables 

The fluctuations of the metric tensor can be classified by their properties under spatial 
rotations into scalar, vector and tensor perturbations. In the linear theory, their evolution is 
decoupled. In the case of scalar perturbations, the perturbed metric of a homogeneous and 
isotropic spacetime can be written as 

ds^ = a\r]) {(1 + 2(p)dr]^ - 2B.idr] dx' - [(1 - 2ij)-fij + 2E^i.jdx'dx^]} , (343) 



117 




Figure 17: Plot of the mean number of photons as a function of the conformal time for 
m = 20, in the case without matter (dashed hne) and for the case with matter creation (full 
line), for A = 1. 



where is the metric of the 3-space. We shall sketch the case of hydrodinamical perturba- 
tions of a perfect fluid with energy-momentum tensor 

T- = (p + p)n"M^-p5"^. (344) 

Following [349] . it is convenient to build, from the four variables appearing in (13431) . two 
gauge-invariant quantities, given by 

^ ^ [{B-E')a'y ^ , a\B-E') 

a a 

In terms of these, the gauge-invariant perturbed EE are 

- m{n^ + ^') + V^* + 3A;^ = ^a^^To^^^)", (345) 

{n^ + ^!')^i = \^6T^^'''^\ (346) 

[(27^' + 7i:2<l> + W + ^" + 2W-A;^ + iv2(<l>-^)]5;--7*^($-^)|fc,- = --a^5Tf^\ (347) 

2 2 2 ■' 

^^For other cases, such as s scalar field, see [174] . 
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where the STp^^^"' are gauge invariant combinations of the STj^, B, and E (see |174j for 
details). 

In the case of hydrodynamical matter, the most general form of the perturbation can be 
written in terms of the perturbed energy Sp, the perturbed pressure 6p, the potential V of 
the 3- velocity v\t,x), and the anisotropic stress a as follows [349] : 



h -{po+Po)a ^V, 

(Po + Po)aVi -5p5ij + aij 



For the case of a perfect fluid, with energy-momentum tensor given by Eqn. fl344p . a^j = 0. 
The pressure perturbation can be split into its adiabatic and entropy components as 

Sp = 5p + 5S = cl5p + t5S. (348) 



dp J g \dS 

Entropy perturbations may be important in the case of two-component systems, such as 
plasma and radiation. 

The gauge-invariant perturbations of the energy-momentum tensor can be expressed in 
terms of the gauge-invariant energy density, pressure, and velocity perturbation: 

5To(^^)° = 5p^^'\ 5T/^^)° = (po + Po)a-'5u\''\ 5Tf^' = -Sp^^'^d), 

with 

5p^9i) = 5p + p'^(^B - E'), 6p^''^ = 6p + p'oiB - E'), 6uf^ = dm + a{B - E% 

From Eqns. fl345l) - (13471) applied to this case, it follows that $ = ^. Using Eqn.(l3lH]), the 

system can be written as 

$" + 3H(1 + c2)$' - c^V^* + [2n' + (1 + 3cl){n^ - A;)]$ = -a^rSS. (349) 

{a^)r = ^a\po+Po)5ufK (350) 

For adiabatic perturbations, Eqn. fl349p yields $, which determines Sp^^^^ through Eqn. fl345p . 
and 6ul^^^ through Eqn. fl350p . 

Eqn. fl349p can be simplified with the change of variables 



/I W-H' + k 
^ = V2V 

yielding 

a 



with 

1 / \ 1/2 / 

1 / On \ ' I 



1 / \ 1/2 / Q \ 1/2 



« VPo +PoJ V ct^Po 
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Vector perturbations 

The most general perturbed metric including only vector perturbations is given by 

-S' 



where the vectors S and F are divergenceless. From their transformation properties, it can 
be shown that 

a' = S' + F' 

(where the dot means derivative w.r.t. conformal time) is a gauge invariant quantity. For 
the perturbations of the stress-energy tensor, we have 

where and tt* are divergenceless. is related to the perturbation of the 4-velocity bu 

= ( v'/a 

Adopting the Newtonian gauge (in which F = 0), from the perturbed EE we get 

-l^AS^ = {p + p)V\ (351) 

- ^V,(a2(5^, + S\^)) = p(vrV + vr^,), (352) 
where A is the spatial Laplacian. From Eqn. fl35ip we get 

for the Fourier modes of V and S. Assuming that there are no anisotropic stresses, as in the 
case of presureless dust, we get from Eqn.( 



Vt{a'S\) = 0. 

Hence S\ = C k/ci^, where C is a constant. From this and Eqn. fl353p . we get 

« (354) 



^The results quoted in this section are taken from [ 174j 
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Note that Vj^ increases for uj = 0, and stays constant for radiation, but S\ always increases 
for decreasing a. 

Tensor perturbations 

These perturbations are built using a symmetric 3-tensor hij which satisfies the con- 
straints 

h\ = /i'^ = 0, 

in such way that the metric for tensor perturbations is 



From the perturbed EE we find (see for instance [391] 

where 5T^^^^^ is the gauge-invariant "pure tensor" part of 5T^j_y. In Fourier space, and 
introducing the rescaled variable hij = eijv/a, we have 

v'^ + (k^ - ^) = 0- (355) 



11.5.2 The quasi-Maxwellian method 

The QM method it has its roots in the formulation of Jordan and his collaborators [190J 
and uses the Bianchi identities to propagate initial conditions. The basic idea is to identify 
gauge invariant quantities from the beginning, using Stewart's lemma [191] , which guarantees 
that the perturbation of an object Q is gauge-invariant if Q is either constant or a linear 
combination of with constant coefficients. In conformally flat models, the Weyl tensor 
(deflned below) is identically zero, so its perturbation is a true perturbation, and not a 
gauge artifact. We shall see below how to obtain a minimum set of variables to completely 
characterize a perturbation, along with their evolution equations. 

Definitions and notation 

The Weyl conformal tensor is deflned by means of the expression 

aP/iui 

where 

ga|3^J.u = ganQpu — daugPfi, (356) 

and 

'2Ma/3fMu = Ra^lg^3u + R^3uga^l — RauQPfi — Rpfidau- (^57) 
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The 10 independent components of the Weyl tensor can be separated in the electric and 
magnetic parts, defined (in analogy with the electromagnetic field) as: 



Eap = -W^^.p,v^v\ (358) 

Hap = -W^^pXv". (359) 

The dual operation was performed with the completely skew-symmetric Levi-Civita tensor 
Va/Bfiu- From the symmetry properties of the Weyl tensor it follows that the operation of 
taking the dual is independent on the pair in which it is applied. 

It follows from these definitions that the tensors Efj^^ and H^i, are symmetric, traceless 
and belong to the three-dimensional space orthogonal to the observer with 4- velocity v^, 
that is: 

E^^^E,^, E^^v^^^Q, E^^gf^'^Q, (360) 

and similar relations for 11^,^. The metric gf^^, and the vector (tangent to a timelike 
congruence of curves F) induce a projector tensor hf^^, which separates any tensor in terms of 
quantities defined along F plus quantities defined on the 3-dimensional space orthogonal to 
v^. The tensor /i^,^, defined on this 3-dimensional space is symmetric and a true projector, 
that is 

V^"^ = - = K^. (361) 
We shall work with the FLRW geometry written in the standard Gaussian coordinate system: 

(is2 = + gijdx'dx^ (362) 

where g^j = —a^{tYfij{x^). The 3-dimensional geometry has constant curvature and thus the 
corresponding Riemannian tensor ^^^Rijki can be written as 

The covariant derivative in the 4-dimensional space-time will be denoted by the symbol 
and the 3-dimensional derivative will be denoted by "||". 

The irreducible components of the covariant derivative of are given in terms of the 
expansion scalar (9), shear (aap), vorticity (cUi^,^) and acceleration (Aa) by the standard 
definition: ^ 

Va;p = Cr^p + -^OKjB + l^ajB + (363) 

where 

= V^-a, 

(364) 



2"'[a"'/3]'"A';''' 
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We also define ^ 

Oap = cTq/? + -OKp. (365) 



Quasi-Maxwellian equations of gravity and their perturbation 

We shall present in this subsection a sketch of the deduction of the equations that govern 
the perturbations in the quasi-Maxwellian formalism. The details of the calculations in this 
section can be found in [170] . Using Einstein's equations and the definition of Weyl tensor, 
Bianchi identities can be written in an equivalent form as 

''^ 2 12^ 
2 6^ 

The quasi-Maxwellian equations of gravity are obtained by projecting these equations {i.e. 
the Bianchi identities are taken as true dynamical equations which describe the propagation 
of gravitational disturbances). The evolution equation for the perturbations for 59., 6(7^,^, and 
Su^, as well as 3 constraint equations, are obtained projecting and perturbing the equation 

which follows from the definition of the curvature tensor. Finally we get two more equations 
by projecting the conservation law T^'^.^ = 0. Adding up, we have a set of twelve equations 
which when perturbed yield (after straighforward manipulations) the coupled differential 
equations needed to give a complete description of the perturbation. In a general case, the 
variables are 

V = {6Eij, 6Hij, 6ujij, 6aij, diTij, 6ai, 6qi, 6p, 66, 6Vo, 6Vk} . 

From now on we will concentrate on the case of scalar irrotational perturbations. As shown 
in [166] . it is useful to develop the perturbed quantities in the spherical harmonics basis. It 
is enough for our purposes to work only with the scalar Q^'^^x^) (with dQ^''^ /dt = 0) and the 
vector and tensor quantities that follow from it, defined by Qf'^ = Q[i \ Q[f = Q^a-j- The 
scalar Q^''^ obeys the eigenvalue equation defined in the 3-dimensional background space by: 

v2gW =A;gW, (366) 

where k is the wave number, and the symbol denotes the 3-dimensional Laplacian: 

V^Q = Y'QmIj = l''Q,^■,J■ (367) 

Since the modes do not mix at the linear order, we will drop the superindex (k) from Q. 
The traceless operator Qij is defined as 

Qij = Qij + —Qlij, (368) 
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and the divergence of Qij is given by 

Q% = -2[e + ^^ Q\ (369) 

Due to Stewart's lemma, the good (since they are gauge-invariant an null in the background) 
objects in the list V are 5Eij, 5Ejj, SiTij, Sai, and 6qi. According to causal thermodynamics 
the evolution equation of the anisotropic pressure is related to the shear through [175 J 

ril,, + Uij = CcTij (370) 

in which r is the relaxation parameter and ^ is the viscosity parameter. For simplicity we 
will take the case in which r can be neglected and ^ is a constant Eg. (15701) then gives 

= ^aij, (371) 

and the associated perturbed equation is 

^n,, = e Saij. (372) 
We shall decompose the four independent and gauge-invariant perturbations as |^ 

k 

5s,, = 5^ sW(t)g;f, 

k 
m 
m 

It can be shown that ip is a function of S and E |170] . It follows that, restricting to the 
case g = (no energy flux) [^, E{t) and S(t) constitute the fundamental pair of variables in 
terms of which the dynamics for the perturbed FLRW geometry is completely characterized. 
Indeed, the evolution equations for these two quantities (which follow from Einstein's equa- 
tions) generate a dynamical system involving only E and S (and background quantities) 



^^In the general case ^ and r are functions of the equihbrium variables, for instance p and the temperature 
T and, since both variations 511^^ and daij are expanded in terms of the traceless tensor Qij., it follows 
that the above relation does not restrain the kind of fluid we are examining. However, if we consider ^ as 
time-dependent, the quantity 511^ must be included in the fundamental set M[a]- 

^^In fact, SEijSE"^^ is the only quantity that characterizes without ambiguity a true perturbation of the 
Debever invariants |170j . 

^^We further assume an equation of state relating the pressure and the energy density, i.e. p ~ Xp, which 
is preserved under arbitrary perturbations. 
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which, when solved, contains all the necessary information for a complete description of all 
remaining perturbed quantities of the FLRW geometry. 
The evolution equations are given by |17U] 

t = -E-]^^T.-k^i), (373) 



(374) 



As mentioned before, if) can be expressed in terms of E and S 
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(1 + A)p^ = 2(^l + ^^ a-2[_Ai? + iAeS + i^S]. (375) 

Thus the set of perturbed equations reduces to a time- dependent dynamical system in the 
variables E and S: 

S = Fi(S,E), E = F2(S,E), (376) 

with 

Fi = - S - (377) 

and 

\e + ^^P+\^o)^ (378) 
where \& is given in terms of E and S by eq. (13751) . so the system [376] can be written as 



E\ _ fa (3\ (E 
a) ~ 

where 
Since 

as _ 
as ~ ~3' 



(379) 
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Except when (1 + A) = 0, see |170) for this case. 
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the system 13791 is not Hamiltonian due to the expansion of the universe. Nonetheless, new 
variables (Q, P) can be introduced in such a way that the system 13791 is Hamiltonian in 
them. Defining 

Q = a'^a, P = a^'E, 

it is easily shown that the otherwise arbitrary powers m and n must satisfy the relation 
m + n = 1 for the variables Q and P to be canonically conjugated. It follows that 

P = MiP + M2Q. 

The choice n = 3A/2 + 2 yields = 0, and P satisfies the equation 

P + fi{t)P = 0, 

with 

i^W = (^^ + V\p + -^ \ 2 \ 2 



5, 2\ 1 



which is equivalent to Eqn. f l282p . 

This method can be extended to vector and tensor perturbations in the FLRW model 
|171j . In the first case, the observable quantities are described in terms of the vorticity 
and the shear, while the electric and magnetic parts of the Weyl tensor suffice for the 



gravitational waves r^H . The three types of perturbation are describable in Hamiltonian 
form, thus paving the way to canonical quantization [172] . which was performed for scalar, 
vectorial, and tensor perturbations using the squeezed states formalism in [172] . In fact, in 
the case of scalar perturbations, the Hamiltonian in terms of the {Q, P) variables (with the 
choice m = is given by 

H=^Q' + ^P^ + ^P, 
2^2 Q 

with 

hi = — , h2 = —— e + ^ a, /13 = 0. 

2a 1 + A\ 6 J ap 

11.6 Relation between the two methods 

The Bardeen variables ^P) are related to the quasi-Maxwellian variables (-E, S). For 
instance, in the case of scalar perturbations the relation between E and $ (for a perfect 
fiuid) is given by [1221 1123] 

E = -A;2$, 

from which the relation for the spectrum given in Eqn. (12941) follows. 



^'^^Perturbations in the Kasner solution were studied in [40 6| . 
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12 Conclusion 



The idea of a bouncing universe has been considered since the early days of relativistic 
cosmology, as shown in this review. However, only a few analytical solutions describing a 
nonsingular universe served as a starting point to build a complete cosmological scenario. 
The main reason for this neglect by the majority of the physics community in the last 30 years 
of the 20th century was the strong influence of the singularity theorems, which led to the 
belief that some sort of singularity was inevitable in gravitational processes. The situation 
should have changed with the recently discovered positive acceleration of the universe since, 
in the realm of GR, the accelerated expansion means that the matter content must satisfy 
the condition p + 3p < 0, which is precisely one of the conditions needed to have a bounce in 
Einstein's gravity. This violation of the SEC was already accepted in the early 80 's in order 
to have a phase of inflationary expansion, and nowadays several systems are known which 
do not satisfy the inequality p + 3p > (see for instance [26]). Hence, there is mounting 
evidence against one of the main theoretical prejudices forbidding bouncing universes in GR. 
Surprisingly, nonsingular models have not arisen the interest that should be expected based 
on the preceding considerations 
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Almost contemporaneous to the discovery of the accelerated expansion was the gradual 
advent of a handful of cosmological models based on nonsingular solutions. These models 
aimed at solving the most stringent problems of the (pre-inflationary) cosmological stan- 
dard model: the initial singularity, the isotropy and homogeneity of the currently-observed 



universe, the horizon problem, the flatness probleni I. and the formation of structure 
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Bouncing universes have partially met these challenges. The singularity is obviously absent, 
and its avoidance requires any of the assumptions listed in Sect Jl.ll which range from the 
violation of SEC (in GR) to quantum gravitational effects. 

As explained in SectiH a phase of accelerated contraction may solve the flatness problem 
in GR, and may also get rid of particle horizons (see for instance |91] ) 
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Finally, the amplification of primordial seeds (a problem prior to the formation of struc- 
ture) in bouncing universes has been intensely debated recently (see Sect JTT]) . The asymp- 
totic behavior of these universes is markedly different from that of the SCM or inflation. The 
universe at past infinity starts to collapse from a flat empty structure-less state, that at past 
infinity can be approximated by Minkowski geometry written in terms of Milne coordinates 



^•^^It may be argued that this lack of interest is due to the fact that the bounce is expected to involve 
scales where quantum effects render GR inapplicable. But this is true also of the singularity theorems, as 
was known already in the early 70's. Moreover, there is no evidence against the possibility of a bounce in 
the classical regime [7S], as follows from some of the models presented in Sect lll.11 see also |425) ). 

^•^•^Note that the flatness problem may in principle not be a problem in gravitational theories other than 
GR (see SectH^). 

^'^''in spite of its historical importance, the so-called monopole problem is not included in this list, since 
there is still room for it to be be considered as a problem of field theory first, and then (perhaps) of the 
standard cosmological model, see for instance |424 i 1426] . 

^"^See however the concerns in |429j about the efficiency of some bouncing models in erasing possible initial 
inhomogeneities. 
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The transmission of the quantum fluctuations from this initial state to the post-bounce 
phase is strongly model-dependent, but there are some models which yield a scale- invariant 
spectrum for the scalar perturbations in the post-bounce phase (see Sect Jll.lTTl) . 

An offspring of the bouncing models are the cyclic universes (see Sect JTOj) . The cyclic 
models also attempt to solve the above-mentioned problems, and also may offer a new view 
on the initial conditions: since by definition, there is neither a beginning nor an end of time 
in these models, there is no need to specify initial conditions. Generically, cyclic universes 
share the problems of the universes that bounce only once. In addition, they must assure 
that the large scale structure present in one cycle (generated by the quantum fluctuations 
in the preceding cycle) is not endangered by perturbations or structure generated in earlier 
cycles, and will not interfere with structure generated in later cycles. One of the latest cyclic 
models, presented in |43U] . claims to have successfully faced these issues (however see |431j ). 

As compelling a scenario may (or may not) seem, the ultimate judge is observation, so 
we can ask if there are any that may point to the occurrence of a bounce. As far as we know, 
there are two possibilities 
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• As discussed in Sect Jll.ll the tensor spectrum of a nonsingular universe has a unique 
feature. As an example, the SPBB models predict a stochastic spectrum of gravita- 
tional waves whose amplitude increases as a function of frequency in some frequency 
ranges (see Sect Jll.lT^ . hence avoiding the bounds due to the CMB, pulsar timing, 
and Doppler tracking [433] . The parameter space of the "minimal" SPBB model |432] 
was limited using LIGO results in [433]. Notice also that nonsingular universes may 
produce vector perturbations (see Sect Jll.l]) . 

• The bounce may cause oscillations, that will be superimposed on the power spectrum 
of scalar perturbations. These oscillations would also appear in the WMAP data, 
linked to the spectrum through the multipole moments which are in turn defined 
through the two-point correlation function of the temperature fluctuations. [390J . Let 
us note however, that such oscillations may be due not only to a bounce, but also to 
transplanckian effects |390j or to non-standard initial conditions in the framework of 
hybrid inflation [427] . 

We would like to close by pointing out that although they do not yet give a complete 
description of the universe, a better understanding of bouncing models in classical GR should 
be attempted since they are inevitably imposed upon us by the apparently observed violation 
of the strong energy condition. It must also be noted that there are at least two more reasons 
to attempt this task. First, the current solution to the problems of the standard cosmological 
models (namely inflation) is successful, but has several problems (see Sectil]). Second, even 
if bouncing models do not succeed in yielding a complete description of the universe (thus 



^"^We have also seen that there are eternal (non-bouncing) universes, that start in a de Sitter regime. 
-'^"^Some bouncing models in GR were severely restricted in [428] . using SNIa data, CMB analysis, nucle- 



osynthesis, and the age of the oldest high-redshift objects. 
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offering an alternative to inffation they may throw hght upon the singularity problem 
(an issue in which inffation has nothing to say). 

Summing up, we have seen in this review that bouncing universes have some attractive 
features, but they are not complete yet: much work is needed to achieve a stage in which 
their predictions can match those of the cosmological standard model. Therefore, we hope 
this review encourages further developments in nonsingular cosmologies. 
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